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USE OF TETRAFLUOROMETHANE TO SIMULATE REAL-GAS EFFECTS 
ON THE HYPERSONIC AERODYNAMICS OF BLUNT VEHICLES 

By Robert A. Jones and James L. Hunt 
Langley Research Center 

SUMMARY 

An experimental and theoretical study was made of the use of tetrafluoromethane 
(CF 4 ) in low hypersonic flows as a means of simulating the flow field over the forebody 
of blunt configurations in hypersonic flight where dissociation of the gas in the shock 
layer gives large density ratios across the shockwave. A pilot model hypersonic CF 4 
blowdown tunnel at the Langley Research Center was built and studied. With this facility 
a normal- shock density ratio of 12.1 was obtained at a Mach number of 6 with stagnation 
conditions of 1650 N/cm2 and 736° K. The results of the investigation indicate that flow 
fields over the forebody of blunt configurations for real-gas conditions resulting in 
normal-shock density ratios up to 17 can be simulated by this method. A computer pro- 
gram for isentropic flow of CF 4 as well as flow across normal and oblique shocks was 
set up. Plots of various flow parameters are presented as a function of free-stream 
Mach number for various stagnation temperatures. A method for computing axisymmet- 
ric nozzle contours which includes a correction for boundary-layer displacement thick- 
ness is described. 


INTRODUCTION 

At supersonic and hypersonic speeds the aerodynamic characteristics, including 
the pressure distribution, shock shape, drag, and stability, of very blunt, low-fineness- 
ratio bodies are determined almost exclusively by the flow field over the forebody. The 
inviscid forebody flow in turn is essentially independent of Mach number and depends 
mainly on the density ratio across the normal shock and the body geometry. (See, for 
example, the discussion of Mach number independence principle in ref. 1 and the discus- 
sion of real-gas effects on the flow field of several high-drag planetary entry bodies in 
ref. 2.) Therefore, the inviscid forebody flow can be simulated in a wind tunnel if the 
normal-shock density ratio of flight can be matched. 

For hypersonic flight of large vehicles in planetary atmospheres with equilibrium 
gas conditions, density ratios as high as 20 are encountered. These large density ratios 
result from chemical dissociation of the molecules in the high-temperature gas behind 



the shock. However, the maximum normal- shock density ratio that can be achieved in 
air wind tunnels, even in those with large specific -energy content, is usually less than 8 
because the small size of the models precludes the attainment of chemical equilibrium 
in the forebody shock layer. Thus, a means of experimentally simulating high shock 
density ratios is needed. 

It has been proposed (refs. 3 and 4) that flow-field calculations for blunt bodies 
where real-gas effects are significant can be greatly simplified by the use of ideal-gas 
relations but with a low value of the specific-heat ratio to obtain the correct shock den- 
sity ratio. It also appears possible that the desired density ratios may be obtained 
experimentally with no dissociation by testing in a relatively low- temperature flow of a 
gas that has a low ratio of specific heats. If the Mach number is sufficiently high to 
apply the Mach number independence principle, then the desired inviscid forebody flow 
field can be simulated. If the Reynolds number is matched also, then the entire forebody 
flow field can be simulated. 

The use of gases other than air in aerodynamic facilities has been considered by 
many previous investigators. For example, Chapman (ref. 5) considered several gases 
and gas mixtures including tetrafluoromethane (CF 4 ) for various applications. He 
showed that for a given subsonic Mach number, Reynolds number, and pressure, a gas- 
mixture wind tunnel having the same specific-heat ratio as air could be smaller and thus 
require much less power. Also it has been shown that by using gases other than air, it 
is possible to match free- stream Mach number at much lower stagnation temperatures. 
The helium tunnels now in operation are examples of this application. 'The use of CF 4 
has also been considered to obtain high Mach numbers at lower stagnation temperatures 
(ref. 6 ) but was found to be somewhat impractical for that purpose. The only previous 
attempt known to the authors to simulate the real-gas effects on flows over very blunt 
bodies at high speeds by the use of a gas with a low ratio of specific heats is that of ref- 
erence 7 in which a mixture of water vapor and oxygen is proposed. 

The purpose of this paper is to (1) discuss the normal-shock density ratio as a 
simulation parameter for blunt-body flows, (2) propose the use of CF 4 as a test gas in 
order to obtain high shock density ratios, (3) present some results of an experimental 
investigation in both a wind tunnel and a ballistic range using CF 4 as the test gas, and 
(4) present equations and charts for the compressible flow characteristics of CF*. 

A small (7.6-cm-diameter test section) CF 4 tunnel with a conical nozzle was built 
and put into operation at the Langley Research Center in October 1965. Later, a 
ballistic- range investigation at the Naval Ordnance Laboratory (NOL) using CF 4 as the 
test gas was supported by NASA (Langley Research Center) to study effects of high den- 
sity ratios on the aerodynamics of three very blunt low-fineness-ratio configurations. 
Some of the results from both these investigations are presented herein. Also, the 
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design characteristics of a facility which would be inexpensive to build and simple to 
operate but would enable the correct simulation of real-gas density- ratio conditions for 
blunt configurations are presented. A few of the results (obtained in the small CF4 wind 
tunnel) that show the sensitivity of the flow field over a tension- shell configuration to high 
shock density ratios have been published in reference 8. Some of the results obtained m 
the NOL ballistic -range investigation are given in reference 9 . 

This report includes appendixes by James L. Hunt, Kathryn A. Smith, Robert B. 
Reynolds, and Lillian R. Boney, of the Langley Research Center, which present calcula- 
tion procedures for evaluating isentropic expansions and flow across shocks in CF4, a 
method for determining hypersonic nozzle contours for CF4, and related computer 
programs. 


SYMBOLS 


A 

A* 

C D ' 

H 

M 

P 




s 

T 

u 


x,y 


inviscid test- section area 

area at sonic throat 

drag coefficient 

enthalpy 

Mach number 

pressure 

radius 

f ree _ s tream Reynolds number (based on base diameter) 

surface distance from stagnation point 

temperature 

velocity 

coordinates 


y 


ratio of specific heats 



A shock standoff distance 

e inverse density ratio, 

6 shock-wave inclination angle 

p viscosity 

P density 

Subscripts: 
b base 

eff effective 

max maximum 

n nose 

0 free-stream stagnation conditions 

1 free-stream static conditions 

2 static condition behind normal shock 

3 stagnation condition behind normal shock 

BLUNT-BODY FLOW- FIELD SIMULATION 
Ideal-Gas Flow Fields 

For high-speed flight in an atmosphere the chemical dissociation of the gas in the 
s hock layer can result in a marked increase in density ratio across the strong bow shock 
compared with a lower speed or lower temperature flow of the same gas where no dis- 
sociation takes place. Under conditions where dissociation is present, the aerodynamic 
characteristics of very blunt vehicles (shock-wave angles larger than 50°) depend pri- 
marily on this shock density ratio (ref. 1) and not on the Mach number if Mj > 4. Three 
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blunt configurations which will be used to illustrate this strong dependence of aerody 
namic characteristics on shock density ratio are shapes being studied by NASA as pos- 
sible Mars probe/lander vehicles. These three configurations are shown m figure 1. 

In order to land a payload on the surface of Mars and make use of aerodynamic braking 
in the thin atmosphere, a large high-drag vehicle is required. Because of the large 
amount of C0 2 ir the Mars atmosphere and the large size of the vehicle, high shock den- 
sity ratios will be encountered even at velocities as low as 3000 m/sec. 

The changes in aerodynamic characteristics which are due to the large shock den- 
sity ratios associated with high-speed flight are primarily the result of changes in sur- 
face pressures acting on the forebody. The surface pressures are affected by a change 
in shock density ratio in two ways: First, the level of pressure at the stagnation point is 

changed, and second, the nondimensional distribution of surface pressure relative to 
stagnation-point pressure is changed. This dependence on shock density ratio can be 
illustrated by considering the flow of a perfect gas about a blunt body. The density-ra 10 
effect on the stagnation-point pressure level is easy to estimate. From conservation of 
mass and momentum across a normal shock, the following equation is obtained: 


P 2 ~ p l 

Pi u i 2 


= 1 - e 


( 1 ) 


which holds true for both ideal and real gases. For high Mach numbers p, can be 
neglected compared with p 2 , and the flow from the shock to the stagnation point can be 
considered incompressible so that equation (1) reduces to 

P 3 _ _ 1 (2) 

p l u l 2 2 

This effect on the vehicle drag coefficient of a given shape can be approximated as 

^ Cp \arge density ratio _ ( 2 ~ 6 Wge density ratio, ( 3 ) 

( C D)small density ratio (2 " e) small density ratio 

For a blunt vehicle in flight with a shock density ratio 1/e of 15, the forebody drag 
would be increased about 5 percent over that measured in a conventional high Mach num- 
ber air tunnel (where the maximum density ratio is 6) because of this increase in 
stagnation -point pressure level. 

The surface pressure distributions for the spherical segment, the 120° cone, and 
the 140° cone have been computed for ideal-gas flow using a one-strip integral method 
(ref. 10) at various shock density ratios and free-stream Mach numbers. These results 
are shown in figures 2 and 3. The pressure distribution for the spherical segment 


5 



(fig. 2) is affected by density ratio mainly at large values of s/r n where the effects on 
drag and moments are largest. Figure 2 shows that the effects of free-stream Mach 
number on pressure distribution are negligible for Mach numbers of 6 or larger. The 
results for the 120® cone (fig. 3(a)) also show a negligible effect of Mach number on sur- 
face pressure. For density ratios somewhat above 10 the sonic point moves forward 
from the base of the cone. Since the program requires the sonic -point location to be at 
the point of maximum radius, no solutions could be obtained for the 120° cone at a den- 
sity ratio of 20. Solutions were obtained for the 140° cone (fig. 3(b)) at a density ratio 
of 20. For this case the flow over the conical portion is almost sonic; this gives the 
sharp decrease in pressure near the sphere-cone juncture (s/r b * 0.08). The results 
shown in figures 2 and 3 indicate that for bodies as blunt as these, the primary factor 
governing the aerodynamic characteristics is the shock density ratio. Therefore this 
shock density ratio must be simulated when such shapes are studied experimentally. 

Other aerodynamic characteristics of blunt bodies which depend primarily on shock 
density ratio are the shock shape and shock standoff distance. The shock shapes of these 
blunt bodies calculated by the method of reference 10 are given in figures 4 and 5. Here 
again it is apparent that Mach number has very little effect. Both the shock shape and 
the standoff distance at the stagnation point are strongly dependent on the shock density 
ratio. 


Real-Gas Flow Fields 

The surface pressure distributions for real-gas equilibrium flow have been calcu- 
lated for all three shapes of figure 1 in reference 2. The results of these calculations 
are presented in figures 6, 7, and 8. The value p 2 / p 1 = 5.6 corresponds to an ideal 
gas at a Mach number of 9 with y = 1.4; the value of p 2 /pj = 16.7 corresponds to a 
condition near peak heating and peak dynamic pressure for entry into an assumed Mars 
atmosphere (model 2 of ref. 11). The blunt-cone pressure distribution of figure 7 was 
made for a 130° cone since that was the smallest cone angle for which the sonic point 
was located at the rearward corner for a shock density ratio of 16.7 and thus the small- 
est cone angle for which solutions could be obtained with the method used in reference 2. 
The trends shown by the data of figures 6 and 7 are in close agreement with the ideal-gas 
calculations shown in figures 2 and 3; this indicates that real-gas equilibrium flow and 
ideal-gas flow give about the same surface pressure distribution at high Mach numbers 
provided the value of the normal- shock density ratio is the same. Therefore, even when 
real-gas effects are present, the density ratio is still the dominant parameter as long as 
the bow shock is strong and the fineness ratio is low. 


Predicted Simulation in CF4 

In order to duplicate high density ratios obtained in flight at high hypersonic speeds 
in a wind tunnel without dissociation, the test gas must have a low ratio of specific heats. 
Tetrafluoromethane (CF4) was selected as the working fluid for this simulation because 
of its low ratio of specific heats, low boiling point, thermal stability, and low vibrational 
relaxation time, and also because it is a readily available, colorless, odorless, nontoxic 
gas that is easily reclaimed; however, no exhaustive search to find the best test gas for 
this particular type of simulation was made. Much of the physical and chemical property 
data for CF4 used in this study was obtained from the manufacturer of CF4. Some of the 
physical properties are given in table I. 

TABLE L- PHYSICAL PROPERTIES OF TETRAFLUOROMETHANE (CF 4 ) 


Molecular weight 88.01 

Boiling point at 1 atm*, °K 145 

Freezing point, °K 89.4 

Critical temperature, °K 227.33 

Critical pressure, atm 36.96 

Critical volume, cm 3 / mole . . . 141 

Critical density, kg/m 3 0.626 X 10 3 

Density, liquid at 193° K, kg/m 3 1.317 x 10 3 

Density, saturated vapors at boiling point, kg/m 2 7.62 

Specific heat, liquid at 193° K, J/kg-°K 1.23 x 10 2 

Specific heat, vapors at 1 atm and 298° K, j/kg-°K 0.707 x 10 2 

Specific-heat ratio at 1 atm and 298° K 1.159 

Heat of vaporization at boiling point, J/kg 136.0 x 10 2 

*1 atm = 101 kN/m 2 . 


The method of reference 12 was used to compute the equilibrium composition of 
CF4 as a function of temperature for several pressures covering the range that would 
normally be used in a wind tunnel, and it was found that no dissociation occurred for tem- 
peratures lower than 1600° K. A more complete description of the properties of CF4, 
including the equations of state, and the computer programs used to calculate conditions 
for isentropic expansion and flow across normal and oblique shocks in CF4 are given in 
appendixes A and C. Some of the flow variables of interest which were computed as 
described in appendix A have been plotted as functions of pertinent parameters for the 
ranges of stagnation chamber pressure and temperature studied. The stagnation enthalpy 
Hq, which is practically independent of stagnation pressure in the range from 1034 to 
1724 N/cm 2 , is given as a function of stagnation temperature Tq in figure 9. Plots of 
U l/|2H0, Yl> Tj/Tq, P ly /P 0 , Pi^/Mi, P\^\ 2 /^ P 2 /PV P 3 /P 0 > P 3 /Pi> y lj y V 
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and P2 u 2 / A i o a S ainst M 1 are P resentec * in figures 10 to 20 for Tq = 478°, 588°, 700°, 

and 811° K. The specific-heat ratio, unit Reynolds number, and Reynolds number behind 
normal shock are given for stagnation pressures of 1034, 1378, and 1724 N/cm^. The 
other flow variables are independent (within the accuracy of the plots) of stagnation pres- 
sure in the range from 1034 to 1724 N/cm2). Figure 16 shows that for a temperature of 
811° K, a shock density ratio of 13.6 is obtained at a Mach number of 7. The calcula- 
tions of appendix A show that for 1400® K, a shock density ratio of 14.5 is obtained. At 
this same temperature a value of 15.6 can be obtained if the Mach number is increased 
to 8 rather than 7. 

Similar calculations have been made by the method of appendix A for test condi- 
tions in a ballistic range. Plots of the density, temperature, and pressure ratios across 
a normal shock as a function of velocity in 100 percent CF 4 are shown in figures 21 to 23. 
The ambient temperature in the range is assumed to be room temperature (294° K) and 
the plots are accurate within 4 percent for range pressures from 30 to 1000 N/m2. Cal- 
culations are not shown for launch velocities greater than 1500 m/ sec because dissocia- 
tion of CF 4 in the shock layer could occur at these higher velocities. Figure 21 indicates 
that shock density ratios up to 17 can be obtained by using CF 4 as the test gas in a bal- 
listic range. 

The shock density ratio and Reynolds number simulation capability of the CF 4 
facility (assuming a 2.5-cm-diameter model and a realistic pressure range) is shown in 
figure 24, and that of several other facilities at the Langley Research Center, which are 
thought to be typical of most facility capabilities, is shown in figure 25. Also shown in 
both figures are two entry trajectories typical of very blunt high-drag vehicles. Both 
these trajectories are based on the assumption of chemical equilibrium in the shock layer. 
The Earth entry trajectory is for an entry angle of 90°, an entry velocity of 5600 m/ sec, 
a ballistic coefficient of 146 kg/ m^, and a diameter of 0.67 m. The Mars trajectory is 
for an assumed atmosphere of 2000 N/m^ (model 2, ref. 11), an entry angle of 90°, an 
entry velocity of 8700 m/sec, a ballistic coefficient of 23.5 kg/m^, and a diameter of 5 m. 
The "knee" of the trajectories corresponds roughly to conditions of peak heating and peak 
dynamic pressure. Figure 24 shows the simulation capability of a CF 4 wind tunnel 
operating in the Mach number range from 3 to 8 and that of a CF 4 ballistic range at 
velocities from 450 to 1500 m/sec. A comparison of figures 24 and 25 shows that with 
CF 4 the ability to simulate shock density ratios which occur with high-speed atmospheric 
entry is increased. Therefore the use of a relatively low velocity flow of CF 4 should 
greatly facilitate the study of real-gas effects on the aerodynamics of blunt entry con- 
figurations (those having shock-wave angles larger than 50°). 

The variation of shock density ratio with shock-wave inclination angle 6 for a 
CF 4 flow at a Mach number of 6 and a stagnation temperature of 811° K is shown in 
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figure 26. Also shown is the shock-density-ratio variation of an ideal gas at a Mach 
number of 6 with a value of y which results in the same normal- shock density ratio as 
CF 4 . These two curves are compared with one for equilibrium air at an altitude of 
70 km and a velocity of 5240 m/sec; for these conditions the free-stream Mach number is 
16.5 and the normal-shock density ratio is 12.7 (same value as for the other two curves). 
This plot shows that for shock waves with 0 greater than about 60°, the CF 4 simulation 
gives very nearly the same shock-density- ratio variation as would occur in flight in the 
earth's atmosphere for equilibrium conditions at a Mach number of 16.5 and an altitude 
of 70 km. It also demonstrates that for shock-wave angles greater than about 60°, the 
CF 4 simulation gives approximately the same result as an ideal gas at the same Mach 
number and with the value of specific-heat ratio adjusted to give the same value of the 
normal-shock density ratio P 2 /Py The CF 4 simulation method should therefore be very 
accurate for bodies blunt enough to have bow shock-wave angles of 60° or greater. At a 
shock-wave angle of 50°, the density ratio in CF 4 is 87.5 percent of that which occurs in 
equilibrium air at the stated conditions. This difference is believed to be within accept- 
able limits for engineering calculations; thus, blunt bodies have been heretofore defined 
in this report as bodies blunt enough to have bow shock waves of 50° or greater. 

PILOT CF 4 TUNNEL 
Apparatus 

A schematic diagram of the 7.6-cm CF 4 tunnel is given in figure 27. The CF 4 
supply consisted of bottles each containing about 32 kg at a pressure of 1440 N/cm2. The 
small bag-type pneumatic accumulator was used together with a check-valve arrange- 
ment to pump the CF 4 into the large accumulator at the pressure desired for a test. 

When the desired pressure was reached in the large accumulator, the high-pressure air 
regulator was set at this value and the valve between the large air bag and air supply was 
opened. As CF 4 was used during a test, the large air bag inflated and thus maintained a 
constant tunnel stagnation pressure for times as long as 3 minutes and eliminated any 
chance of contamination of the CF 4 . 

The heat exchanger consisted of several long coiled stainless- steel tubes immersed 
in a bath of molten lead. Energy was supplied by electrical strip heaters which were 
also immersed in the lead. The heat exchanger, the line from the heat exchanger to the 
stagnation chamber, and the stagnation chamber itself were all heated by electrical strip 
heaters and wrapped with insulation. The lead bath could be maintained at any tempera- 
ture up to 812° K but the line and stagnation chamber metal temperatures were limited 
to a maximum of 534° K. This system could supply CF 4 at stagnation conditions up to 
2000 N/cm^ and 750° K. 
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The nozzle was conical with a 5° half-angle, a 1.32-mm-diameter throat, and a 
7.6-cm-diameter test section. The tunnel exhausted into a 1200-m 3 vacuum sphere which 
could be pumped to pressures as low as 25 N/m 2 . No attempts were made to reclaim the 
CF4, although it should be possible to do so easily since it is a heavy gas with a relatively 
high critical- temperature (227.33° K) and a relatively low critical pressure (36.96 atm). 

Experimental Results and Discussion 

Several different measurements have been made in this pilot facility. The primary 
purpose of these measurements was to demonstrate that high normal-shock density ratios 
could be obtained with this facility and to verify the theoretical test conditions calculated 
by the method discussed in the appendixes. A secondary purpose was to study the oper- 
ating characteristics so that a more nearly optimum facility could be designed. Measure- 
ments of the stream pitot pressure, test- section wall static pressure, stagnation chamber 
pressure and temperature, shock standoff distances, and shock shapes for several differ- 
ent configurations were made for a range of total temperatures and pressures. 

Nominal test- section conditions . - Although tests were made for a wide range of 
stagnation temperatures and pressures, a preponderance of data were obtained at stagna- 
tion conditions of 736° K and 1650 N/cm 2 . These conditions are referred to subsequently 
as the nominal test conditions. The corresponding test-section conditions as determined 
from the measured stagnation pressure ratio across the normal shock and the measured 
ratio of free-stream static to free-stream stagnation pressure in conjunction with the 
real-gas calculation procedure of appendix A (figs. 11 to 19) are presented in table II. 

TABLE II.- NOMINAL TEST-SECTION CONDITIONS 


Condition 

Value calculated by method of appendix A for - 

p 3 /p 0 = 1.13 x 10-3 * 

Pl /p 0 = 2.25X10-5* 

Mj 

6.06 

6.20 

P2/ p l 

12.2 

12.3 



1.190 

1.205 

^2/ y l 

0.925 

0.900 


^Measured value at Pq = 1650 N/cm 2 and Tq = 736° K. 


Density ratio measurement .- Since the primary concern of this report is the simu- 
lation of forebody flow fields of blunt configurations by matching the normal- shock density 
ratio without dissociating the test gas, a more direct method for verifying this ratio was 
desired. 
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Normal- shock density ratios in the test section were determined from the shock 
standoff distance measured for a sphere by using the following relation: 

£=0.76e (4) 

This relation is an empirical fit to the theories of references 13, 14, and 15 and has been 
shown to predict very accurately the shock standoff distance for spheres. A test-section 
normal- shock density ratio of 12.1 was determined with this method for the nominal stag- 
nation test conditions (p Q = 1650 N/cm 2 and Tq = 736° k). This value compares favor- 
ably with those from the real-gas calculations (12.2 and 12.3) given in table II. These 
values are more than twice the normal-shock density ratio which can be obtained in an 
air wind tunnel operating as an ideal -gas facility. 

Effective test-section conditions .- As discussed previously in the section "Blunt - 
Body Flow-Field Simulation," it should be possible to calculate approximately the flow 
field about very blunt bodies using ideal -gas relations for some effective Mach number 
and specific-heat ratio provided that the correct normal-shock density ratio is obtained. 
One method for determining the values of effective Mach number and effective specific - 
heat ratio is to solve the following two ideal -gas equations simultaneously using the value 
of the normal-shock density ratio determined from the measured shock standoff distance 
of a sphere and the measured ratio of stagnation pressure behind the normal shock to the 
total pressure in the stagnation chamber. 


e ( y eff ~ 1 ) M eff + 2 
( y eff + 1 ) M eff 


y eff 


( y eff + >) M eff 

y eff ~ 1 

y eff + 1 

JL 

y eff _1 

( y eff ‘ 1 ) M eff + 2 


2y eff M eff " ( y eff " *) 



The solution 
p 3 /P0= 1.13 
heat ratio of 


of these equations for the nominal test conditions (e = 0.0826, 
x 10 -3 ) gives an effective Mach number of 5.78 and an effective specific- 
1.12. A comparison of calculated pressure distributions and shock shapes 


11 



based on ideal -gas theory at these effective conditions with experimental data is needed 
for verification of this method. 

Operational characteristics .- The nominal stagnation temperature and pressure 
(Tq = 736° K and Pq = 1650 N/cm2) resulted in approximately the largest value for the 
normal-shock density ratio that was obtained in the pilot facility. (See table n.) Since 
the specific-heat ratio of CF 4 decreases with an increase in temperature (fig. 11), a few 
attempts to obtain higher shock density ratios by testing at higher temperatures were 
made. However, the ratio of test-section area to throat area of the pilot CF 4 tunnel was 
fixed at a value of 3300 and the decrease in y associated with higher temperatures also 
caused a reduction in test-section Mach number. (See fig. 28.) The decrease in Mach 
number with a decrease in y offset the desired increase in shock density ratio so that 
a value of 12.1 was the highest realized in this pilot facility. The variation of area ratio 
with Mach number is shown in figure 28 for CF 4 at two different stagnation temperatures. 
This figure shows that if the Mach number is 6 at 589° K, an increase in stagnation tem- 
perature to 811° K while the area ratio is kept constant will result in a decrease in Mach 
number to 5.3. 

Some tests were made at stagnation temperatures as low as 275° K and stagnation 
pressures of about 1650 N/cm2. At these conditions the Mach number was about 8 and 
the normal-shock density ratio was about 10. Even at this low stagnation temperature, 
no condensation of CF 4 in the test-section flow was evident. Other tests were made to 
determine the vacuum required to operate the tunnel since CF 4 gives a very low test- 
section static pressure (fig. 13). At the nominal test conditions the pilot facility required 
a sphere pressure of 30 N/m2. 

Comparison of shock shapes in air and CF 4 .- A comparison of shock shapes for the 

120° blunt cone configuration obtained in the pilot CF 4 tunnel and the Langley Mach 8 
variable-density hypersonic tunnel is shown in figure 29. The CF 4 shadowgraph was 
taken at the nominal test conditions where the normal-shock density ratio was approxi- 
mately 12.1. The schlieren in the Mach 8 tunnel was taken in air at a Mach number of 8 , 
a specific-heat ratio of 1.4, and a normal-shock density ratio of 5.56. There are two 
notable differences in these shock shapes. At the higher shock density ratio in CF 4 , the 
standoff distance is considerably less. Also in CF 4 the shock over the conical portion of 
the body is almost linear and thus indicates a supersonic or near- supersonic flow in the 
shock layer downstream of the spherical nose as compared with the highly curved shock 
and corresponding subsonic flow in air. Although no force or pressure-distribution data 
have been obtained in this pilot facility (the small size of the model, 15-mm diameter, 
makes such measurements difficult), it is believed that the results discussed herein 
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demonstrate the feasibility and desirability of this type of facility for simulating real-gas 
effects on the aerodynamic characteristics of blunt high-drag entry vehicles. 

Design calculations for contoured nozzle .- As a result of these exploratory tests 
a Mach 6 nozzle with a 16.4-cm-diameter test section has been designed for CF 4 at stag- 
nation conditions of 811° K and 1724 N/cm2. a real-gas program was used to calculate 
the contoured axisymmetric nozzle, and a correction for boundary -layer displacement 
thickness was included. The details of this calculation procedure are discussed in 
appendix B. 


BALLISTIC -RANGE TESTS 

Ballistic ranges have the capability of obtaining high normal-shock density ratios 
in the actual gas of a planetary atmosphere, but to do this, they must operate with high 
free-stream pressures in order to achieve equilibrium composition of the gas in the shock 
layer because of the very small size of the models. These high pressures result in very 
rapid damping of the model motions and limit the stability data that can be obtained. If 
CF 4 is used as the test gas and the models are flown at much lower speeds so that no dis- 
sociation occurs, the desired high shock density ratio can be obtained at any test- section 
pressure and data can be obtained at conditions giving optimum motion of the model for 
data reduction. 

With these considerations in mind, an exploratory investigation in a ballistic range 
using CF 4 as the test gas was undertaken by the Naval Ordnance Laboratory with NASA 
support. This was done as part of an extensive investigation of the three shapes shown 
in figure 1. Most of the data were obtained in air; however, several launches were made 
in CF 4 to study the aerodynamics of these shapes at high shock density ratios. For these 
CF 4 experiments small amounts of air (5 to 10 percent) were present and the concentra- 
tion along the length of the range varied from almost pure CF 4 to mixtures of about 
90 percent CF 4 and 10 percent air. Since the properties of these mixtures may be some- 
what different from pure CF 4 , the present results should be used with some reservation. 

Launches were made at or near 460 and 1500 m/sec in CF 4 . Spheres as well as the 
three high-drag configurations were tested, and measurements of the shock standoff dis- 
tance for the spheres were used to verify the actual shock density ratio. The highest 
shock density ratios obtained in air at launch velocities of 5000 m/sec were about 9, 
whereas ratios as high as 17 were obtained in the CF 4 tests at launch velocities of about 
1500 m/sec. 

Drag is one of the most accurately measured quantities in a ballistic range test. 
When the drag force is nondimensionalized by the dynamic pressure i PjUj 2 to obtain 
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a drag coefficient, an error due to uncertainty in gas composition is introduced. This 
error is due to the uncertainty in the density of the gas in the test section. Since the 
pressure is kept constant by pumping and since air leaks into the test section between the 
time of measurement of the gas density and the time of launching, the density is actually 
somewhat less than the value used for data reduction. Therefore the drag-coefficient 
values computed from the measured data are somewhat less than the actual values at the 
time of launch. 

The measured drag coefficients of the three shapes in air and CF 4 are shown in 
figures 30 to 32. The drag coefficient is always higher for the high density ratios 
obtained in CF 4 . This increase is particularly large for the tension shell (fig. 32). The 
difference in shock shape for the tension shell in air and CF 4 is also more pronounced. 
Shadowgraphs of the flow field of the tension shell in air and in CF 4 at approximately the 
same Mach number are shown in figure 33. The flow over the tension shell is sensitive 
to boundary- layer separation which can alter the shock shape as well as the forces. A 
discussion of the Reynolds number and wall-temperature effect on separation and on the 
shock shape for this configuration is given in reference 8 . The results of reference 8 
and tests made in the ballistic range at different Reynolds numbers indicate that the 
boundary- layer flow is attached for both shadowgraphs shown in figure 33. Shadowgraphs 
of the flow over the 120° cone in both air and CF 4 are shown in figure 34. 

The increases in drag at the higher shock density ratios in CF 4 shown in fig- 
ures 30 to 32 are in agreement with the trends based upon the theoretical considerations 
discussed previously in the section entitled "Blunt-Body Flow-Field Simulation." It is 
apparent that useful studies of real-gas effects (high shock density ratios) on the aero- 
dynamics of blunt configurations in ballistic ranges can be made by using CF 4 as the test 
gas. 


CONCLUDING REMARKS 

An exploratory experimental and theoretical investigation of the use of tetrafluoro- 
methane (CF 4 ) as a test gas in wind-tunnel and ballistic- range facilities to simulate cer- 
tain flow-field characteristics of very blunt vehicles in high-speed flight where dissocia- 
tion takes place in the shock layer indicated the following: 

1. The effects of dissociation in the shock layer of very blunt vehicles on the aero- 
dynamic characteristics, such as shock shape, shock standoff distance, stagnation-point 
pressure level, surface pressure distribution, and forces and moments, can be simulated 
by matching the density ratio across the normal shock. 

2. Comparison of shock shapes of a 120° blunt cone configuration obtained in the 
pilot CF 4 tunnel at a shock density ratio of 12.1 with the shock shape obtained in a Mach 8 
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air tunnel at a shock density ratio of 5.6 and with shock shapes calculated by a one- strip 
integral method for ideal gases at a range of shock density ratios from 5 to 20 demon- 
strates the strong dependence of the flow field on shock density ratio and the validity of 
simulating high- shock-density- ratio flows by using CF 4 . 

3. The use of CF 4 as the test gas in both wind tunnels and ballistic ranges is a 
practical and relatively simple method of obtaining test conditions with large normal- 
shock density ratios. 

4. Tests made in a small pilot model CF 4 wind tunnel confirmed real-gas calcula- 
tions for test conditions in CF 4 . Measurements of pitot pressure and free- stream static 
pressure were in good agreement with calculations of isentropic expansion and flow 
across shock waves in CF 4 . Experimentally determined normal-shock density ratios 
obtained by measuring the shock standoff distance for a sphere were in agreement with 
calculated values. 

5. Values of normal-shock density ratio as high as 12.1 were obtained in this pilot 
CF 4 wind tunnel at a stagnation pressure of 1650 N/cm2, a stagnation temperature of 
736° K, and a free-stream Mach number of 6.1. 

6 . Exploratory tests using CF 4 in a ballistic range demonstrated that shock density 
ratios as high as 17 can be obtained at launch velocities of only 1500 m/sec, whereas the 
highest shock density ratios obtained on models of the same size in air was about 9 at 
launch velocities of 5000 m/sec. 

7. Measurements of drag for three different blunt configurations in the ballistic 
range indicated a definite increase at the high shock density ratios obtained in CF 4 . 
These results were in agreement with the trends of theoretical predictions. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., March 12, 1969, 

129-01-03-07-23. 
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APPENDIX A 


CALCULATION PROCEDURE FOR EVALUATING ISENTROPIC EXPANSIONS 
AND FLOW ACROSS NORMAL AND OBLIQUE SHOCKS 
IN TETRAFLUOROMETHANE 

By James L. Hunt, Kathryn A. Smith, 
and Robert B. Reynolds 

Tetrafluoromethane (CF 4 ) is a nonlinear polyatomic molecule (ref. 16). Its nine 
vibrational energy modes display anharmonic oscillator characteristics having four 
fundamental frequencies, two of which are threefold degenerate and one of which is two- 
fold degenerate. The acute temperature dependence of the vibrational energy of this 
molecule reveals itself in the specific heats at constant pressure and volume and their 
ratio y. Thus, an isentropic expansion in CF 4 is a variable y expansion, being neither 
calorically nor thermally perfect. The dependence of y on temperature along with its 
small (compared with temperature) dependence on density is vividly displayed in the dis- 
continuity in y that occurs across strong shocks. 

A calculation procedure for evaluating isentropic expansions and flow behind normal 
and oblique shocks in CF 4 is presented herein in the order in which it was programed for 
the digital computer. In this analysis the CF 4 was at all times assumed to be in thermo- 
dynamic equilibrium. At 298° K and 10.1325 N/cm2 the overall vibrational relaxation 
time of CF 4 is approximately 8.2 x 10"^ sec (ref. 17). The general theory for the relaxa- 
tion of a system of harmonic oscillators predicts that the relaxation time varies inversely 
with pressure to the first power and decreases monotonically with temperature (ref. 18). 
For polyatomic gases this temperature dependence is given, according to the present con- 
sensus, as e T ' 1//3 (refs. 19 and 20). For the temperature regime of interest in this 
report only a small fraction of the vibrational energy is stored in the anharmonic modes 
(0.025 percent at 100° K and 2.9 percent at 1000° K (ref. 16)). 

In each case for all the flow models studied, the use of the thermodynamic - 
equilibrium assumption was justified by a comparison of the vibrational relaxation dis- 
tances (velocity of flow multiplied by a local relaxation time calculated from the harmonic 
relation for the dominant vibrational mode) with flow lengths (radius of model forebody, 
nozzle Mach cone length, etc.). For example, 

(1) The relaxation times and distances for isentropic expansion in the Mach 6 , CF 4 
nozzle (design method, design conditions, and nozzle coordinates are presented in 
appendix B) at various axial stations are as follows: 
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Axial 

station, 

cm 

Mach 

number 

Relaxation time, 
sec 

Relaxation distance, 
cm 


0 

3.2 x 10' 9 


0 

1.0 

6.6 x lO- 9 

2.2 x 10" 4 

.04 

1.3 

8.3 x 10" 9 

3.1 x 10" 4 

2.9 

3.1 

2.5 x 10' 7 

1.9 x 10‘ 2 

7.6 

4.5 

7.0 x 10" 6 

6.3 x 10- 1 

65.9 

6.0 

2.2 x 10' 4 

2.2 x 10 1 


The fact that the axial length of the Mach cone in this nozzle (50.4 cm) is more than 
twice the relaxation distance at a Mach number of 6 (22.5 cm) indicates that equilibrium 
exists in the free stream at the test- section location. 

(2) Behind a normal shock at a Mach number of 6 (for the table HI nozzle) the relax- 
ation time is 3.3 x 10 - 6 sec and the relaxation distance is 0.026 cm. According to equa- 
tion (4), the shock standoff distance on a 3.2-cm-diameter sphere (average model base 
diameter tested) for the nominal test-section design conditions is 0.0958 cm, approxi- 
mately 3.7 times the relaxation distance. 

(3) Behind a normal shock in a ballistic range with an ambient pressure of 

2666 N/m 2 and temperature of 294° K for a projectile velocity of 1524 m/sec, the relax- 
ation time is 2.2 x 10' 7 sec and the relaxation distance is 1.7 X 10~ 3 cm. The smallest 
model fired at this velocity for the ballistic- range test discussed in this report had a 
base radius of 1.58 cm. A sphere with this radius would have a shock standoff distance 
(eq. (4)) of 0.0709 cm, more than 40 times the relaxation distance. 

SYMBOLS 

The following symbols are used in this appendix in addition to some of those defined 
in the body of the paper: 

A 2 ,A 3 ,A 4 ,A 5 ,A 6 constants in equation of state (eq. (Al)) 
a speed of sound 

a 4 constant in equation for specific heat at zero pressure (eq. (A4)) 

BgjBg.B^B^Bg constants in equation of state (eq. (Al)) 
b constant in equation of state (eq. (Al)) 
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b 4 constant in equation for specific heat at zero pressure (eq. (A4)) 

C 2 >C 3 ,C 5 constants in equation of state (eq. (Al)) 

c 4 constant in equation for specific heat at zero pressure (eq. (A4)) 

Cp specific heat at constant pressure 

c v specific heat at constant volume 

c v * specific heat at zero pressure 

d 4 constant in equation for specific heat at zero pressure 

g acceleration due to gravity 

J unit conversion factor 

j mechanical equivalent of heat 

K constant in equation of state (eq. (Al)) 

R gas constant 

5 entropy 

T c critical temperature 

V specific volume 

a constant in equation of state (eq. (Al)) 

6 flow deflection angle 

Ejj functional representation of conservation of mass (eq. (All)) 

functional representation of conservation of momentum (eq. (A12)) 
fijq- functional representation of conservation of energy (eq. (A13)) 
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APPENDIX A - Continued 


Subscripts: 


i point in isentropic expansion 

r reference 

N normal 

t tangential 


GENERAL EQUATIONS 


The equation of state of CF 4 and an expression for its viscosity were obtained from 
E. I. du Pont de Nemours & Company. The equation of state is 

RT A 2 + B 2 T + C 2 e" KT A3 + b 3 t + C 3 e KT + A 4 + B 4 T 
P = V - b + ( V - b)2 + (V - b) 3 (V - b) 4 


A 5 + B 5 T + C 5 e 


-KT 


(V - b) 5 


+ ( a 6 + B 6 T ) e 


aV 


where p is given in lb/in2, v in ft3/lb, and T in °R and where 

R = 1.219336 x 10- 1 
A 2 = -2.162959 
B2 = 2.135114 x 10 -3 
C 2 = -1.8941131 x 10 1 
A 3 = 4.404057 x IO- 3 
B 3 = 1.282818 x 10" 5 
C 3 = 5.39776 x 10 _1 
A 4 = 1.921072 x 10“ 4 
B 4 = -3.918263 x 10 -7 
A 5 = -4.481049 x 10 -6 
B 5 = 9.062318 x 10 -9 
C5 = -4.836678 x 10 -5 
Aq = 5.83882 x 10 7 
B 6 = -9.26392 x 10 4 
b = 1.5000 x IO- 3 
K = 9.76798 x 10" 3 
a = -6.61199 x 10 2 


(Al) 
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The following expression was used for the viscosity of CF 4 : 






1.22 
\ + 0.22 


(A2) 


where jj. is in micropoise and T c is the critical temperature (409.5° R). 
The specific heat at constant volume is 


c v ~ 



(A3) 


A third-degree expression in temperature was obtained for c v * from a curve fit of 
points calculated from two overlapping fourth-degree expressions for c v * presented 
in reference 16 for a temperature range of 540° to 9000° R and a third-degree expres- 
sion (obtained from E. I. du Pont de Nemours & Co.) for a temperature range of 180° to 
1170° R. Each of these three expressions for c v * were curve fits of points calculated 
from partition functions. The equation for specific heat at zero pressure c v * used 
herein was fitted over a temperature interval of 180° to 2500° R and is 

c y * = a 4 + b 4 T + c 4 T2 + d 4 T 3 (A4) 

where c v * is in Btu/lb-°R and 

a 4 = 1.90458084 x 10 -2 
b 4 = 3.00892783 x 10 " 4 
c 4 = -1.30237441 x 10 ‘ 7 
d 4 =1.96802894 x 10 " 11 


Evaluating the integral term of equation (A3) by the equation of state (Al) gives 


= c v * - JTK 2 e" KT + 


C 3 


v - b 2(V - b ) 2 4(V - b ) 4 

where J = 1.850505 x 10 _1 Btu-in 2 /(ft 3 -lb). 

The specific heat at constant pressure is 

[(3p/3T ) v ] 2 
(3p/3V) T 


(A5) 


c p = c v - TJ 


(A 6 ) 


ISENTROPIC EXPANSION 


Consider the isentropic expansion of CF 4 in a nozzle where the stagnation chamber 
pressure p 0 and temperature T Q are specified. (The equation of state specifies Vq.) 
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APPENDIX A - Continued 
The governing statement in the expansion is that the entropy remains constant: 

'V 0 

[£ E. 


Sq - s r - 


T ° Sl£I + j 


h~\ 

%] dV = Si - s r = 

^3T/ T y 


rT 


* c v dT 


+ J 


'Vt 


1 (!f) v dV (A7) 


where the subscript r represents reference conditions taken herein as 

T r = 820° R 

S r = 0.848 for <p r = 1 psia 

I V r = 100 ft 3 /lb 

This point was obtained from reference 21 where the datum enthalpy- entropy point was 
taken as zero at absolute zero for CF 4 in the solid state. 

The subscript i represents the tth point in the expansion. Here, the temperature 
to which the flow is to be expanded isentropically is designated Tj. 
modynamic variables T, and S 0 = S, completely specify a poin in 
which all other thermodynamic variables at this point, and hence flow quantities, may 

calculated. 

For simplicity of form, it is advantageous to express the thermodynamic variables 
at a point in the expansion as a function of Ti and Vi rather an 1 an 1 
value of Vi is calculated from the following integrated form of equation (A ). 

Si - S r . a 4 log e £ ♦ b 4 ( Ti - T r ) ♦ ^(T,* - T r 2) ♦ ^ - T r ’) 

+ JK (e- KT i-e- KT t) 


C 3 


C 5 


2(v r - b)2 4(V r - b) ' 


Vi - b> 


+ JR 10g 6 ( ^ 


, -KTA, , , \ ■>(% - KC 3 e ' KT ‘) 

r 1 

1 .1 

- J^B 2 - KC 2 e _ b v r - bj 


2 

j(Vi-b) 2 (V r - b)^J 

/ KT-\ 

j(b 5 - KC 5 e ‘Jr 1 l 1 

- JB, 

r 1 . 

1- 1 .1 


4 K V i' b ) 4 ( Vr ' b ) 4 


[3(Vi - b Y 3(V r - b y\ 



+ JB 6 ( 


«Vi «V r \ 


(A 8 ) 


a 


a 
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The enthalpy Hi expressed as a function of Ti and Vi is 


^ - H r = 


r Ti 


c v + V 


dT + 


VavL + 1 UtA 


dV 


and becomes in integrated form 


(A9) 


, T- 2 - T 2 t 3 f-n 4 ^4 

‘ H r ' « 4 (T, - T r ) + b 4 + c< Ii_L 


+ d 4 - A Tr + jfe KT ‘(KTj ♦ ll - e' KTr (KT r + l)l C ? + Ea__ + c 5 ~ 

4 L ' V F ^J[Vr-b 2(v r - b)* 4(v r - b) 4 


+ JV t 


vfrsFi-'r)* 


B 2 ( Tj - T r ) * C 2 (e’ KT i - + B 3 (T,-T r ),C 3 (e- KT i.e- KT r) + B 4 ( Tj - T r ) 

(V r - b) 3 


( V r - b ) 2 


(V r - b) 4 


B 5( T i - T r ) + C 5 (e- KT l - e - KT r) 

+ ^ , Ti" +B 6 e r (Ti^T r s 

(V r - b) 5 V 

+ 2(A 2 .B 2 T i + C 2 e- KT i)Lj_ s + 


- 




r 2(V r -b) 


2 V, - b 


+ 3(a 3 ♦ B 3 Tj + C 3 e' KT i)f_! + 


-1 1 6- 

A3 , 


L 3 ( V r ‘ b ) 4(v r -b) 4 3 (v*-b) 3 4(Vi-b) 4 

3 


2(Vi - b) : 

5 ( A 5 + B 5 T i + c 5 e KTl )| 


l )fl 1 + _ b . j b 1 

2(V r - b) 2 3(V r - b) 3 2(Vj - b) 2 2(v* - b)j 


fjy^ 

r-^T, 

1*1 

( v i - f 


+ A 6 ae > + BgTjae * + jvj 


BdTj aVi a \, 

3_L^ +Ae 1 + B R T i e £yV| 

(Vj-b) 4 6 6 1 


4(v r - b) 4 5(v r - b) 5 4(V, - b) 4 5(V, - b) 5 

a Q'V'r a 

Agae r + BgT|Cye 


- JV, 


B 4 T i . ffV r _ _ aV~ 

— + Age + BgTje 


L( V r - by 


JV r 2 

-4B 4 Ti 

2 

(V r - b) 5 ' 

' 1 

1 \ 

7i‘ b v r -b) 


KC 3 e' KT i 

r 1 i i 

. B r i ii 

V*V"* 

r i 

2 

(Vj-b) 3 (V r - b) 3 

]3(Vj-b) 3 3(V r -b) 3 

4 

\& 

i 


4(* ffVi -‘“ V 4 (A10) 


where H r = 200 Btu/lb at reference conditions (ref. 21). The total enthalpy Hq can be 
computed from equation (A10) by substituting the stagnation temperature and volume in 
place of the values at a given point in the expansion. 

For a given point in the expansion (assuming T t and Vi have been calculated 
for this point), the pressure can be obtained by using equation (Al) and the viscosity 
Mi, by using equation (A2). Other flow quantities of interest to this investigation at points 
in the expansion are given in the following equations: 
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Limiting velocity: 

= ^2gj(H 0 - H sv ) 

where H sv is the enthalpy of saturated vapor extrapolated (from points presented in 
ref. 21) to zero pressure and is approximately 2.27 x 10 5 joules/kg; g is the accelera- 
tion of gravity; and j is the mechanical equivalent of heat. 

Local velocity: 

Ui = \j2gj (H 0 - Hi) 

Local specific heat ratio: 

r i = ( C p)i/(° v )i 

where (c p ). and (c v ). are given by equations (A6) and (A5), respectively. 

Local speed of sound: 



Local Mach number: 

M i = u i/ a i 

The other flow parameters used in this report are easily obtained from the foregoing 
relations. 


FLOW BEHIND NORMAL AND OBLIQUE SHOCK WAVES IN CF 4 

In the previous section, CF 4 was expanded isentropically from specified stagnation 
conditions to a given Mach number. In this section, the flow behind two-dimensional, 
infinitesimally thin, normal and oblique shock waves when the process is always in a state 
of equilibrium will be considered. 

Static Conditions Behind Shock 

The following diagram illustrates the problem being considered and the notation 
that will be used throughout this analysis: 
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where 


(“l,i) sln 9 ' "w 
("l,i) cos e ' u l,t,i - "2,1,1 


The governing conservation equations are as follows: 
For mass, 

i.g. 


For momentum, 


For energy, 


v l,i v 2,l ' " N>1 


, 1 ("l.N.0 2 . - . 1 ( U 2,N,l) 2 . „ 

Pl,i + i V7“' P 2.‘ + B V 21 

9 

1,1 2gj 2. 1 2gj 


(All) 


(A12) 


(A13) 


Substituting u ? N . = E n ,V 2 4 from the conservation equation for mass into the conser 
vation equations for momentum and energy and noting that all thermodynamic variables 
have been previously expressed as functions of temperature and specific volume give 


P 






(A14) 


and 




(Tl,i,Vl,i,«i)] 2 (v 2;1 ) 2 

~~ 2gj 



(A15) 


where i» V 2,i) 

shock (eq. (A10)) and 


is the functional representation of the static enthalpy behind the 

d /t V n A is the functional representation of the static pressure 

*V2,i» 2,y 
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behind the shock (equation of state, eq. (Al)). Thus for given stream conditions desig- 
nated by expansion points in the previous section, equations (A14) and (A15) for a given 
shock-wave inclination angle (either normal or oblique shocks) are expressed in terms 

of two unknowns T„ . and V 0 .. These two equations were solved simultaneously for 

^,i <s,i 

the two unknowns with a digital computer program. 

The thermodynamic and flow quantities behind a particular shock wave are stipu- 
lated herein (once the independent variables T 2 i and V 2 i are known) by the following 
equations: 

Velocity: 

"2.1 ■ |“2,N,i ) 2 + ("2,t,i) 2 J 1/2 

where 

u 2,N,i = %( T l,i> V l,i> e i) v 2,i 
and 

2,t, i l,t,i 

Ratio of specific heats: 

c p( T 2,i> V 2,l) 
r2>i = Cv(r 2| i,V 2 .l) 

Speed of sound: 


a 2,i = 


'8P \ 


?P/T 


1/2 


J 2,i 


Mach number: 


M 


u 2,i 


2,i 


The functional relations expressed in these equations are given in the section of this 
appendix entitled "General Equations." 


Total Conditions Behind Shock 

With the static conditions behind a given shock wave determined as outlined in the 
previous section, the total conditions behind this shock can be obtained by bringing the 
flow to rest isentropically. The governing statement in this compression process is that 
the entropy (stipulated by the static condition behind the shock) remains constant; that is, 
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S 2 ,i = S(T 2 )i ,V 2ji ) = Constant = (A16) 

where T 3 j and V 3 j are the respective stagnation temperature and specific volume 
behind the ith shock. 

A second equation in terms of T 3 j and V 3 j which stipulates the final point in 
this compression process (the point at which the velocity is zero) is obtained from the 
fact that the total enthalpy is constant on either side of the shock. Thus, 

H 0 = H(T 0 ,V 0 ) = H(T 3 ,V 3 ). (A17) 

Equations (A16) and (A17) were solved simultaneously for the two unknowns with a digital 
computer program. Once T 3 i and V 3 i were determined, p 3 . was calculated from 
the equation of state: 

P 3 ,i = p(t 3 ,v 3 ). 


COMPUTER PROGRAMS FOR EVALUATING ISENTROPIC EXPANSIONS AND 
FLOWS ACROSS NORMAL AND OBLIQUE SHOCKS IN CF 4 


FORTRAN IV Control Data series 6000 computer programs were developed for 
evaluating isentropic expansions and flows across normal and oblique shocks in CF 4 in 
the format presented in the previous sections of this appendix. The abstracts and listings 
for these programs are given in appendix C. Inputs common (need to be input with first 
case only) to these programs are 


Symbol 

Machine n 

R 

R 

A 2 

A2 ^ 

a 3 

A3 

a 4 

A4 

a 5 

A5 

A 6 

A 6 

b 2 

B2 

b 3 

B3 

b 4 

BIGB4 

b 5 

B5 


Comments 

Gas constant, 1.219336 X 10" ^ (ft3)(psia)/lb-°R 


Constants in equation of state 
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Symbol 


Machine name 


Comments 


C 2 

C2 

'N 

C 3 

C3 


C 5 

C5 


b 

BLIT/ 

a 

AL 


-K 

XKK J 

a 4 

Z4 


b 4 

B4 


c 4 

C4 

/ 

d 4 

D4 


J 

XJ 


g 

GEE 

j 

AJA 

T c 

TZ 


T r 

TREF 

P r 

PREF 

S r 

SREF 

H r 

HREF 


Constants in equation of state 


Constants in curve fit of c v * 

1.850505 x 10" 1 Btu-in 2 /(ft 3 -lb) 

32.3 ft/sec 2 

778 ft-lb/Btu 

409.5° R 

820° R 

1.0 psia 

0.848 Btu/lb-°R 
200 Btu/lb 


Program A: Isentropic Expansion and Static Conditions Behind Normal and 

Oblique Shock Waves in Tetrafluoromethane 

The primary inputs to this program are the stagnation conditions for a particular 
isentropic expansion Pq and Tq and a tabulation of values of Ti (Ti < Tq) in a 
decreasing sequence which stipulates the stations to which the flow is to be expanded. 

For each of these Tj stations the flow quantities and thermodynamic properties are 
calculated and tabulated. 

For the shock-wave calculations, a minimum Mach number and a shock-wave angle 
6 are specified. At all Tj stations for which the local Mach number exceeds the spec- 
ified value, static properties behind a shock at the input 6 are calculated and tabulated. 
The program is set up so that the shock portion is activated when the stream Mach num- 
ber exceeds 3. 

Input 


Symbol 

Machine name 

Comments 

Po 

pee''! 

Total conditions from which the flow is expanded 

T 0 

TEE J 


Ti 

TEI 

Table listed in descending order 
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Input 


Symbol 

Machine name 

Comments 


IEND 

Number of entries in TEI table 


T 2,i 

TC"l 

"Estimates tables" correspond to values of Tj for Mi ? 3 

V 2,i 

VC J 




A 

Initial estimate ^ 



B 

Upper limit 



DELTX 

Iteration increment 



El 

Relative error criterion 

> Used to find Vq and Vj 


E2 

Absolute error criterion 



MAXI 

Maximum number of iterations 


6 

THETA 

Shock angle 


Output 





Symbol 

Machine name 

v 0 

VO ^ 


p 0 (Input) 

P0 


H 0 

HO 

> 

s 0 

SO 


^0 

GAMMA^ 


T r (Input) 

tref') 


p r (Input) 

PREF J 


Tj (Input table) 

TEMPERATURE 

Vi 

VOLUME 

Pi 

PRESSURE 

Hi 

ENTHALPY 

Mi 

MACH NO 

u i 

VELOCITY 

( c p)i 

CP 


( c v)j 

CV 


n 

GAMMA 


a i 

ACCELERATION 


[speed 

of sound] 

Pi 

VISCOSITY 

( R i/«)i 

REYNOLDS NO/F 

(Pl/ p o)i 

P1/P0 


(Tl/T 0 )i 

T1/T0 



Comments 


Total conditions 


Reference conditions 


~\ 


Stream static properties listed in corre- 
sponding sequence 


j 
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Output 

Symbol 

V 2 ,i 

T 2ji 

P 2 ,i 

u 2 ,i 

M 2 ,i 


H 


2 ,i 


Machine name 

VOLUME2 

TEMPERATURE2 

PRESSURE2 

VEL0CITY2 

MACH N02 

ENTHALPY2 


'\ 


p 2 ,i 

RH02 


( c v) 2 ,i 

CV2 


( c p) 2 ,i 

CP 2 


y 2 ,i 

GAMM2 

a 2 ,i 

SPEED 

( T 2/ T o)i 

T2/TCP 

( V l/ V 2 )i 

V1/V2 

( T 2/ T l)l 

T2/T1 

( p 2/Pl)i 

P2/P1^ 

u l,N,i 

uin' 


U M,i\ 

U1T 


u 2 ,t,ij 



u 2,N,i 

U 2 N^ 





Comments 


Static conditions behind shock at specified 8 


Static condition ratios across shock 


Velocity components parallel and normal 
to shock 


Program B: Ballistic -Range Normal Shock 

The normal-shock portion of the preceding program was used with input modifica- 
tions in evaluation of the normal-shock static -density, static -temperature, and static- 
pressure ratios given in figures 21 to 23 as functions of velocity in conjunction with 
ballistic -range tests conducted in CF 4 . The inputs to this program are range static tem- 
perature and pressure along with projectile velocity. 


Input 



Symbol 

Machine name 

Comments 

Ti 

TEfl 

Ballistic -range static conditions 

Pi 

PEI J 

> 

u i 

UI 

Projectile velocity J 


Corresponding sequence 


Number of sequential entries 


IEND 
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Input 

Symbol 

T 2,i 

V 2 ,i 


Machine name 



Comments 

"Estimates table" entries correspond to Ti, and ui 


Output 


Symbol 

Machine name 

Tj (Input) 

TEMPERATURE j 

Pi (Input) 

PRESSURE 

Vi 

VOLUME 

Uj (Input) 

VELOCITY 

Mi 

MACH NO 

«i 

ENTHALPY 

( c p)i 

CP 

( c v)i 

cv 

y i 

GAMMA 

a i 

ACCELERATION 


[speed of sound] 

Mi 

VISCOSITY 



T 2,i 

TEMPERATURE2 

P2,i 

PRESSURE2 

V 2 ,i 

VOLUME2 

u 2,i 

VELOCITY2 

M 2)i 

MACH N02 

=2,1 

ENTHALPY2 

^2,i 

RH02 

( c v) 2 ,i 

CV2 

( C P)2,i 

CP2 

y 2,i 

GAMMA2 

a 2,i 

SPEED OF SOUND 

p 2/ p l 

V2/V1~) 

T 2/ T l 

T2/T1 > 

P2/Pl 

P2/P1 J 


Initial estimate 
Upper limit 
Iteration increment 
Relative error criterion 
Absolute error criterion 
Maximum number of iterations 


Comments 


\ used to find values of Vi 


> Ballistic -range static conditions 


^Static conditions behind normal shock 


Ratio of static conditions across shock 
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APPENDIX A - Concluded 
Program C: Total Conditions Behind Shock 

For a given shock, the static temperature ^-2,i static specific volume V^i 
behind the shock determined from either of the preceding programs are used as inputs 
along with the stagnation temperature and pressure, in the case of the expansion problem, 
and the static temperature, static pressure, and projectile velocity, in the case of the 
ballistic -range problem, for a third computer program. This program was set up in the 
format presented in this appendix in the section entitled "Total Conditions Behind Shock 
to determine these conditions. 


Input 



Symbol 

Machine name 

Comments 

P 0 ’ p l,i 

PEE''| 

Total conditions for expansion case; static conditions for 

r *’ 0 » rp l,i 

TEE > 

ballistic -range case 

VoV M 

VEE J 



IEND 

Number of entries in following tables 


XMACH 



TTOT 

flnput zero for expansion case 

U 1 

Uzl 

\lnput velocity for ballistic- range case 

T 2,i 

TEI 


<3 

to 

VEI 


( p i/ p o)i 

PRAT 

Input 1 for ballistic -range case 

v 3 ,i 

VC 

Estimates 

T 3,i 

TC 

Estimates 

Output 



Symbol 

Machine name 

Comments 

T 3 

T2T 


V 3 

V2T 


p 3 

P2T 


S 2 

S2 

Entropy behind shock 

p 3 / p 0 

P2T/P01 

For ballistic- range case these ratios have the same value 

Of 

p 3/ p l 

P2T/P1J 

because of "ad hoc” substitution 

t 3 /t 0 

T2T/T0 

Becomes T 3 /T 1 for ballistic- range case 

v 3 /v 0 

V2T/V0 

Becomes V 3 /V 1 for ballistic -range case 

% 

HO 


s 0 

SO 
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HYPERSONIC NOZZLE CONTOURS FOR TET RAF LUOROME THANE 
By James L. Hunt and Lillian R. Boney 

Inviscid axisymmetric contours for a Mach 6 nozzle for tetrafluoromethane (CF4) 
at various stagnation conditions and maximum turning angles of the flow were calculated 
with the use of the computer program presented in reference 22. These calculations 
involved an application of the method of characteristics in which the real-gas variations 
of the thermodynamic properties of CF4 in an isentropic expansion were required inputs. 

As presented in reference 22, the viscous part of the program (Part in) that calcu- 
lates the boundary-layer displacement thickness along the nozzle which must be added to 
the inviscid contour to obtain the physical-wall contour is restricted to a gas composed 
of diatomic or linear polyatomic molecules. The vibrational energy modes of these mol- 
ecules display harmonic oscillator characteristics having only one fundamental frequency 

These "gas type" restrictions are imposed on the viscous part of the program 
through the expression for the local static enthalpy, which is written in reference 22 as 

h = ^ t + ^57t77 . (B1) 

where 0 is the characteristic molecular vibrational temperature and T is the static 
temperature. This equation includes real-gas effects due to vibration but no real-gas 
effects due to pressure; however, provisions are made in the program to account for 
real-gas effects due to pressure on enthalpy by multiplying equation (Bl) by a correc- 
tion factor expressed as a polynomial in pressure. 

Tetrafluoromethane is a nonlinear polyatomic molecule, the vibrational energy 
modes of which display anharmonic oscillator characteristics having four fundamental 
frequencies. Thus, in order to make the boundary -layer displacement-thickness part of 
the program applicable to tetrafluoromethane, the expression for the static enthalpy in 
the program must be made to comply with this molecular description. 

Since provisions are already made in the program to account for real-gas effects 
due to pressure through a correction factor rather than an equation of state, the enthalpy 
H must be expressed as a function of temperature only and may be stated as 

H = E + RT (B2) 

where R is the gas constant and the internal energy E is 

^ = ^tran + ^rot + ®Vib 
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From considerations of kinetic theory, the translation-energy contribution is 

Etran - f RT (B3) 

The rotational contribution for nonlinear polyatomic molecules is 


Erot = | RT 


(B4) 


In evaluating the vibrational contribution to the internal energy of CF 4 , anharmonics 
may be neglected safely for lower vibrational levels, that is, at lower temperatures for 
which the neglect of the interaction of vibration and rotation is also permissible. Calcu- 
lations based on expressions presented in reference 16 for both the harmonic and anhar- 
monic energy terms show that the harmonic assumption (neglecting anharmonics) induces 
an error in the total vibrational energy of only 2.56 x 10 " 2 percent at 100° K and 2.88 per- 
cent at 1000° K. Since the stagnation temperature in the nozzle being designed was only 
811° K, and in order to simplify the programing procedure, all vibrational contributions 
were assumed to be of a harmonic nature. The harmonic oscillator partition function 
for CF 4 is (ref. 16) 

Qvib = n(l-e- e ‘/ T )- S ‘ (B5, 

i=l 

where the degeneracies 


g for i = 

1, 2, 3, and 4 are 

g l = 1 

0 X = 2341° R 

(1300° K) 

g 2 = 2 

0 2 = 1127° R 

(626° K) 

g 3 = 3 

O 3 = 3322° R 

(1845° K) 

g 4 = 3 

0 4 = 1653° R 

(918° K) 


Now 


^vib 


= RT' 


d ( Iog e Qyib) 
dT 


(B 6 ) 


The vibrational contribution to the internal energy is 

4 
E 


r, V g i 0i 

Svib = R A, e Q i/ T - 1 

i=l 


(B7) 


Thus, the expression for the enthalpy of tetrafluoromethane, neglecting anharmonics and 
Van der Waal effects, becomes 

4 

(B 8 ) 


S= 4T + 
R 


V g i 0i 

- 1 


33 



APPENDIX B - Continued 


All the changes (approximately 21 equations and "if" statements) made to the vis- 
cous part (Part HI) of the nozzle design program of reference 22 to adapt it for use herein 
with CF 4 stem from the enthalpy expression in equation (B 8 ). The form of the turbulent 
skin-friction law as used in reference 22 was left unchanged. 

Part m of the FORTRAN IV IBM 7090 computer program presented in reference 22 
which computes the displacement thickness along a nozzle and applies it to the inviscid 
results to yield a physical contour, as modified to comply with the molecular description 
of CF 4 , is given in program D of appendix C. The flow properties at the edge of the 
inviscid region are supplied by Parts I and II of the parent program, the inputs and out- 
puts having the same form as Part III. (For additional information, see ref. 22.) 

Table in gives the coordinates of a Mach 6 nozzle with a test- section diameter of 
16.4 cm designed with this program for CF 4 at a stagnation temperature of 811° K, a 
stagnation pressure of 1724 N/cm2, and a turning angle of 18°. A sketch of this nozzle 
is shown in figure 35. 
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TABLE HI.- CF4 NOZZLE COORDINATES 

[Mach number, 6; stagnation pressure, 1724 N/cm2; stagnation 
temperature, 811° K; turning angle, 18^ 


in. 


0.0000 

.1000 

.2000 

.3000 

.4000 

.5000 

.6000 

.7000 

.8000 

.9000 

1.0000 

1.1000 

1.2000 
1.3006 

1.4000 

1.5000 
1.6000 

1.7000 
1.8000 

1.9000 

2.0000 

2.1000 

2.2000 

2.3000 

2.4000 

2.5000 
2.6000 

2.7000 
2.8000 

2.9000 

3.0000 


cm 


0.0000 

.2540 

.5080 

.7620 

1.0160 

1.2700 

1.5240 

1.7780 

2.0320 

2.2860 

2.5400 

2.7940 

3.0480 

3.3020 

3.5560 

3.8100 

4.0640 

4.3180 

4.5720 

4.8260 

5.0800 

5.3340 

5.5880 

5.8420 

6.0960 

6.3500 

6.6040 

6.8580 

7.1120 

7.3660 

7.6200 


in. 


0.0637 

.0656 

.0680 

.0718 

.0771 

.0840 

.0930 

.1045 

.1196 

.1397 

.1640 

.1913 

.2205 

.2511 

.2827 

.3148 

.3475 

.3805 

.4137 

.4471 

.4806 

.5142 

.5480 

.5817 

.6154 

.6491 

.6822 

.7148 

.7470 

.7788 

.8101 


cm 


0.1618 

.1666 

.1727 

.1824 

.1958 

.2134 

.2362 

.2654 

.3038 

.3548 

.4166 

.4859 

.5601 

.6378 

.7181 

.7998 

.8827 

.9665 

1.0508 

1.1356 

1.2207 

1.3061 

1.3919 

1.4775 

1.5631 

1.6487 

1.7328 


y (inviscid) 

Mach 

number 

in. 

cm 

0.0625 

0.1588 

1.00 
















.1027 

.2609 

2.295 




.1616 

.4105 

2.920 

.2458 

.6243 

3.361 










.4046 

1.0277 

3.859 







.566 

1.4376 

4.202 













.724 

3.2766 

4.430 

1 .786 

1.9964 

4.500 
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TABLE m.- CF 4 NOZZLE COORDINATES - Continued 


y (inviscid) 


Mach 

number 


3.1000 

3.2000 

3.3000 

3.4000 

3.5000 

3.6000 

3.7000 


4.2000 

4.4500 

4.7000 

4.9500 

5.2000 

5.4500 

5.7000 

5.9500 

6.2000 

6.4500 

6.7000 

6.9500 

7.2000 

7.4500 

7.7000 

7.9500 

8.2000 

8.4500 

8.7000 

8.9500 
9.2000 

9.4500 

9.7000 

9.9500 
10.2000 
10.4500 


7.8740 0.8410 


8.1280 

8.3820 

8.6360 

8.8900 

9.1440 

9.3980 


3.9500 10.0330 


10.6680 

11.3030 

11.9380 

12.5730 

13.2080 

13.8430 

14.4780 

15.1130 

15.7480 

16.3830 

17.0180 

17.6530 

18.2880 

18.9230 

19.5580 

20.1930 

20.8280 

21.4630 

22.0980 

22.7330 

23.3680 

24.0030 

24.6380 

25.2730 

25.9080 

26.5430 


.8715 

.9017 

.9314 

.9606 

.9895 

1.0180 

1.0876 

1.1548 

1.2198 

1.2829 

1.3440 

1.4032 

1.4605 

1.5160 

1.5699 

1.6220 

1.6728 

1.7223 

1.7703 

1.8170 

1.8623 

1.9063 

1.9491 

1.9906 

2.0310 

2.0703 

2.1085 

2.1458 

2.1822 

2.2176 

2.2520 

2.2855 

2.3181 


4.4! 
4.6: 
4.7: 
4.8' 
4.9507 


5.2586 


2.2936 


5,271 






APPENDIX B - Continued 

TABLE HI.- CF 4 NOZZLE COORDINATES - Continued 
~ ! I y (inviscid) 


Mach 

number 


in. 

cm 

in. 

cm 

10.7000 

27.1780 

2.3498 

5.9685 

10.9500 

27.8130 

2.3806 

6.0467 

11.2000 

28.4480 

2.4106 

6.1229 

11.4500 

29.0830 

2.4397 

6.1968 

11.7000 

29.7180 

2.4680 

6.2687 

11.9500 

30.3530 

2.4956 

6.3388 

12.2000 

30.9880 

2.5223 

6.4066 

12.4500 

31.6230 

2.5484 

6.4729 

12.7000 

32.2580 

2.5737 

6.5372 

12.9500 

32.8930 

2.5982 

6.5994 

13.2000 

33.5280 

2.6222 

6.6604 

13.4500 

34.1630 

2.6456 

6.7198 

13.7000 

34.7980 

2.6684 

6.7777 

13.9500 

35.4330 

2.6907 

6.8344 

14.2000 

36.0680 

2.7123 

6.8892 

14.4500 

36.7030 

2.7334 

6.9428 

14.7000 

37.3380 

2.7539 

6.9949 

14.9500 

37.9730 

2.7738 

7.0455 

15.2000 

38.6080 

2.7932 

7.0947 

15.4500 

39.2430 

2.8120 

7.1425 

15.7000 

39.8780 

2.8303 

7.1890 

15.9500 

40.5130 

2.8481 

7.2342 

16.2000 

41.1480 

2.8653 

7.2779 

16.4500 

41.7830 

2.8820 

7.3203 


42.4180 

43.0530 

43.6880 

44.3230 

44.9580 

45.5930 

46.2280 

46.8630 

47.4980 


2.8982 

2.9139 

2.9292 

2.9440 

2.9583 

2.9721 

2.9855 

2.9985 

3.0111 


2.330 5.9182 5.605 


7.3614 

7.4013 

7.4402 

7.4778 

7.5141 

7.5491 

7.5832 

7.6162 

7.6482 


2.679 6.805 
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TABLE III. - CF4 NOZZLE COORDINATES — Concluded 


X 

y 

y (inviscid) 

Mach 

in. 

cm 

in. 

cm 

in. 

cm 

number 

18.9500 

48.1330 

3.0232 

7.6789 



_ _ _ _ 



19.2000 

48.7680 

3.0350 

7.7089 

— 

_ _ 


19.4500 

49.4030 

3.0463 

7.7376 



_ _ 


19.7000 

50.0380 

3.0571 

7.7650 

— _ _ 

_ _ _ _ 



19.9500 

50.6730 

3.0677 

7.7920 

— 

_ _ _ _ 



20.2000 

51.3080 

3.0779 

7.8179 

— 

_ _ _ _ 

— — ^ — 

20.4500 

51.9430 

3.0878 

7.8430 

_ 

_ _ _ _ 


20.7000 

52.5780 

3.0973 

7.8671 

— . 

_ _ _ _ 


20.9500 

53.2130 

3.1064 

7.8903 


_____ 


21.2000 

53.8480 

3.1152 

7.9126 



_ _ _ _ 


21.4500 

54.4830 

3.1237 

7.9342 


_ _ _ _ 


21.7000 

55.1180 

3.1319 

7.9550 

2.777 

7.054 

5.944 

21.9500 

55.7530 

3.1397 

7.9748 

— 

_ 

w _ 

22.2000 

56.3880 

3.1473 

7.9941 

— 

_ _ _ _ 


22.4500 

57.0230 

3.1545 

8.0124 

— 

_ 



22.7000 

57.6580 

3.1615 

8.0302 


_ _ _ _ 


22.9500 

58.2930 

3.1682 

8.0472 



_ _ _ _ 


23.2000 

58.9280 

3.1746 

8.0635 

— - - 

_ 



23.4500 

59.5630 

3.1807 

8.0790 

— 

_ _ _ _ 

««i « m » 

23.7000 

60.1980 

3.1868 

8.0945 

_ . _ . 

, 



23.9500 

60.8330 

3.1926 

8.1092 

— 

_ _ _ _ 



24.2000 

61.4680 

3.1956 

8.1168 

— « _ 

_ _ _ _ 

v w M 

24.4500 

62.1030 

3.2010 

8.1305 

— 

_ _ _ _ 


24.7000 

62.7380 

3.2062 

8.1437 

_ . 

_____ 



24.9500 

63.3730 

3.2112 

8.1564 


_ _ _ _ 

— m 

25.2000 

64.0080 

3.2160 

8.1686 

— 



^ ^ 

25.4500 

64.6430 

3.2206 

8.1803 

— - - 

_ _ _ _ 

^ — m # — 

25.7000 

65.2780 

3.2250 

8.1915 

— 

_ 

_ w ^ ^ 

25.9500 

65.9130 

3.2291 

8.2019 

2.801 

7.115 

5.999 
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DIGITAL COMPUTER PROGRAM LISTINGS 

Page 

Program A - Isentropic Expansion and Static Conditions Behind Normal and 

Oblique Shock Waves in Tetrafluoromethane 40 

Program B - Ballistic-Range Normal Shock Program 56 

Program C - Total Conditions Behind Shock 68 

Program D - Hypersonic Nozzle Contours for Tetrafluoromethane 79 
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Program A 

NASA-LANGLEY RESEARCH CENTER 


EAR E 1263 


COMPUTER PROGRAM ABSTRACT 


01 TITLE 0T PROGRAM CHARACTERS MAXIMUM' 

Isentropic Expansion and Static Conditions Behind normal 
and Oblique Shock Waves in Tet rafluoromethane 


10/68 


PARENT PROGRAM 

I 0? I3SITE jo? 16 PROGRAM NO 


FtfR6M 


WHO TO CONTACT ABOUT THE PROGRAM 

1 O', ?B _ SIT L I O', 31 ORCN CODE IcT 


K. A. Smith 


DATES 

[OS 34 COMPLETED 


10/66 


10/68 


LR 11.170 


iK Sfc REVISION CODE 


I — I a revision 

I \ B CANCELLATION 


[cs 39 PROJECT NO 


RGL-168 


Q A THIS PROGRAM 

l CD A UN oER DEVELOPMENT ISNOTfOR 

tcr SHARING 

1 I B OPERATIONAL 

OD C COMPLETEO I 


TIME AND COST FOR DEVEl .QPWENT 

OS hi MACHINE IDS 69 COMPUTER *TVPE I OS 74 TOTAL COST 


[5? 6061 62 63 6S 66 6? 68 1 


g ABSTRACT 

The re 1st i on s for isentro p ic ex pa nsions and sta tic condi* — i — 

tions beh ind norma l and ob lique Bhoc^VBves_in. equilibrium, — i 
tet rafluoromethane. axe programed^to yield_. the static thermo- 4- 
.dynamic properties of CF 4 in high-speed, one-di m ens ion al flow. 4 - 
The pr imar y inputs to this program are the stagnation conditions. — 
a. given expansion and a tabulation of- t e mpe n a t ures t o. vhi eh — | — 
. the_f!oy._is to he expanded .and processed_thxough_a specified— 
.shock. 
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Program A - Continued 


PROGRAM E 1 263 ( I NPUT , OUTPUT , TAPES = I NPUT # TAPE6 = OUTPUT t TAPEd ) 

SIBFTC ISEERN DECK, LIST 

COMMON A 2 « A3»A4 ,A5#62*B3tt34 * 35 ♦ C2 ♦ C3 ♦ C4 # C5 # XKK «R , TREF ♦ VREF • 0L I T # X J 
COMMON D4 , 24 , ESS * TUrF ♦ CLOD • GEE t A J A , KX « P2 I < 30 ) i U2 I ( 30 > * H2 I ( 30 ) 
COMMON PT2 ( 30 ) ,RH0 { 30 ) * TE2 1(30) * VE2 1(30)# CV2 < 30 > «CP2 (30 ) t G AM2 ( 30 ) 
COMMON ACC2 ( 30 ) « XMACH2 ( 30 ) « RORAT ( 30 ) * TR A T ( 30 ) ♦ PRAT (30) * P TR AT ( 30 ) 
COMMON TH 

COMMON VI N ( 30 ) , V2N ( 30 ) i VI T ( 30 ) 

COMMON B I Gb4 , B6 , A6 , Al i TZ 
COMMON/ CL UX/TC t VC 
COMMON/DESP/VOODOO, ZULU 
EXTERNAL FOP 
DIMENSION T C ( 30 ) < VC ( 30 ) 

D I MENS ION PI OR AT ( 3C ) , T 1 ORAT ( 30 > ♦ REN ( 30 ) 

COMMON /RATS /PEE « TEE, V I S2 C 3 0 ) , PT20RT ( 30 ) * P 1 T 2R T ( 30 ) • T20RT ( 30 ) 
EQUIVALENCE ( C3 ♦ G3 ) ♦ <C5#XJ3># (C2*Q3) 

PO (XL# AT ) =1 ./EXP (XL* AT ) 

CIMEN5I0N TE I ( 30 ) # VE I ( 30 ) # PE I ( 30 ) # V I ( 30 ) #GAM ( 30 ) ♦ SPEED ( 30 ) 

C I MENS I ONXMACH ( 30 ) # CP ( 30 ) * C V ( 30 > « H < 30 ) t V I S ( 3 0 ) 

NAMELIST /DAT A/ A2 ,A3,A4,A5#A6#B2*83#BIG64«d5#B6#B4,AL# 
1C2#C3#C4,C5# XKK , XJ *R « BL I T # D4 , PEE * TEE • GEE # A JA # NUMB # I NT *Z4 # A *3 * 
2DELTX » I END * TE I #E1 , E2 , M A X I # TREF # PREF * TZ ♦ T C # VC ♦ THETA 
99 READ (5# DATA ) 

WR I TE ( 6 # DAT A ) 

WR I TE ( 6 # 88 ) A 2 , A3 , A4 , A 5 , B2# S3 # B4 .85 
QQ FORMAT ( 10H0A2 5 62 5//SE13.5) 

WRITE (6 #89 )C2,C3#C4iC5«XKK#XJ#R#dL IT 

89 FORMAT ( 20H0C2 5 K JR B/8E13.5) 

WR I TE ( 6 # 90 )D4 , PEE # TEE #GEE# A JA , NUMB « I NT « Z4 

90 FORMAT ( 23H0D4 PC TO G J N I /5E 1 3 • 5 # 2 I 1 3 * E 1 4 • 5 ) 

KX = 1 

TH = THET A*Q • C 1 74 532925 

V00D00 = TEE 

ZULU=PEE 

CALL ITR2(SP,A,B,0ELTX,F0P#E1 f E2#MAxI , I CODE ) 

1 VEE = SP 
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Program A - Continued 


V OODOO= TREF 
ZULU=PREF 

CALL ITR2 (SP , A , B « DELTX , FOP. El itZ.MAXl . ICOOt) 

2 VREF = 5P 

WR I TE (6.91 >VREF 
91 FORMAT (5H0VRFF/E1 A .5 ) 

CVO=CEV (VEE.TEE ) 

CP0=CEP (VEE. TEE J+CVO 

3 GAMO=CPO/CVO 

4 ESS=ENT (VEE.TEE ) 

CALL EQU3tTEE.VEE.ANS3) 

5 HO = ANS3+200 . 

UZ=SQRT (2.0*GEE*A JA*H0 ) 

VIS0=149.5* (TEE/TZ )*SQRT (TEE/TZ )*1 . 22/ < TEE/TZ+ .22 ) 
VISO=0.672E-7*V ISO 
VE=SQRT <2.*GEE*A JA*ho ) 

WRI TE (6.55 )VE 
55 FORMAT (3H0VL///E1 5.5) 

DO 111 N= 1 « I END 
TUFF=TE I (N) 

CALL EQUA2(A,B. DELTX, El .E2.MAXI ,ANS) 

6 VEI (N)=ANS 

PEI <N)=PRES<VEI (N ) * TE I (N) ) 

CALL EQU3 ( TE I ( N ) . VE 1 ( N ) , ANS3 ) 

7 H (N ) =ANS3+200. 

VI (N)=SQRT <2.*GEE*AJA*A6S <H0-H(N) ) ) 

CV(N)=CEV(VEI (N).TEl <N) ) 

CP (N ) = CEP(VEI ( N ) » TE I (N) )+CV(N) 

8 GAM (N ) =CP ( N ) /CV (N ) 

( N ) =SQRT (GAM(N)-*PART (VEI (N).TEI (N))*4533» C56 ) 
XMACH ( N ) =V I ( N ) /SPEED ( N ) 

VIS (N) = 149.5* (TEI ( N ) /TZ > *SQRT ( TE 1 (N)/TZ>* 1 .22/(TEl 

9 VIS(N)=0.672E-7*V IS(N ) 

PIC RAT ( N ) =PE 1 ( N ) / PEE 

T 1 ORAT ( N ) =TE I (NJ/TEE 
10 REN ( N ) = V 1 ( N ) / ( VE 1 ( N ) *V 1 S ( N ) ) 

IF(XMACH(N) • GT .3. .AND. IE ND .LT.31 ) CALL SHOCK ( H ( N ) , V 

2 ( N ) . TE I (N) ) 

111 CONTINUE 
KQ=KX-1 

WRITE(6.52)VEE.PEE,H0.ESS»GAM0,TREF,PRtF, VlSO 

52 FORMAT ( 3H 1 VO . 1 5X , 2HP0 « 15X.2HH0. 15X.2HS0, 15X, 5HGAMMA 

2 « 4HPREF , 15X.4HMU-0//8E16.6) 

]jg(5,5g,)(TEI ( N ) « VE 1 ( N ) , PE I (N),H(N),XMACH(N)»N=1, 
56 FORMAT ( 1 2H0 TEMPER AT URE , 1 5X , 6HV0LUME , 1 5X , 8HPRESSURE , 


(N)/TZ+.22) 

I (N) » PE I (N) , VE I 

, 15X.4HTREF, 15X 
I END ) 

1 5X .ShENTHALP Y . 
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215X.7HMACH NO/ // ( 5E20 • 6 > ) 

WR I T E ( 6 . 54 ) (VI ( J ) . CP ( J ) * CV ( J } . GAM { J } , SPEED (J ) . V I S ( J ) . J = i * I END ) 

54 FORMAT ( 9 H 1 VELOC I T Y . 20X « 2HCP . 2 OX .2HC V . 20X . 5HGAMMA . 1 5X . 1 2HACCELERAT I 
20N . 1 5X . 9HV I SCOS I T Y/// (6E2 1.6)) 

WRITE (6.60 ) (PI OR AT (N) . TIORAT(N) I END) 

60 FORMAT( 1HI .7X.5HP1/P0. 1 5X.5HT1/T0. 10X. 14HQ5YN0LDS NO/FT/ 

1 ( 3E20.6 ) ) 

I F ( I END.LT .31 ) WR I TE < 6 * 57 ) ( VE2 1 (K ) . TE2 I (< ) . P2 I ( K ) . U2 I ( K ) . XMACH2 ( < > ♦ 

2H2 I (K ) t K = 1 .<0 ) 

I F ( I END.LT .31 ) WH I TE < 6 . 58 > ( RHO (K ) . CV2 (»< ) . CP2 ( K ) «GAM2 (M • ACC2 (K ) • 
IKsliKQ) 

I F ( I END.LT. 31 ) WRI TE (6 . 59 ) ( T20RT ( X ) ♦ RORA T ( K ) . THAT <K ) . PkAT (K) iKsi, 

1KQ ) 

57 FORMAT ( 8H 1 V0LUME2 «20X. 1 2HTEMPERATURE2 . 1 2X . 9HPRES5URE2 t 1 2X. 9HVEL0C I 
2 TY2. 12X.8HMACH N02. 1 2X . 9HENTHALPY2/// ( 6E2 l .5 ) ) 

58 FORMAT ( 1H1 . 8X . 4 HRH02 . 2 OX . 3HC V2 . 16X.3HCP2. 1 6X * 6HGAMM2 . 16X. 

1 1 4HSPEED OF SOUND// ( 5E 2 1 . 5 ) ) 

59 FORMAT ( i H 1 . 7X . 5HT2/T0 . 2 1 X . 5HV 1 / V2 . 2 OX , 5HT2/T 1 * 1 6X. 5HP2/T 1 /// 

1 (4E21 .5 ) ) 

I F ( I END.LT .-3 1 ) WR I TE ( 6 . 62 ) ( V 1 N < K ) . V 1 T ( < ) . V2N ( K ) . VIS2(K)«K = 1«KG) 

62 FORMAT < 1 H I . 9X . 3 HU 1 N . l 8X « 3 H’J 1 T . 1 8X . 3HU2N . 1 8X . 3HMU2/// ( 4E2 1 .5 ) ) 

WR 1 TE (8 ) HO . VEE .GEE . KQ . I END . UZ . V I SO 

WR I TE (8 ) ( (GAM ( I ) ,P 10RAT ( I ) . V I ( I ) . Tl ORAT ( F > . REN ( I ) . VE 1(1). 

1 XM ACH ( I ) ) * I = 1 ♦ I END ) 

WR I TE (8 ) ( (GAM2 ( I ) .RORAT ( I ) . VE2 I ( I ) . U2 I ( I ) ♦ V I S2 ( I ) ) . I = 1 . KQ ) 

ENDFILE 8 
GO TO 99 
END 

FUNCTION FOP(V) 

C EQN OF STATE-USED TO F I ND V AS A FUNCTION OF P+T FOP 

COMMON A2 . A3 , A4 . A5.B2 . 83.tJ4 .05 .C2 . C3 . C4 . C5.XKK .R . TREE . VREF .8LI T . XJ 
COMMON D4 . Z4 . ESS . TUFF .CLOD. GEE . A J A . KX . P2 I <30)*U2I ( 30 ) . H2 I (30) 

COMMON PT2 (30 ) .WHO ( 30 ) . TE2 1(30) . VE2I (30 ) .CV2 ( 30 ) .CP2 ( 30 ) ♦ G A M 2 ( 30 ) 
COMMON ACC2 ( 30 ) . XMACH2 ( 30 ) . RORAT ( 30 ) . THAT (30) .PRAT(30) .PTH ATI 30) 
COMMON TH 

COMMON VI N ( 30 ) . V2N ( 30 ) . V 1 T ( 30 ) 

COMMON 8 I GB4 .86.A6.Al.TZ 
PO (XL. AT ) = 1 ./EXP (XL* AT ) 

C OMMON/DE SP/ VOODOO. ZULU 
T = VOODOO 
P = ZULU 


POT = P0 ( XK< . T ) F0 P 

1 FOP = P-(R#T)/( V-8L I T ) - ( A 2 +8 2* T+C2*P0T )/(V-BLlT)+*2-(A3+63*T+ FOP 

1 C3*P0T ) / ( V-BL IT ) * * 3- ( A4 f-B I GB4*T ) / (V-BL I T ) **4- ( A5+B5* T + C5 *POT )/ FOP 

2 ( V-BL I T >**5- ( A6-B6*T ) *EXP ( Al_*V ) FOP 
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RETURN 

END 

SUBROUTINE EQU3(TE.VE. A NS 3 ) 

ENTHALPY 

COMMON A2.A3.A4 . A5 . B2 .B3.B4.B5. C2 i C3 » 
COMMON D4 . Z4 »ESS» TUFF .CLOD. GEE iAJAiKX 
COMMON PT2 (30)i RHO ( 30 ) * TE2 I ( 30 ) . VE2 I ( 
COMMON ACC2 ( 30 ) . XMACH2 ( 30 ) i RORAT ( 30 ) « 
COMMON TH 

COMMON V 1 N ( 30 ) f V2N ( 30 ) * VI T ( 30 ) 

COMMON BIGb4tB6tA6 « Al*TZ 
EQUIVALENCE <C3 . G3 ) • (C5.XJ3 ) . (C2«D3 ) 


EQU3 

C4 . C5 . XKK . R . TREF ♦ VREF <dLIT«XJ 
. P2I (30) * U2 I (30 ) ♦h2I (30) 

30 > * CV2 ( 30 > .CP2 ( 30 > .GAM2 ( 30 ) 

THAT (30 ) « PRAT ( 30 > « PTW AT (30 ) 


PO* ( XL . A T ) = 1 ./EXP (XL* AT ) 

3 A NS 3 = 24* ( TE-TREF ) +84* (TE**2-TREF**2)/2.+C4*< TE**3-T*EF**3 ) /3. +D4 * 
2 ( TE**4-TREF**4 ) /4 •+ (PO (XKK , TE ) * ( XK<*TE * 1 • ) "PO ( XKK . TREF ) * (XKK* TREF 
3+1 • ) )*(XJ*D3/(VE-BLIT)+XJ*G3/ (2** (VE-BL 1 T)**2 ) +X J*X J3/ ( 4 . * (VE- 
4BL I T ) **4 ) ) +X J* (VE*R* (TE-TREF )/ ( VE-BL I T ) +VE* ( B2* ( TE-TREF ) +C2* ( PO 
5 (XKK * TE )-PO ( XKK * TREF ) ) > / ( VE-BL I T ) ** 2 ) +X J* ( VE* ( B3* ( TE-TREF ) +C3* £ PO 
6(XKK,TE)-P0(XKK,TREF) ) )/< VE-BL I T ) **3 ) +XJ* (VE* < 65* ( TE-TREF ) +CE* (PO 
7 ( XKK * TE ) -PO ( XKK .TREF ) ) ) / (VE-BL I T ) **5 ) - ( R* TREF * (ALOG < ( VE-BL IT )/( VR 
8EF-BLI T.) )+BL IT/ ( VREF-BL I T J-BLI T/ (VE-BL IT ) ) )*X J- (2.* ( A2 + 32*TREF + C2 
9*P0 (XKK .TREF ) )* ( 1 ./(VREF-dL I T ) +BL I T/ ( 2 • * ( VRLF-BL I T ) **L )-!•/( VE— 


01 

1 


ABL I T ) -BL IT/ ( 2 . * ( VE-BL IT)** 2) ) )*XJ-(3.*( A3+o3* TREF + C3*P0 ( XKK . TREF 
B))*<l«/(2**< VREF-BL I T ) **2 ) +BL I T/ ( 3 • * ( V«EF-BL I T ) +*3 >- 1 . / (2. * ( VE- 
CBL I T ) **2 ) — dL I T/ ( 3# * ( VE-BLl T )**3> ) ) *X J- ( 4 • *A4* ( 1 • / ( 3 • * ( VREF-8L I T ) 
0**3 >+BL I T/ ( 4 •* ( VREF-BL IT)**4)-1 ./(3.*( VE-BL I T ) **3 ) -BL I T/ ( 4 •* ( VE- 
EBL IT ) ** 4 ) ) )*XJ 

ANS3= ANS3- < (5.* (A5+B5*TREF+C5*P0 (XKK. TREF )))*(l./£4.*( VREF-BL I T ) 
2**4 ) +BL I T / (5 •* ( VREF-BL IT>**5)-1 ./(4.*(VE-BLIT)**4) -BL I T/ (5 •* < VE- 
3BL I T ) **5 ) ) ) *XJ+TREF*R*ALOG ( ( VE-BL 1 T ) / ( VREF-BL IT) )*XJ-TREF* (B2-XKK 
4*C2*P0 ( XKK . TREF ) ) * ( 1 * / ( VE-BL I T ) -1 • / ( VREF-BL IT)) *XJ-TREF* (B3-XKK* 
5C3*P0 ( XKK , TREF ) )/2. * ( 1 • / ( VE-BL I T >**2-1 •/ ( VREF-BL I T ) **2 ) *XJ-TREF* 

6 (B5-XKK*C5*P0 ( XKK . TREF ) )/4 . * ( 1 • / < VE-BL I T > **4- 1 • / ( VREF-oL I T ) **4 ) 

7* XU 

EXAL = EXP ( AL* VE ) 

alvr=al*vref 

EXAR=EXP ( ALVR ) 

CONTINUE 

TEMP = XJ*VREF* ( ( B I GB4* ( TL-TRLF )/( VRE»- -BL 1 T ) **4 ) +B6*EXAK* (TE-TREF ) 

DUM = B I GB4*TE/ ( VE-BL IT > **4 

DUM1 =A6*ExAl+B6*TE*LXAL 

TEMP=TEMP+XU*VE* ( 0 • 5*VE* ( "4 • *DUM/ ( VE-BL I T ) + AL*DUM1 )+ DUM+DUM1 ) 

DUM = B I GB4*TE/ (VREF-BL I T ) **4 

DUM1 =A6*EXAR+B6*TE*EXAR 

TEMP = TEMP-XJ*VREF* (0.5* VREF* ( -4 .*DUM/ ( VREF-BL I T ) +AL*DUM1 ) + 


EQU3 


)EOU3 

EQU3 

E0U3 

EUU3 

EQU3 

EQU3 

EQU3 
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l DUM+DUM 1 ) EQU3 

TEMP=TEMP+XJ* ( B I GB4*TE* ( 1 . / ( 3. * ( VREF-BL I T ) **3 ) - 1 • / ( 3. * ( VE-BL I T ) E QU3 

1**3) ) + B6* TE* (EXAL-EXAR ) /AL ) EQU3 

5 ANS3=ANS3+TFMP E0U3 

RETURN 
END 

SIB^TC ENTROP DECK, LIST 

FUNCTION ENT ( VFE , TEE ) 

C ENTROPY ENT 


COMMON A 2 , A3 » A 4 »A5,d2i63,B4 *65*C2,C3,C4, C5f XKK , R « TREF , VREF ,dL I T ♦ X J 
COMMON D4 ,Z4 ,ESS , TUFF ♦ CLOD# GEE ♦ A JA , KX,P2 1 ( 30 ) * U2 I (30 ) t M2 I ( 30 ) 

COMMON PT2 ( 30 ) , RHO ( 30 ) « TE2 I (30 ) * V E 2 l ( 30 ) * CV2 ( 30 ) *CP2 ( JO ) *GAM2 < 30 ) 

^ COMMON ACC 3 ( 30 ) .XMACH2 (J0)» RORAT ( JO > • TRAT < 30 > * PRAT (30>*PTwAT(30) 

'* COMMON TH 

COMMON VI N ( 30 ) « V2N ( 30 ) , V 1 T ( 30 ) 

COMMON B I GB4 ,B6 «A6, Al* TZ 
EQUIVALENCE (C3.G3) . (C5,XJ3 ) , (C2,D3> 

PO ( XL , AT ) = 1 ./EXP (XL* AT ) 

1 ESS 1 =Z4*AL0G ( TEE/ TREF ) +64* ( TEE- TREF ) + C4* ( TEE* *2“ TREF * *2 ) /2 • 

ESS2=ESS 1 +D4* < TEE** 3-TREF**3 ) /3 • + X J*D3*XKK/ ( VREF-BL I T ) * ( PO ( XKK , TEE 

2 )-P0 (XKK, TREF ) ) 

E S S 3 = E S S 2 + X J *G 3 * XK K / ( 2 • * ( VREF-BL I T > **2 ) * ( PO ( XKK , TEE ) -PO ( XKK , TREF ) ) 

ESS4 = ESS3 + XJ*XJ3*XKK/{4.* (VREF-BL I T ) **4 ) * (PO ( XKK , TEE ) -PO ( XKK , TREF ) 

2 ) 

ESS5 = ESS4+ (R*ALOG( (VEE-BLIT }/ (VREF-BL IT ) ) )*XJ 

ESS6 = ESS5- (B2-XKK *C2*P0 (XKK , TEE ) >* ( 1 . / ( VEE-BL IT) — 1 • / ( VREF-BL IT)) *X 
2 J 

ESS7 = ESS6- ( B3-XKK*C3*P0 ( XKK , TEE ) ) /2. * ( 1 ./ (VEE-BL I T ) **2- 1 */ (VREF-BL 
2 I T ) **2 ) *X J 

ENT = ESS 7- ( B5-XKK*C5*P0 ( XKK , TEE ) )/4 .* ( 1 ♦ / ( VEE-BL 1 T ) **4- 1 ./(VREF-BL I 
2T ) **4 ) *X J 

EXAL=EXP ( AL*VEF ) ENT 

alvr=al*vref 

EXAR = EXP ( AL VR ) 

3001 CONTINUE 

2 ENT =ENT + X J* ( B I GB4* ( 1 . / ( 3.* ( VREF-BL IT)**3)-I./(3.* (VEE-BL I T ) **3 ) ) ENT 

1 +BE * ( EXAL-EXAR )/AL) ENT 

RETURN 

END 

SIBFTC PARTIL DECK • L I ST 
FUNCTION PART (V,T ) 

C PARTIAL P W.R.T. RHO PART 

COMMON A2 , A3 , A4 ,A5, B2 i B3 ,B4 , 85 , C2 , C3 , C4 , C5 .XKK ,R * TREF , VREF , QL I T , X J 
COMMON D4 , Z4 , ESS * TUFF , CLOD , GEE , AJA , KX ,P2 I (30 ) « U2 I (30) * H2 I (30) 

COMMON PT2 ( 3 0 ) , RHO ( 30 ) ♦ TE2 I ( 30 ) , VE2 I ( 30 > « CV2 ( 30 ) , CPE ( 30 ) » GAM 2 ( 30 > 
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COMMON ACC2 ( 30 ) • XMACH2 ( 30 ) * POP AT ( 30 > < T R A T (30) * PRAT < 2 0 > * P TP AT < 30 ) 

COMMON TH 

COMMON V 1 N ( 30 ) < V2N ( 30 > * V 1 T ( 30 ) 

COMMON 8 I GB4 ,B6 « A6 % Al ♦ TZ 
EQUIVALENCE ( C3 . G3 ) * ( C5 ♦ X J3 ) ♦ ( C2 ♦ D3 ) 

PO ( XL « AT ) = I • /E XP ( XL* AT ) 

\ PART=P*T*V**2/ (V-bLlT)**2 + 4.*A4*V*-*2/(V-aLlT )**b+2.*V**2*(A2+b2* 
2T+C2*P0 ( XKK f T ) ) / ( V-bL I T )**3+3. *V**2# ( A3 +d 3*T +C3 *PQ ( XKN iT ) ) / ( V- 6L 
3JT)**4+5.*V**2* < A5+b5*T + Cb*P0 < XKK . T ) >/ < V-BL I T ) **6 
EXAL=EXP ( AL*V } PART 

2 PART=PART-(V**2)*(-4 # *3I Gb4*T/ ( V-bL I T ) **5+A6*AL*EXAL+b6#T*AL*EXAL ) PART 


RETURN 


END 

SIBFTC CEEV DECK .LIST 
FUNCTION CEV(VE.TE) 

COMMON A2.A3.A4 »A5i B2 ♦ 63 « B4 .B5.C2.C3.C4.C5.XKK.R. TREF * VREF .BLIT iXJ 
COMMON D4 .Z4 tESStTUFFfCLOD.GEE. AJA*KX.P2 I ( 30 ) * U2 I < 30 ) « H2 1(30) 
COMMON PT2(30) »RhO(30 ) . TE2 I (30) . VE2 I (30) . CV2 ( 30 ) • CP2 ( 3C ) * GAM2 ( 30 ) 
COMMON ACC2 ( 30 ) « XMACH2 ( 30 ) . ROR AT ( 30 ) . TRAT < 30 ) * PRAT ( 30 ) * PTR AT ( 30 ) 
COMMON TH 

COMMON V1N(30> »V2N( 30 > *V1T (30 ) 

COMMON B I GB4 . B6 . A6. AL • TZ 

EQU I VALENCE <C3 *G3 ) . (C5»XJ3) , ( C2 *03 ) 

PO ( XL t AT } = 1 ./EXP (XL* AT ) 

1 CE\, =Z4 + B4*TE + C4*TE**2 + D4*TE**3-XJ* ( D3*TE*XKK**2*P0 ( XK*. . TE ) )/ < VE- 
2BL I T ) -X J* (G3*TE*XKK**2*P0(XKK,TE ) ) / < 2 . * < V£-bL i T ) **2 ) -X J* <XJJ*TE* 
3X<K**2*P0 ( XKK . TE ) ) / ( 4 • * ( VE-8L I T ) **4 ) 


RETURN 
END 

FUNCT l ON PRES ( VE * TE ) 

EQUATION OF STATE 
COMMON A2 * A 3 * A4 .A5.B2.B3*b4 .B5*C2.C3.C4.C5.XKK.R.TREF. VREF . 3L ITiXJ 
COMMON D4 ,Z4 .ESS. TU FF . CLOD . GEE * A J A ♦ KX * P2 I ( 30 > * U2 I ( 30 ) * H2 I ( 30 ) 
COMMON PT2 (30) iRHO ( 30 > . TE2 1(30) *VE2I (30 ) .CV2 ( 3C ) .CP2 (30 ) *GAM2 (30 ) 
COMMON ACC2 ( 30 ) .XMACH2 ( 30 ) * RORAT ( 30 ) . TRAT (30 ) * PRAT ( 30 ) * P TR AT ( 30 ) 
COMMON TH 

COMMON V 1 N ( 30 ) . V2 N ( 30 ) * V 1 T (30) 

COMMON BIGB4.B6.A6.Al.TZ 
PO ( XL . AT ) = 1 ./EXP (XL* AT } 

POT = P0 ( XKK * TE ) 

PRES=R*TE/ (VE-BL 1 T )+ < A2+B2 *TE +C2 *POT )/ ( VE-bL I T ) **2 + ( Aj+B 3*TE 
1 +C3*P0T )/ (VE-BL IT ) ** 3+ ( A4+B I G34*TE ) / ( VE-BL I T )**4+(A5 + bio*TE + Cb*POT ) 
2/ ( VE-BL I T ) **5+ ( A6+B6*TE )*EXP(AL*VE > 

RETURN 

END 


PRES 


PRES 

PRES 

PRES 

PRES 
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FUNCTION CEP < VE * TE ) 

C CP CEP 

COMMON A 2 » A3 . A4 .A5.B2.B3.B4 « B5 » C2 tC3tC4i C 5 .XKK .R .TREF * VREF , 6 L 1 T . X J 
COMMON D 4 .Z 4 .ESS. TUFF * CLOD . GEE .AJA.KX.P2I ( 30 ) . U2 I < 30 ) * H2 I < 30 ) 

COMMON PT2(30).RHO(30)«TE2I(30). VE2 I < 30 ) . CV2 ( 30 ) « CP2 ( 30 > • GAM2 ( 30 ) 
COMMON ACC2 ( 30 ) .XMACH2 ( 30 ) • RORAT ( 30 ) * TRAT ( 30 ) » PRAT <30)»PTRAT(30) 

COMMON TH 

COMMON V 1 N ( 30 ) . V2N ( 30 ) . VI T ( 30 ) 

COMMON B I GB4 1 B 61 A 61 AL 1 TZ 

EQU I VALENCE (C3.G3) . (C5.XJ3 ) . (C2.D3) 

PO (XL. AT )=1 ./EXP(XL*AT ) 

EXAL=EXP <AL*VE > CEP 

1 DEK=R*TE/ ( VE-BL I T ) **2+2.* ( A2+B2*TE+C2*P0 ( XKK , TE ) )/ ( VE-BL I T ) **3+3. 

2*(A3+B3*TE+C3*P0<XKK.TE))/(VE-BLIT)**4+4«*A4/<VE-BLIT>**5+5.* 

3 ( A 5 + B 5 *TE + C5*P0 ( XKK . TE ) ) / ( VE-BL I T ) * *6 
DEN=DEN-4.*BIG64*TE/ (VE-BL I T ) **5+A6*AL*EXAL+B6*AL*TE *EXAL CEP 

XNU=R/( VE-BL IT )+ (B2-XKK*C2*P0 (XKK.TE ) ) / ( VE-dL I T ) **2+ ( B3-XKK*C3* 

2P0 ( XKK , TE ) )/ (VE-BL I T ) ** 3+ ( B5-XKK*C5*P0 < XKK . TE ) )/ (VE-BL I T )**5 
XNU =XNU+BIGB4/(VE-BLIT >**4+36*EXAL CEP 

XNUM = TE*X J *XNU**2 

2 CEP =XNUM/DEN 
RETURN 

END 

SIBFTC EQUA2 

SUBROUTINE EQUA2 ( A, B , DELTX « E 1 . E2 .MAX I . ANS ) 

EXTERNAL FOFX 

CALL ITR2(ANS,A,S,DElTX.F0FX,E1 * E2 . MAX I . I CODE > 

IF ( I CODE .EQ. 0 ) GO TO 100 
I F ( I COOF.EQ. 1) WRl TE ( 6 . 1 0 ) 

IF < I CODE. EQ. 2 ) WR I TE ( 6 . 1 1 ) 

IF ( I CODE . EQ. 3 ) WR I TE < 6 » 1 2 > 

I F ( I CODE. EQ. 4 ) WRI TE ( 6 . 1 3 ) 

10 FORMAT ( 25H0MAX NUMBER OF ITERATIONS) 

11 FORMAT ( 1 2H0DELX=0 OR -) 

12 FORMAT (8H0N0 ROOT) 

13 F ORM A T ( 7H0 A GT B) 

10O RETURN 

END 

FUNCTION FOFX(VX) 

C ENTROPY T 6 GET V _ F0F: 

COMMON A2.A3.A4 .A5.B2 .63. B4 . B5 .C2.C3.C4.C5. XKK .R .TREF . VREF »dL I T . X J 
COMMON D4 .Z4 .ESS. TUFF .CLOD. GEE. AJA.KX.P2 I ( 30 ) . U2 I ( 30 ) . H2 I ( 30 ) 

COMMON PT2 ( 30 ) . RHO ( 30 ) » TE2 I ( 30 ) * VE2 I (30)«CV2(30)«CP2(30).GAM2(30) 
COMMON ACC2 ( 30 ) . XMACH2 ( 30 ) .RORAT <30 ) . TRAT (30 ) .PRAT < 30 ) . PTRAT ( 30 ) 
COMMON TH 
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COr 'ION V 1 N ( 30 ) iV2N( 30 ) t V 1 T ( 30 ) 

COMMON B I GB4 ,B6 , A6, Al . TZ 
EQUIVALENCE (C3.G3) , (C5.XJ3 ) , (C2.D3) 

PO (XL. AT ) =1 ./EXP(XL*AT ) 

1 TX=TUFF 

AND = ESS- Z4*AL0G < TX/TREF >-B4* ( TX-TWEF ) -C4* < TX**2-TREF**2 ) /2 . 

ANE= AND-D4* ( TX**3-TREF**3 ) / 3 . - X J*D3*XK<* (PO ( XKK . TX >-PO ( XKK . TR£F ) ) / 

2 ( VREF-BL IT) . . . Vr .. . 

ANF = ANE- XJ*G3*XKK*(PO( XKK , TX ) -PO ( XKK . TRtF ) > / < 2.* < VRE' -BlI T > **2 ) 

ANG= ANF— X J* X J 3 * XKK * (PO (XKK.TX ) -PO ( XKK . TREE ) )/<4.*(VREF -cL I T ) **4 ) 

ANH= ANG-X J*R*ALOG ( ( VX-BLI T )/ ( VREF-BL IT) ) + (B2-XKK*C2*P0 ( XKK ,TX j )# < 1 
2./ (VX-BL I T )-l ./ ( VREF-BL I T ) ) *X J 

AN I = ANH + (B3-XKK*C3*P0 (XKK.TX) ) * < 1 ./ ( VX-8L IT ) **2-1 ./ (VREF-BL I T ) **2 ) 

2/2. *XJ 

FOFX = AN 1 + (B5-XKK*C5*P0 (XKK.TX ) )* ( 1 . / ( VX-BL I T ) **4- 1 ./( VHEF-BLl T ) **4 
2)/4.*XJ 

2 FOFX=FOFX-XJ+ (B IGB4* ( 1 ./( 3 .*( VREF-BL I T ) **3 ) - 1 ./(3.* ( VX-BLI T )**3 ) ) FOFX 

1 +B6* (EXP ( AL*VX ) -EXP ( AL*VREF ) ) /AL ) FOFX 

RETURN 

END 

SIBFTC SHOC DECK, LIST 

SUBROUTINE SHOCK ( H , VEL . P . V . T ) 

COMMON A2.A3.A4.A5.B2.B3.B4. B5 . C2 « C3.C4 .C5.XKK .R. TREF . VREF »3L I T , XJ 
COMMON D4 . Z4 . ESS . TUFF . CLOD. GEE « A J A . KX . P2 1 ( 30 > . U2 I (30 . H2 I < 30 ) 

COMMON PT2 (30 ) , RHO ( 30 ) . TE2 I < 30 > , VE2 I ( 30 ) .CV2 ( 30 > . CP2 (30 ) . GAM2 ( 30 ) 
COMMON ACC2 (30 ) « XMACH2 ( 30 ) . RORAT (30 ) ,TRAT (30 ) .PRAT ( 30 > .PTRAT ( 30 ) 

COMMON TH 

COMMON VI N ( 30 ) . V2N ( 30 ) . VI T ( 30 ) 

COMMON B I GB4 ,B6 . A6. AL . TZ 
COMMON/CLUX/TC . VC 
COMMON/DOPE/A ,B »C .PRAY 
COMMON /PRTS/ PxTV.PVTV 

COMMON /RATS/PEE. TEE, V I S2 <30 ) . PT20RT (30 ) .PI T2RT(30> . T20RT < 30 ) 

DIMENSION I PI VOT ( 2 ) , I NDEX (2 .2 ) 

DIMENSION TC <30 ) ♦ VC ( 30 ) 

D I MENS I ON HOPE ( 25 ) . BEG ( 25 ) . TEMP ( 25 ) 

DIMENSION AMAT ( 2, 2 ) .BMAT (2 ) 

EQUIVALENCE (C3 ,G3 ) , ( C5 .X J3 ) . ( C2 .D3 ) 

EXTERNAL FONC.FOND 

3 I C-hfK = 0 
ELM ./I 44 . 

K = KX 

VLN=VEL*S IN ( TH ) 

A=H+VLN**2/50082. 

B=P+ ( EL/32 • 1 74 ) *VLN**2/V 
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Program A — Continued 


C = VLN/V 

UT=VEL*COS CTH ) 

V 1 N ( K >=VLN 
V i T ( K ) = UT 

8000 FORMAT (5H1SH0C. 5X.2HA=E 15.8. 5X.2HB=E 15.8. 5X.2HC=E 15. 8. 
12HV=E1 5.8.5X .4HVEL=El 5.8 ) 

WR I TE ( 6 * 8000 ) A.B.C.V.VEL 
MAXI =30 
El = .0001 
E2 = 0 *000 1 
PRAY = TC ( K ) 

HOPE ( 1 ) = VC (< ) 

77 DO 59 IND=1 .30 
C H ( T . V ) 

HTV = FOND (HOPE ( 1 ) ) 

C P(T.V) 

PTV = FONC (HOPE ( 1 ) ) 

WR I TE (6.90 00 ) HTV.PTV 
9000 FORMAT ( 6E 1 7.8 } 

C PX(T.V) 

PXTV=FONCX (HOPE ( 1 ) ) 

C PV(T.V) 

PVTV = FONCV ( HOPE ( 1 ) ) 

WR ITE (6 . 9000 ) PXTV.PVTV 
C HX (T.V ) 

HXTV=FONDX (HOPE ( 1 ) ) 

C HV(T.V) 

H VI V = F OND V ( HOPE ( 1 ) ) 

WRITE (6.90 00 ) HXTV.HVTV 

AMAT <2.1 )=HXTV 

AMAT (2.2 )=HVTV 

AMAT (1.1 ) =PXT V 

AMAT (1.2 )=PVTV 

BM AT (2 ) = ( — H T V ) 

BM AT ( 1 ) = (-PTV) 

WR ITE ( 6 . 2000 ) BMAT.AMAT 

CALL MAT INV ( AMAT .2.3MAT . 1 . DETERM ♦ I P I VOT . I NDEX . 2 . I SC ALE ) 
WR I TE ( 6 . 2000 ) BMAT 

WR I TE ( 6.2 000 ) HOPE ( 1 ) .PRAY . PT V « HT V 
2000 FORMAT (6E2 0.8 ) 

DELt=ABS(BMAT( 1 J/PRAY ) 

DFL2 = ABS (BMAT (P )/HOPF( 1 )> 

I F ( ( DELI .LE.E1 ) • AND. (DEL2 .LE.E2J) GO TO 60 
PRA Y=PRAY+BMAT ( 1 > 

HOPE ( 1 ) =HOPE (1 ) +t3MAT ( 2 ) 
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Program A - Continued 


59 CONTINUE 

C TOO MANY ITERATIONS HERE 

WRITE (6 *2001 ) 

2001 FORMAT ( 1 5H0MAX ITERATIONS) 

RETURN 

60 CONTINUE 

WR I TE ( 6 t 1 000 ) PRAY « HOPE ( 1 ) 

1000 FORMAT ( 7H ROO T S= E 1 5 . 8 « E20 • 8 ) 

1 FORMAT ( 5H HOPE * 20X * 3HBEG • 2 0 X , 4 HTE MP/3E22 • 6 ) 

VE2 I (K ) =HOPE ( 1 ) 

TE2 I (K ) =PRAy 

P2 I C K ) = PRES ( VE2 I (K ) « TE2 I (K ) ) 

UN=C*VE2 I (K ) 

V2N (K ) = UN 

U2 I (K ) =SQRT ( (UN* *2 ) 4- (UT**2 ) ) 

CALL EQU3(TE2I <K > • VE2I (K ) ,ENTH ) 

H2I ( K ) = ENTH + 200.0 
RHO (K ) - 1 . /VE2 I (K) 

CV2 (K ) = CE V (VE2 I (K > , TF2 1(0) 

CP2 (K ) =CEP ( VE2 I (K ) , TE2 I (K ) ) + CV2 (K ) 

GAM2 (< ) =CP2 ( K )/CV2 (K ) 

ACC2 (< ) =SQRT ( ABS ( GA M2 <K ) ) *A6S ( PART ( VE2 I ( K ) « TE2 1(0)) *4633. 056 ) 

XMACH2 (O =U2 I (K )/ACC2 (K ) 

RORAT (K ) =VE2 I ( K ) / V 
TRAT (0=TE2I (K ) /T 
PRAT <K ) =P2 I (K ) /P 
RORAT (K ) =V/VE2 I (K) 

T 20RT ( K ) = TE2 I ( K ) /TEE 

V I S2 < K ) = ( 149.5* (TE21 ( K ) / T Z ) * SORT < TE2 I (K)/TZ)*1 .22/(TE2I (O/ 

1 TZ+ .22 ) ) 

VIS2 <0=VI S2 (O *0 • 672E- 7 
5 KX=KX+1 
RETURN 
END 

t> I BFTC FOND DECK 

FUNCTION FOND(VE) 

C ENTHALPY FOR NEWTON-R APHSON FOND 

COMMON A2 , A3 • A4 ♦ A5» 62 * B3 *o4 * B5 * C2 « C3 t C4 * C5 * XKK ,W ♦ TREF ♦ VRbF »dLITiXJ 
COMMON D4 . Z4 » ESS . TUFF « CLOD * GEE * A JA i <X « P2 I (30 } « U2 I (30) *H2I (30) 

COMMON PT2 (30 ) i RHO ( 30 ) « TE2 I (30)i Vt2 I (30 > « CV^d ( 3C ) * CPE ( 30 ) • GAM2 ( 30 ) 
COMMON ACC2 ( 30 ) t XMACH2 ( 30 ) « RORAT ( 30 ) * THAT ( 30 ) « PRAT ( 30 > * P TRAT ( 30 ) 

COMMON TH 

COMMON VI N (30 ), V2N ( 30 ) ♦ VI T ( 30 ) 

COMMON B I GB4 *B6«A6*Al»TZ 
COMMON/DOPE/ A , B * C ♦ PRAY 
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Program A — Continued 


EQUIVALENCE ( C3 . G3 ) ♦ ( C5 . X J3 ) . ( C2 * D3 ) 

PO (XL. AT) = 1 ./EXP (XL* AT ) 

TE = PR A V 

I H6 = Z4* (TE-TREF )+B4* ( TE**2- TREF **2 ) /2 . + C4* ( TE**3-TREF**3 ) /3 • + D4* 

2 ( TE**4-TREF **4 ) /4 • + (PO ( XKK « TE ) * < XKK*TE+ 1 • )-PO( XKK . TREF ) * ( XKK* TREF 
3+1 . ) )* (XJ*D3/( VREF-BL IT ) +X J*G3/ ( 2 . * ( VREF-BL I T ) **2 ) +X J*X J3/ (4.* 

4 ( VREF-BL I T ) **4 ) ) + XJ* ( VREF *R * ( TE-TREF ) / ( VREF-BL IT ) + VREF* < 32 * < T£> 
5TREF )+C2* (PO ( XKK ♦ TE ) -PO (XKK ♦ TREF ) ) )/< VREF-BL I T ) **2 >+XJ* < VREF* (B3* 

6 ( TE-TREF ) +C3* < PO ( XKK , TE ) -PO ( XKK , TREF ) ) ) / ( VREF-BL I T ) **3 ) + XJ* (VREF* 

7 ( B5* ( TE-TREF ) + C5* (PO ( XKK *TE)-PO(XK< ♦ TREF ) ) }/ (VREF-BL I T )**5 )- (R*TE* 
8 ( A LOG ( ( VE-BL I T ) / ( VREF-BL IT)) +3L IT/(VREF-3LlT)-BLIT/(VE-BLlT ) ) )*XJ 
H7 = H6- ( 2 •* ( A2 + B2*TE+C2*PO ( XKK * TE ) ) * ( 1 •/ (VREF-BL I T ) +BL I T/ ( 2 • * < VREF- 
2BL IT ) **2 ) - 1 */ ( VE-BL I T ) -3L I T/ ( 2 •* (VE-BL I T ) **2 ) ) > *XJ 
H8 = H7- ( 3. * ( A3 + B3*TE + C3*P0 ( XKK «TE))*(l./(2.*< VREF-BL I T ) * * 2 ) +3L IT/ (3 
2. * (VREF-BL I T )**3 )- 1 •/ (2.* ( VE-QL I T ) **2 ) -3L I T/ ( 3 • * ( VE-BL I T ) **3 ) ) ) *X J 
H9 = H8- ( 4.*A4* ( 1 ./ ( 3.* ( VREF-BL I T ) **3 )+BL I 7/ (4 .* (VREF-BL I T ) **4 ) - 1 •/ ( 
23. * ( VE-BL IT )**3)-BLIT/(4.*( VE-BL T T ) **4 ) ) ) *X J 
H 1 0 = H9- ( ( 5.* ( A5 + B5*TE + C5*P0 ( XKK «TE)))*(l*/<4.*< VREF-BL I T )**4 ) +BL I T 
2/ ( 5 • * ( VREF-BL I T ) **5 ) - 1 ./(4.*(VE-BLIT)**4)-BLIT/(5.*(VE-BLIT)**5) ) ) 
3*XJ 

H 1 1 =H1 0 + TE*R*ALOG ( ( VE-BL IT )/ ( VREF-BL I T ) > *XJ 

H 1 2 -H 1 1 -TE* ( B2- XKK* C2 *PO ( XKK , TE > ) * ( 1 • / < VE-BL I T )- l •/ ( VREF-BL IT) ) * X J 
H 1 3 = H 12- TE* ( B3-XKK*C3*PO ( XKK . TE ) >/2 •* ( 1 ♦ / ( VE-BL I T ) ** 2- 1 •/ ( VREF-BL I 


2 T ) **2 )*XJ 

HI 4 = H1 3-TE* (B5-XKK*C5*P0 ( XKK , TE ) ) /4 • * ( 1 ./ (VE-BL I T ) **4- 1 • / ( VREF-BL I 
2 T ) **4 ) *X J 

FOND=H 1 4+C**2*VE**2/ (2.*GEE*AjA )-A +200. 

EXAL=EXP ( AL*VE ) 
alvr=al*vref 


EXAR=EXP(ALVR) 

3001 CONTINUE 

TEMP = XJ*VREF* ( ( E I GB4* ( TE-TREF ) / ( VREF-BL I T > **4 ) +86*EXAR* ( TE-TREF ) )EQU3 


DUM = B I GB4* TE/ ( VF-BL I T ) * *4 EOUJ 

DUM1 = A6*EXAL+B6*TE*EXAL EQU3 

TEMP= TEMP+ X J* VE* ( 0 • 5* VE* < -4 • *DUM/ ( VE-BL I T ) + AL*DUM 1 )+ DUM+DUM1 ) EQU3 

DUM = BIGB4*TE/< VREF-BL I T ) **4 EQU3 

CUM 1 =A6*EXAR+B6*TE*EXAR EQU3 

TEMP = TE MP- X J * VREF* ( 0 . 5* VREF* ( -4 . *DUM/ ( VREF-BL I T ) + AL*DUM 1 ) + EQU3 

1DUM+DUM1 ) EQU3 

T EMP = TEMP+ X J* ( B I GB4 *TE* ( l • / ( 3* * < VREF-BL I T ) **3 ) - 1 . / ( 3 • * < VE-BL I T ) EQU3 

1**3) )+B6*TE* (EXAL-EXAR )/AL ) EQU3 

ANS 3 = AN 53 + TE MP EQU3 

2 FOfO = FOND+TEMP FOND 

RETURN 
END 
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Program A - Continued 


SIBFTC FONC DECK 

FUNCTION FONC ( VE ) 

C EQ. OF STATE FOR NEWTON-RAPHSON FONC 

COMMON A2*A3*A4*A5*B2*B3*64*B5*C2*C3*C4*C5* XKK *R * TREF ♦ VREF tBL I T * X J 
COMMON D4 * Z4 * ESS * TUFF * CLOD* GEE * A JA * KX * P2 I (30) * U2 I < 30 > * H2 I < 30 ) 

COMMON PT2 ( 30 ) * RHO ( 30 ) * T£2 I ( 30 ) * VE2 I ( 30 ) * CV2 ( 30 ) *CP2 ( 30 ) * G A M 2 ( 30 ) 

COMMON ACC2 ( 30 > * XMACH2 < 30 ) • RORAT < 30 ) * TR AT ( 30 ) * PRAT ( 3 0 ) * PTR AT ( 30 ) 

COMMON TH 

COMMON V1N(30) *V2N<30 > *V1T (30 ) 

COMMON BIGB4«B6*A6*AL.TZ 
EQUIVALENCE ( C3 * G3 ) ♦ ( C5 * X J3 ) « ( C2*D3 ) 

COMMON/DOPE/ A *B ,C * PRAY 
PO ( XL * AT ) = 1 • /EXP ( XL* AT ) 

EL = 1 ./I 44 . 

TE=PRAY 

POT = P0 ( XKK * TE > PRES 

I F ONC =R*TE/ ( VE-BL I T ) + ( A2+B2*TE+C2*POT ) / ( VE-BL I T ) **2+ ( A3+B3*TE+ FONC 

1 C3*P0T ) / (VE-BL I T >**3+ ( A4+B I GB4*TE ) / ( VE-BL I T >**4 + ( A5+B5*TE+CS*P0T > FONC 
2/ (VE-BL I T ) **5 + ( A6+B6*TE ) *EXP ( AL* VE ) + (EL /GEE >*C**2*VE-P FONC 

RETURN 
END 

FUNCTION FONDX(VE) 

C PARTIAL H W.R.T* T FONDX 

COMMON /PRTS / PXTV.PVTV 

COMMON A2*A3*A4 * A 5 * 62 * 83 * d4 *B5*C2*C3*C4* C5* XKK *R , TREF * VREF * BL I T « X J 
COMMON D4 * Z4 * ESS * TUFF * CLOD * GEE * A J A * KX « P2 I (30 > *U2 I (30 ) t H2 I ( 30 ) 

COMMON PT2 (30 ) * RHO ( 30 ) * TE2 1(30) * VE2 1(30)* CV2 ( 30 ) * CP2 ( 30 ) * GAM 2 ( 30 ) 

COMMON ACC2 ( 30 ) • XMACH2 ( 30 ) * RORAT ( 30 ) * TRAT (30 ) * PRAT ( 30 ) * PTR AT (30 ) 

COMMON TH 

COMMON V1NC30) * V2N ( 30 ) *V1T (30 ) 

COMMON BIGB4*36*A6*Al*TZ 
COMMON/DOPE/A *B *C * PRAY 

EQUIVALENCE ( C3 * G 3 ) * ( C 5 * X J3 ) * ( C2 * D3 ) 

PO (XL* AT ) = 1 ./EXP ( XL* AT ) 

1 XE=PRAY 

F0NDX=Z4 + B4*TE + C4*TE**2-XJ* ( ( D3*TE* ( XKK**2 ) *PO (XKK , TE ) )/ ( VE-BL IT)) 

1 - XU* ( ( G3*TE* ( XKK** 2 ) *P0 (XKK * TE ) )/ ( 2. 0* < VE-BL I T ) **2 ) ) 

2- XJ*< (XJ3*TE*(XKK**2)*PO(XKK,TE) ) / ( 4 • 0* ( VE-BL I T ) **4 ) ) 

3+VE*PXTV+D4*TE**3 

alvr=al*vref 
EXAR=EXP ( ALVR ) 

3001 CONTINUE 

EXAL=EXP ( AL*VE ) 

DUM = B I GB4/ ( VE -BLI T ) **4 
DUM 1 =B6*E X AL 


FONDX 

FONDX 

FONDX 
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Program A — Continued 


c 


TEMP=X J* VE *(DUM + DUM1+VE*(-4.*DUM/(VE-BLIT)+AL*DUM1 > > 

TEMP=TEMP-0.5*XJ*VREF*VREF*C-4.*BIGB4/<VREF-BLIT)**5+66*AL*EXAR> 

1 FONDX=TEMP+FONDX+XJ* ( B I GB4* C 1 • / (3#* ( VREF-BL I T ) **3 ) l • / < * 

1 BL I T ) **3 > ) + (56/AL )* (EXAL-EXAR ) ) 

RETURN 

FND 

FUNCTION FONCX(VE) 

PARTIAL OF P ( V 4 T ) W • R * T • T 

COMMON/PR TS/ PXTV*PVTV _____ wQ „ TT 

COMMON A2.A3.A4.A5,B2.B3.34.S5.C2.C3.C4,C5.XKK,R,TREF.VREF.i3LI < 


COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 


A2.A3,A4.A5,B2.83.34.S5.C2.C3.C4,C5.XKK,R,TREF.VREF,0LIT.XJ 

D4 .Z4 iESSi TUFF . CLOD .GEE.AJA.KX.P2 I < 30 ) » U2 I (30). M 2 1(30) 

PT2 ( 30 ) » RhO ( 30 ) . TE2 I ( 30) . VE2 I ( 30 ) « CV2 ( 30 ) . CP2 ( 30 ) «GAM2(30> 
ACC2 ( 30 > . XMACH2 < 30 > . RORAT ( 30 ) . TRAT ( 30 ) . PRA T ( 30 ) « PTRAT ( 30 ) 


COMMON V 1 N ( 30 ) « V2N ( 30 ) t VI T ( 3C ) 

COMMON B I GB4 4 36 4 A6 4 AL 4 TZ 
EQUIVALENCE ( C3 4 G3 I 4 (C54XJ3) 4 { C2 4 D3 ) 

COMMON /DOPE /A »B 4 C 4 PRAV 
PO (XL 4 AT ) = 1 «/EXP ( XL* AT ) 

EL = 1 * / 1 4 4 • 

1 TF = PR A V 

FONCX=R /(VE-BLlT)+(82-C2*XKK*P0(XKK,TE> ) / ( VE-BL I T ) * *2 
1 + ( B3-C3*XKK*P0 ( XKK, TE > ) / ( VE-BL I T ) *#3 + ( B5-C5* XKK*PO ( XKK . TE ) ) 

2/ ( VE-RL I T )**5 P 

2 FONCX=FONCX+B I GB4/ ( VF-BLI T ) **4 +B6 »EXP ( AL* VE ) 

RFTURN 

FND 

FUNCTION FONDV(VE) F 

PARTIAL H W.R.T. V 

COMMON/PRTS/ PX TV.PVTV ____ al IT V i 

COMMON A2,A3.A4.A5,B2.S3.34,B5.C2,C3.C4.C5.XKK.R,TREF.VREF.3LIT.XJ 

COMMON D4,Z4,ESS,TUFF,CL0D.GEE.AJA.KX.P21 ( 30 ) .U2 I ( 30 > . H2 I < 30 ) 

COMMON PT2 ( 30 ) .RHO (30 ) , TE2 I ( 30 ) . VE2 I ( 30 ) . CV2 (30) .CP2 (30) .GAM2(30) 
COMMON ACC2(30) .XMACH2( 30) .RORAT (30) .TRAT (30) .PRAT (30 ) .PTRAT (30) 


COMMON TH 

COMMON VI N ( 30 ) 4 V2N ( 3C ) 4 V 1 T ( 30 ) 

COMMON B I GB4 * B6 4 A6 4 AL 4 T Z 
COMMON /DOPE/ A , 3 4 C 4 PRAY 

EQUIVALENCE ( C 3 4 G3 ) 4 (C54XJ3) 4 ( C2 4 D3 ) 

PO ( XL 4 AT ) = 1 #/EXP ( XL*AT ) 

] te=ppay 

FONDV=VE*PVTV +TE*PXTV +<C**2)*VE/< GEE*AJA ) 
FXAL=EXP < AL*VE ) 

DUM=B I GB4*TE/ ( VE-BL I T ) **4 
CUM 1 = A6*EXAL 


FONDV 

FONDV 

FONDV 
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Program A - Continued 


DUM?=B6*TE*EXA|_ FONDV 

TEMP = 2.0*XJ* VE* ( -4. *DUM/ ( VE-BL I T ) +AL* ( DUN 1 +DUM2 ) >+X J* <2. 0* ( OU«i+ FONDV 

1 DUM2 )+DUM 1 ) FONDV 

2 F0NDV=F0NDV+TEMP+0*5*X J*VE*VE* ( 20.*DUM/ (VE-BL I T ) **2 + < AL**2 ) * FONDV 

1 (DUM1 +DUM? ) ) FONDV 


RETURN 

ENL 

SIBFTC FONCV LIST 

FUNCTION FONCV(VE) 

C PARTIAL OF P(VtT) W • R # T • V FONCV 

COMMON /PRTS/ PXTV»PVTV 

COMMON A2, A3 , A4 , A5,52 ,53 *B4 , B5 * C2 , C3 , C4 , C5,XKK ,R , TREF , VREF *3LI T , X J 
COMMON D4 , Z4 ,ESS , TUFF , CLOD , GEE , A J A «<X» P2 I ( 30 > , U2 I ( 30 > , H2 I ( 30 ) 

COMMON PT2 (30 ) iRHO < 30 > , TE2 I ( 30 ) * VE2 I ( 30 ) * CV2 ( 30 ) ,CP2 ( 30 ) * G AM2 ( 30 ) 

COMMON ACC2 ( 30 ) ,XMACH2 ( 30 ) * RORAT (30 ) * TRAT (30 ) * PRAT ( 30 ) * PTR AT ( 30 ) 

COMMON TH 

COMMON VI N (30 ) , V2N ( 30 ) * V 1 T < 30 ) 

COMMON RlG34*R6*A6fALtTZ 
EQUIVALENCE ( C3 , G3 ) , (C5,XJ3> * (C2,D3) 

COMMON/DOPE /A , B , C t PRAY 
PC (XL* AT ) = 1 ./EXP (XL* AT ) 

EL = 1 • / 1 44 » 

TE=PRAY 

POT = PO ( XKK * TE ) FONCV 

1 F ONC V = ( -R*TE )/ (VE-BL I T ) **2-2 *0* ( A2+B2*TE+C2*POT )/ (VE-bL I T ) **3 FONCV 

1 - 3.0*(A3+B3*TE+C3*POT >/ ( VE-BL I T >**4-4# 0* ( A4 +B I GB4*TE ) / ( VE-BL I T ) **5F0NCV 

2- 5.0* ( A5 + B5*TF + C5*P0T )/ (VE-BL I T ) **6+AL* ( A6+B6*TE ) *E XP ( AL*VE ) FONCV 

3+ ( EL/GEE ) *C**2 FONCV 

RFTURN 

END 

$DATA A2*-2. 1 62959, A3 =4.04 057E -3 , A4 = l *92 1 0 72 E-4 , A 5 = -4 . 4 8 1 049E -6 , 

A6= 5 • 83882E7 , B I GB4 = -3 . 9 1 3263E-7 , AL=“6*61 199E2, B6"-9 • 26392E4 , 

XKK = *976E— 2,B2=2. 1 351 14E-3,B3=1 • 2828 1 8E-5 , B4 = 3 . 00892783E-4 * 

B5 = 9 . 0623 18E-9,XJ=* 185053, 

C2 -- 18.941 131 , C3=0. 539776, C4=-l • 3023744 1 E-7 , C5=-4 . 836678E- 5 , 

R- 0 • 1219336, BL I T = 0 • 15C-2,D4=1 .96802394E-1 1 ,PEE=0.1026, 

TEE = 260 , GEE = 32. 1 74 , A JA = 778 • 3 , A = • 008,5=12000* DELTX = 2 , 

TRrF=F20 ,PREE= 1 ,TZ=409. 5,74=1 • 904SRC84E-2 , 

TE I =5*532 , El = .0000 1 , E2= .00 00 1 tMAXI =30 ♦ NUMB = 5 , I NT =0i 
IFNO=30 , 

THFT A = 90 , 

PEE = 2000 , 

TEE = 1 460 , 

MAX 1=30, 
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Program A - Concluded 


T £ I = 1 4 4 0 * 1420* 1 380 * 1 340 * 1 3 i 0 ♦ l 280 * 1 240 • 1 190*1140*1090* 1040*985*940 * 900*360* 
820* 760* 740* 700 * 660* 620* 580* 540* 500* 470* 440*410*380*360*340* 

VC(1 ) = 0 • 6 * 1 .02*1 .5*2. 0*2. 5*3. 0*3. 75*9. 0*14. 0*20 #0*25. 0*31 .0*50.0*80.0*100.0* 
150. 0*200*0*250 .0*32 5.0* 

TC { 1 ) = 1 O* 1 447S 
SDATA 
P c_cr = 2000* 

TEE = 1 260 * 

M Ax I = 30 * 

TF I = 1 245 * 1 230 * 1 200 * 1 160*1 120*1080*1020*1000* 960 * 920*880* 840 * 800 * 750 • 700 * 650 * 
600* 575* 550* 525* 500* 475* 450* 425*400*375*350*325*300*275* 

VC(1 )=0.4* 0.58*0.89*1 .39*2.2*3*6*4.5*6.2*8.0*10.0*14. 0*19*0* 25 .0*34 .0*47.0* 
65.0*91 .0*128.0* 190.0* 

TC < 1 ) = 1 9* 1 240S 
SDATA 
PEE=2C00 * 

T FF= 1 060 * 

MAX 1 =30 * 

TE 1 = 1030* 1035* 1020* 1000*965*930*980*830* 785 * 745 « 700*675 *640*600*575*550* 

525* 500 * 4 75* 4 50*425*400*375*350*325*300*275*250* 225*200* 

VC(1 >=0.3*0.44*0 .689 *0.875 t 1.12*1*44*1. 83 *2.46*3. 26 *4. 34 *5. 84 *7.91*10. 83* 
15.0*20.9,30.0*42.0*70.0*91.0* 

TC ( 1 ) = 1 034 , l 034 , 1 7* 1 033S 
SDATA 
PFE=2000 * 

TFE = 360 * 

MAX I =30 * 

T5I=845*830*900* 76 0*720*680* 64 0*600* 560* 520 *500*490*450*400*375* 350 * 325 * 

300 *275*250* 230*2 15*200* 190* 180* 170* 160* 150* 140*130* 

VC ( 1 ) =0.298 * 0 .307, 0 .465 * 0 .793* 1 .05* 1 .42 , 1 .92 *2*65 * 3.71 * 5.26 * 7.59*9.00 * 1 1 . 1 7* 
13.00*15.48* 18*34*21 .eo* 26 . 1 *31.5*38.3* 

TC(1 )= 82 0*819, 917, 816* 815*915*14 *8 14S 
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DATES 

Ob SCI INITIATED OS S4 COMPLETED 
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TIME AM D COST FOR DEVEl OPMENT 


L-, SAMANMONTHS W 64 MAC HINT |0S fc9 COMPUTER TYPE W 74 TOTAL COST 
HOURS 1 DOLLARS' 


6C €1 6? 63 1 64 6S 66 67 68 


ELITE MARGIN j PICA MARGIN 


o ABSTRACT L_ 

D NUMBER O — “ . , j -i r • 

06 The relation s across a n ormal s hock wave in e qul li 4- 

08 thermodynamic properties of one -dimensional flow- _The.-inpnts 

_ are_^a^e ^a/fcic^eii^e r ^ure^ jtatiejEfiil .projectile. 

10 1 velocit y, * f 




APPENDIX C - Continued 


Program B - Continued 


PROGRAM BL ISTC ( I NP'jT « OUTPUT * TAPE5= I NPUT * T APE 6 = OUTPUT ) 

COMMON A2 * A3 * 44 * A5*B2*B3*54 *35 * C2 * C3*C4*C5* XKK iR» TREF * VREF *3L I T * X J 

COMMON C4 * Z4 * ESS * TUFF * CLOD * GEE * A J A « KX * P2 I < 30 ) * U2 I <30 ) * ^2 I ( 30 ) 

COMMON PT2(30)*RHO(30)*TE2I (30)*VE2l (30)* CV2 ( 30 ) *CP2 (30 ) *GAM2(30 ) 

COMMON ACC2 ( 30 ) * XMACH2 (30) * RORAT (30 )*TRAT (30) * PRAT (30 ) *PTRAT (30 ) 

COMMON °IG34*36*A6*Al 

COMMON /CLUX/ TC*VC*EI*E2 

COMMON/ DELS ° / VO C D 00 « ZULU 

EXTERNAL FOP 

0 I MENS ION TC ( 30 ) * VC ( 30 ) 

EQUIVALENCE <C3 *G3 ) * ( C5 *XJ3 ) * (C2 * 03 ) 

0 I MENS I ON TEI (30)*VEI (30)*PEI ( 30 ) • V I (30)*GAM(30)*SPEED(30) 

C I MENS I ONXMACH ( 30 ) * CP ( 30 ) * CV ( 30 ) *H < 30 ) ♦ V I S ( 30 ) 

NAMELIST /DATA/A2*A3*A4*A5*A6*52*33*5IGQ4*35*B6*B4*AL* 

1 C2 * C3 * C4 * C5 * XKK * X J*R* BL I T • D4 * PEE * TEE • GEE* A J A * NUMB * I NT * Z4 * A #B * 
2DELTX* I END* TEI *F1 ♦ E2 * MAX I * TREF * PREF * TZ * TC * VC ♦ THETA , V I *PEI * A1 
PO ( XL * AT ) = 1 */EXP ( XL* AT ) 

99 READ ( 5 * DAT A ) 

WR I TE ( 6 * DATA ) 

KX= 1 

CEL XI = DEL TX 
V OODOO = TREF 
ZULU=PREF 

CALL I TR2 (SP*A*B*DELTX* FOP « E 1 *E2*MAXl * I CODE ) 

VREF = SP 

WRITE (6*91 ) VREF 
91 FORMAT ( 5H0VREF/EI 4.5 ) 

CO 111 N= 1 ♦ I END 
V00D00 = TE I (N ) 

ZULU=PE I (N ) 

DELTX = DELX 1 

CALL I TR2 ( SP * A 1 * 5 * DELT X * FOP * E 1 * E2 * MAX I * I CODE > 

WR I TE ( 6 * 1 00 ) I CODE 
1 00 FORMAT ( 7H0 I CODE= I 5 ) 

VE I (N ) =SP 

WR I TE ( 6 *9 ) VE I ( N ) 

9 FORMAT ( E20 *8 ) 

CALL EQU3 (TEI ( N ) * VE I (N ) * A NS 3 ) 

H ( N ) = ANS3+200. 

HO=H(N)-MVl ( N ) * V I (N) ) / ( 2 • *GEE + A JA ) 

CV (N ) = CEV ( VE I (N ) * TE I ( N ) ) 
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APPENDIX C - Continued 


Program B - Continued 


CP(N)=CEP(VEI < N ) • TE I < N ) ) +CV ( N > 

GAM(N)=CP(N)/CV ( N ) 

SPEED < N ) =SQPT ( GAM (N )*PART (VEI (N ) «T£I (N ) ) *4633 #056 ) 

XMACH (N ) = V I ( N ) /SPEED < N ) 

VIS(N)=149*5*(TEI (N)/TZ )*SQRT (TEI(N)/TZ>* 1 • 22/ ( TE I < N > /TZ+ • 22 ) 

IF (XMACH (N ) .GT.2. .AND. 1 t£ND. LT . 1 5 ) C ALL SHOCK ( H ( N ) . V I ( N ) , PE I ( N ) , VE I 
2 ( N ) * TE I <N ) ) 

1 1 1 CONTINUE 
KQ = KX— 1 

WRITE(6*52)VPE* PEE %H0«ESS*GAM0«TREE# PREF 
52 FORMAT (3H1 VO. 15X.2HP0. lbX.2HH0. 1 5X . 2HS0 . 1 5X , 5HGAMMA , 15X.4HTREF, 1 5X 
2 , 4HPREF// 7E 16*6 ) 

WRI TE(6t56) (TEI (N) • VEI (N) *PEI (N ) «H (N ) *XMACH (N ) «N= 1 • I END > 

^56 FORMAT ( 1 2 HO TEMPER AT URE * 1 5X * 6HV0LUME * 1 5X * 3HPRESSURE ♦ 1 5X * 8HENTHALP V • 
2 I 5X * 7HM ACH N0///(5E20*6 ) ) 

WR I TE ( 6 * 54 ) (VI ( J > • CP ( J ) • CV ( J ) * GAM ( J ) i SPEED (J)*VIS(J) »J=1 « I END ) 

54 FORMAT ( 9H 1 VELOCITY *20X*2HCPt20X*2HCV*20X.5HGAMM A, 15X* 1 2H ACCEL.ER A T I 
20N* 15X*9HVlSC0SITY///(6E21 .6) ) 

I F ( I END * L T • 1 5 ) WRI TE <6« 57) ( VE2 1 <K ) « TE2 I <K ) *P2 I ( K ) * U2 I. ( K ) « XMACH2 ( K ) • 
2H2I (K) *K=1 *KQ) 

I F ( IEND.LT *15) WR I TE (6*56) ( RHO ( K ) # C V2 ( K ) * CP2 (K ) * G AM2 (K ) » 

2 A C C 2 (K ) *K= 1 tKO ) 

I F ( I END* LT *15) WRITE (6*59) (RORAT (K *) « THAT (K) * PRAT (K ) , <= 1 *<Q ) 

57 FORMAT (8H1 VOLUMES «20X * 1 2HTEMPER ATURE2 * 1 2X * 9HPWESSURE2 * 12X*9HVEL0C I 

2T Y2 * 1 2X * 8HMACH NO 2 ♦ l 2X • 9HENTHALPV2/// ( 6E 2 1 • 5 ) ) 

5B FORMAT ( 5H 1 RH02 i20X» 3HCV2 • 1 6 X * 3HCP2 * 1 6X « 6HGAMM2 * I 6X * 1 4HSPEED OF SOU 

1 ND/// ( 5E2 1 *5 ) ) 

crc> FORMAT ( 6H 1 V2/V 1 *2 OX * 5HT2/T 1 *16X*5HP2/Pl///< 3E2 1*5)) 

GO TO 99 
END 

SlBFTC SHOC DECK, LIST 

SUBROUTINE SHOCK ( H * VEL « p « V * T ) 

COMMON A 2 * A 3 * A 4 • A5« 32 « 33*b4 • 95 ♦ C2 • C3 • C4 . C5* XKK ,R * TREF . VREF ^Ll T ,XJ 
COMMON D4 * Z4 * ESS t TUFF * CLOD * GEE » AJA«KX *P2 I ( 30 ) « U2 I ( 30 )« H2 I < 30 > 
COMMON PT2 (30 ) *RHO( 30 ) * TE2 1(30) * VE2 1(30) ♦ CV2 < 3C ) ,CP2 (30 > * GAM2 ( 30 ) 
COMMON ACC2 ( 30 ) t XMACH2 ( 30 ) * RORAT ( 30 ) ,TkAT ( 30) ♦ PRAT ( 30 ) * PTR AT ( 30 ) 
COMMON e I G34 *56 * A 6 * Al 
COMMON /CL'JX/ TC*VC*E1*E2 
COMMON /DOPE /A * 9 * C * PRAY 
COMMON/PRTS/ P XTV * PVT V 

common /Rats/pee * tee * pt 20RT (30) *PiT2RT ( 30 ) * T20RT ( 30 ) 

D I MENS ION I P I VOT ( 2 ) * INDEX (2*2) 

D I MENS ION TC (30 ) * VC ( 30 ) 

D I MENS! ON HOPE (25 ) * BEG ( 25 ) * TEMP < 25 ) 

DIMENSION AMAT ( 2 * 2 ) *BMAT (2 ) 
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APPENDIX C - Continued 


Program B - Continued 


FQL I VALENCE (C3 <G3) « < C5 *XJ3 ) , (C2*D3 ) 

EXTERNAL F ONC « c OND 

I chf:< = c 

EL = 1 • /I 44 ♦ 

K = KX 

A=H+vEL**2/(50082. ) 

B = P+ (EL/32. 1 74 >*VEL**2/V 
C=VEL/V 

60 0 3 FORMAT ( 5H 1 SHOC , 5X , 2HA = E 1 b • 8 * 5X ♦ 2HB =E 1 5 • 8 , 5X i 2 HC = E 1 5 • 8 ♦ 
12HV= FI 5.8 »5X*4HVEL=E1 5.8 ) 

WR I TE (6 . 5000 ) A .9«CtV* VEL 
MAXI =30 

pray = tc (< ) 

HOPE ( 1 ) = VC (K ) 

77 CO 59 I ND = 1 . 30 
C H ( T . V ) 

HTV = FOND (HOPE ( 1 ) ) 

C P ( T • V ) 

PT V =F ONC ( HOPE ( 1 ) ) 

WRITE (6. 9000 ) HTV.PTV 
9000 FORMAT < 6E 1 7*8) 

C PX(TtV) 

PX T V = F ONCX (HOPE ( 1 ) ) 

C PV(T,V) 

PVTV = FONCV (HOPE ( l ) ) 

WR I TE ( 6 t 90 00 ) PXTV , PVTV 
C HX(T.V) 

HXTV = FONDX (HOPE ( I ) ) 

C HV (T .V ) 

HVTV=FONDV (HOPE ( I } ) 

WRITE (6 *9000 ) HXTV.HVTV 

AMAT ( 2# 1 ) =HXTV 

AMAT (2.2 } = H V T V 

AMAT (1.1) =PXTV 

AMAT ( 1 i 2 ) = P V T V 

BM AT ( 2 ) = ( -HTV ) 

9MAT (1 )= (-PTV) 

WRITF(6.2000) B MAT, AMAT 

CALL MAT INV( AMAT , 2 ♦ BM AT , 1 , DETERM , I P I VOT * I NDEX ♦ 2 . I SCALE) 
WR 1 rF(6,2C00 ) BMAT 

WR I TF (6 . 20 00 ) HOPE ( 1 ) , P RA Y , PT V « HT V 
POC'*' FORMAT ( SE2 0.8 ) 

CFL1 =ABS (BMAT ( 1 )/PRAy ) 

DFL2 = A3S (BMAT ( 2 )/HOPE ( I ) ) 

IF(( DELI .LE.E1 ) . AND. (DEL2 . LE • E2 ) ) GO TO 60 
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Program B - Continued 


PRAY=PRAy+BMAT ( 1 ) 

HOPF ( 1 ) = HOPF U )+UMAT ( 2 ) 

39 CONTINUE 

C TOO MANY ITERATIONS HERE 

WR I TE ( 6 , 200 1 ) 

2001 FORMAT < 1 5H0MAX ITERATIONS) 

RETURN 
60 CONTINUE 

WRITE(6« 1000 ) PRAY, HOPE ( 1 ) 

1000 FORMAT ( 7H ROOTS=E 1 5 . 8 , E20 . 8 ) 

1 FORMAT ( 5H HOPE « 20X , 3HBEG, 20X « 4HTEMP/3E22 .6 ) 

VE2 1(0 =HOPE ( 1 ) 

TE2I (K)=PRAY 

P2 1 (K >=PRES < VE2 I (K) ,TE2I (<) ) 

U2I ( K ) = C*VE2 I <<) 

CALL E0U3 (TE2 I <K ) , VE2 I ( O ,ENTH ) 

H2I (K)=ENTH + 200.0 

RHO (K ) = I . /VE2 1(0 

CV2 (O =CEV ( VE2 1(0, TE2 1(0) 

CP2 (K ) =CEP ( VE2I (O , TE2 I (O )+CV2 (< ) 

GAM2 (<)=CP2 (0/CV2 ( <) 

ACC2 ( K ) =SQRT (ABS(GAM2(K) ) * ABS ( PART ( VE2 1 (O.TE2I (O ) ) *4633.056 ) 
XMACH2 ( K ) =U2 I (< )/ACC2 (K ) 

RORAT ( < ) = VE2 I (K )/V 
TRAT ( K )=TE2I (K )/T 
PRAT (K ) =P2 1 CO/P 
RORAT (K ) =V/VE2 I (< ) 

T2CRT(K ) =TE2 I (K > /TEE 

KX=KX+1 

RETURN 

END 

SIBFTC FOND DECK 

FUNCTION FOND(VE) 

COMMON A2,A3,A4,A5, 02 ,B3,B4 ,B5, C2 ,C3,C4,C5,XKK,R, TREF , VREF ,6L I T , X J 
COMMON D4 , Z4 ♦ ESS, TUFF , CLOD, GEE, AJA , KX , P2 I ( 30 ) « U2 I ( 30 ) » H2 I ( 30 ) 
COMMON PT2 ( 30 > , RHO ( 30 ) , TE2 I ( 30 ) « VE2 I ( 30 ) , CV2 ( 30 ) , CP2 ( 30 ) . GAM 2 ( 30 ) 
COMMON ACC2 ( 30 ) ,XMACH2 ( 30 ) .RORAT (30 ) , TRAT (30 ) , PRAT (30 ) , PTH AT (30 ) 
COMMON 3 I GB4 ,36 , A6. AL 
COMMON/DOPE/A , B ,C .PRAY 

EQU I VALENCE (C3»G3)» <C5,XJ3), ( C2 , D3 ) 

PO (XL , AT ) = 1 ./EXP ( XL*AT ) 

TE=PRAY 

H6 =Z4* (TE-TREF )+S4* ( TE**2-TREF**2 ) /2 . +C4* ( T£**3-TREF **3 )/3.+D4* 

2 ( TE**4-TREF**4 ) /4 . + < PO (XK<»TE)*(XKK*TE+1 • ) -PO ( XKK , TREF ) * (XK<*TREF 
3+1 . ) )* <XJ*D3/< VREF-BLl T ) + X J*G3/ ( 2 . * ( VREF-BL I T ) **2 ) + X J*XU3/ (4 . * 
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APPENDIX C - Continued 


S IBFT 

C 


Program B - Continued 


4 (VREF-BL IT ) **4 ) ) + XJ* { VREF*R* ( TE-TREF )/ < VREF-BL I T )+vREF* (82* (TE- 
5TREF )+C2* (PO (XKK , TE ) -PO (XKK , TREF ) ) )/ ( VREF-BL I T ) **2 ) + XJ* < VREF* ( 63* 

6 (TE-TREF ) +C3* (PO (XKK, TE >-PO ( XKK , TREF ) ) )/ ( VREF-BL IT ) **3 ) + XJ* ( VREF* 
7(B5* (TE-TREF )+C5* ( PO ( XKK * T£ )-PO (XKK, TREF > ) )/( VREF-BL IT )**5>- (R*TE* 
8 ( A LOG ( ( VE-BL IT )/ (VREF-BL I T ) ) +BL I T/ ( VREF -BL I T ) -BL I T/ ( VE-BLI T ) ) ) *XJ 
H7 = H6- ( 2* * ( A2 + B2*TEFC2*PO ( XKK , TE ) ) * ( 1 • / ( VREF-BL IT)+BLIT/(2.*(VREF- 
2BL I T ) **2 )- 1 •/( VE-BL I T ) -BL I T/ ( 2 . * (VE-BL I T ) **2 ) ) ) *X J 
H8 = H7- ( 3 • * ( A3 + S3*TE + C3*PO (XKK «TE))*<l./(2.*( VREF-BL I T )**2 >+3L I T/ ( 3 
2. * ( VREF-BL I T )**3 > -1 ./( 2 .* (VE-BL I T > **2 )-BL I T/ ( 3#* < VE-BL I T ) * *3 ) ) ) *X J 
H9 = H8- ( 4 • *A4* ( 1 * / ( 3 * * ( VREF-BL I T ) **3 5+BL I T/ ( 4 • * < VREF-BL I T ) **4 ) - 1 ./ ( 
23** (VE-BL IT )**3)-BLIT/(4«* ( VE-BL I T ) **4 ) ) ) *X J 
HI 0 - H9— ( (5 •* ( A5 + B5*TE + C6*P0 ( XKK . TE > ) ) * ( 1 • / ( 4 • * < VREF-BL I T )**4 )+BL I T 
2/ ( 5 • * ( VREF-BL I T ) **5 ) — 1 . / ( 4 . * ( VE-BL I T ) **4 ) -BL I T/ ( 5 • * ( VE-BL I T ) **5 ) ) ) 
3*X J 

HI 1 =H1 0 + TE*R*ALOG( (VE-BL IT )/ (VREF-BL IT ) >*XJ 

H 1 =H 1 i - TE* < 3 2 - XKK * C 2 * p 0 < XKK , TE > )* ( 1 • / < VE-bL I T ) - 1 • / ( VREF-BL IT) ) *XJ 
HI 3 = H 1 2 — TE * ( 03-XKK*C3*PO ( XKK , TE ) ) /2 •*(!•/ (VE-BL I T ) **2-1 • /< VREF-BL I 
2T ) **2 ) *XJ 

HI 4 = H1 3-TE* (B5-XKK*CF*P0(XKK,TE ) )/4.* ( 1 ./ (VE-BL I T ) **4-1 ./(VREF-BL I 


2 T ) ** 4 )*XJ 

FOND = Hl 4 + C**2*VE**2/ ( 2 . *GEE*A JA ) -A +200* 

EXAL=EXP ( AL*VE ) EQU3 

EXAR=EXP ( AL*VREF ) EQU3 

TEMP = XJ*VREF* ( (3 1 GB4* ( TE-TREF )/ ( VREF-BL I T ) **4 ) +B6*EXAR* ( TE-TREF ) )EQU3 
DUM = B I GB4*TE/ ( VE-BL I T ) **4 * EQU3 

DUM1 =A6*FXAL+B6*TE*EXAL EQU3 

TEMP=TEMP+XJ*VE* (0.5*VE* (-4.*0UM/ ( VE-BL I T )+AL*DUMl ) + DUM+DUM1 ) EQU3 

CUM = B I GB4* TE/ ( VREF-BL I T )**4 EQU3 

DUM1 =A6*FXAR+B6*T E*EXAR EQU3 

TEMP = TEMP- XJ* VREF* (0.5* VREF* ( -4 .*DUM/ ( VREF-BL I T )+AL*DUM 1 ) + EQU3 

I DUM+0UM1 ) EQU3 

TEMP=TEMP+XJ*(BIGB4*TE* ( 1 ./(3.* (VREF-BL I T)**3)-l./(3.* (VE-BL IT ) EQU3 
1**3) ) +B6* TE* (EXAL-EXAR ) /AL ) EQU3 


FOND=FOMD+TEMP 

RFTURN 

END 

C FONC DECK 

FUNCTION FONC(VF) 

EO. OF STATE FOR NE WTON-R APHSON FONC 

COMMON A2 , A3 « A4 ,A5. B2 *83,64 , B5i C2 * C3 , C4 ,C5i XKK ,R » TREF , VREF ,3L I T , XJ 
COMMON D4 , Z4 ,ESS« TUFF ♦ CLOD % GEE ,AJA»KX,P2I (30 ) *U2 I ( 30 ) • H2 I ( 30 ) 

COMMON PT2 ( 30 ) * RHO ( 30 ) » TE2 1(30) « VE2 I ( 30 ) « CV2 ( 30 )• CP 2 ( 30 ) tGAM2 (30 ) 
COMMON ACC2 ( 30 ) * X MACH 2 ( 30 ) * RORAT ( 30 ) * TR A T ( 30 ) , PRAT (3o ) * P TR A T ( 30 ) 

COMMON BIGB4 ,06 *A6, AL 

EQU I VALENCE { C3 *G3 ) * ( C5 *XJ3 ) t ( C2 « D3 ) 
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APPENDIX C - Continued 
Program B - Continued 


COMMON/DOPE/A .BtCi PRAY 
PC (XL. AT )=1 ./EXP(XL*AT ) 

EL = 1 ./I 44 . 

II T=PoIxKK,TE) FONCV 

FONC =R* TE/ ( VE-BL I T ) + < A2+B2*TE+C2*P0T )/ (VE-BL I T )**2 + ( A3+B3*TE+ FONC 

lC3*P0T >/ ( VE-BLI T>**3+ (A4 + BI GB4*TE )/ (VE-BL IT ) **4 + ( A5+B5*TE+C5*P0T ) FONC 
2/ (VE-BL I T )**5 + ( A6+B6*TE ) *EXP ( AL*VE ) + ( EL /GEE ) *C**2*VE-B FONC 

RETURN 
END 

FUNCTION FOP(V) 

EON OE STATE-USED TO FIND V AS A FUNCTION OF P+T 

COMMON A2.A3.A4.A5. B2 .63.B4 .B5.C2.C3.C4.C5. XKK .R.TREF.VREF.3LIT.XJ 
COMMON D4 , Z4 .ESS. TUFF . CLOD ♦ GEE . A J A . KX . P2 I < 30 > . U2 I < 30 ) . B2 I < 30 > 

COMMON PT2 <30 ) ,RHO( 30 > . TE2 1(30) , VE2 I ( 30 ) . CV2 < 30 ) , CP2 < 30 > . GAM2 ( 30 ) 
COMMON ACC2 ( 30 ) . XMACH2 (30).RORAT(30).TRAT(30). PRAT ( 30 > »PTR AT ( 30 > 

COMMON B I GB4 .36 . A6 . AL 
PO (XL . AT ) = 1 ./EXP < XL* AT ) 

COMMON/DESP/VOODOO . ZULU 
T= VOODOO 

P=ZULU FQp 

POT=PO (XKK . T ) rno 

FOP = P-(R*T )/(V-BLIT)-(A2 + B2*T + C2*P0T ) / ( V-BL I T ) **2~ ( A3+B3*T + FOP 

1 C3*P0T ) / ( V-BL I T ) **3- ( A4+b I GB4*T ) / ( V-BL I T ) **4- ( A5+B5* T + C5*POT ) / FOP 

2 ( V-BL IT)**5-(A6-B6*T)*EXP(AL*V) FOP 

RETURN 
END 

FUNCTION CEV(VE.TE) w0trI _ . T VI 

COMMON A2.A3.A4 .A5.B2 .B3.B4 .B5.C2 .C3.C4.C5. XKK ,R , TREF « VREF .BLIT.XJ 
COMMON D4 ,Z4 ,ESS. TUFF , CLOD . GEE .AJA.KX» D 2I ( 30 ) . U2 I ( 30 > . H2 I ( 30 ) 

COMMON PT2 (30) .RHO(30 ) .TE2I (30 ) . VE2 I (30 ) . CV2 ( 30 ) . CP2 ( 30 > .GA^2(30 ) . 

COMMON ACC2 < 30 ) . XMACH2 ( 30 ) . RORAT ( 30 ) , TRAT ( 30 ) .PRAT ( 30 > . PTRAT ( 30 ) 

COMMON BIGB4.B6.A6.AL 

EQUIVALENCE (C3.G3) . (C5.XJ3 ) . (C2.D3) 

PO ( XL . at ) = 1 ./EXP ( XL* AT ) 

CEV=Z4+B4*TE+C4*TE**2+D4*TE**3-XJ* ( D3*TE*X<K**2*PO < XKK , TE ) )/ ( VE 

2BL IT )-XJ* <G3»Tt*X!«**2*P0<X<K.T£ ) ) / ( 2 . * < VE-BL I T ) **2 ) -XJ* ( X J3*TL* 
3XKK**2*PO (XKK.TE ) )/(4.* ( VE-BL I T ) **4 ) 

RETURN 

END 

FUNCTION PRES(VE.TE) PRES 

EQUATION OF STATE w0[rl - , T v , 

CO^ MON A2,A3,A4,A5,62.B3,34,B5.C2.C3.C4,C5,XKK,R,TREF,VREF.dLIT,XJ 

COMMON D4 . Z4 .ESS . TUFF .CLOD. GEE » A JA , KX . P2 I ( 30 ) . U2 I ( 30 ) « H2 I < 30 ) 

COMMON PT2 ( 30 , , RHO ( 30 ) . TE2 I ( 30 ) . VE2 I ( 30 ) . CV2 ( 30 ) . CP2 ( 30 ) . GAM2 ( 30 ) 
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APPENDIX C - Continued 


Program B - Continued 


c 


c 


COMMON ACCS ( 30 ) « XMACH2 ( 30 ) i ROPAT (30 ) * TRAT (30 ) *PRAT (30 ) * PTR AT (30 ) 
COMMON E I G64 *36 « A6* AL 
PO ( XL « AT > = 1 ./EXP (XL* AT ) 

POT =PO ( XKK * TE > 

PRES =R*TE/ (VE-BLI T)+( A2 +B2*TE+C2 *POT ) / ( VE-BL I T ) **2 + ( A3+B3*T£ 

1 +C3*POT )/ ( VE-BL I T ) **3 + ( A4+BIG84*TE)/( VE-BL IT ) **4+ ( A5+BB* TE +C5*P0 T ) 
?/ (VE-BL I T )**5 + ( A6+B6*TE ) *EXP ( AL*VE ) 


PRES 

PRES 

PRES 

PRES 


RETURN 

END 

FUNCT I ON CEP ( VE * TE ) 

CP 

COMMON A2< A3iA4 * A 5 « B2 <B3«64*B5*C2«C3*C4*C5f XKK * R « TREF * VREF ♦ BL I T « X J 
COMMON D4 i Z4 « ESS * TUFF « CLOD * GEE » AJA »KX»P2 I (30) * U2 I ( 30 > « B2 1(30) 
COMMON PT2 (30) * RHO ( 30 ) « TE2 1(30) *VE2I (30 ) * CV2 ( 30 ) *CP2 (30 ) * GAM2 ( 30 > 
COMMON ACC2 < 30 ) ♦ XMACH2 ( 30 ) * ROR A T ( 30 ) <TRAT (30 ) * PRAT ( 30 > « P TR AT ( 30 > 


COMMON 5 IGS4 •B6»A6#AL 

EQUIVALENCE (C3*G3> * <C5*XJ3 ) , (C2*D3> 

PO ( XL i AT ) = 1 ./EXP(XL*AT) 

OEN=R#TE/ (VE-BL IT )**2+2#* ( A2+B2 *TE+C2*P0 ( XKK « TE ) )/ (VE-bL I T ) **3+3. 

?* ( A 3 +B3*TE + C 3* PO ( XKK « T E ) )/(VE-BLlT )**4+4.*A4/( VE-BLI T ) **5+5. * 

3 ( AB + BB* TE+C5+P0 ( XKK * T E ) )/(VE-BLlT>**6 
R I TE ( 6 * 1 00 ) DEN * TE ♦ VF 

EXAL=EXP ( AL* VE ) CEP 

DEN = DEN-4 . *0 I GB4* TE/ ( VE-BL I T ) **5 + A6*AL*EX AL+B6*AL*TE*EXAL CEP 

WR I TE (6*100) EXAL tDEN 

XNU = R/ (VE-BLlT)+(B2-XKK*C2*P0(XKK«TE))/< VE-BL IT)* *2+ (B3-XKK*C3* 
2P0(XKK«TE) >/( VE-BL IT ) ** 3+ < B5-XKK*C5*P0 ( XKK f TE ) ) / ( VE-BL I T ) * *5 
V/R I TE (6 ♦ 1 00 ) XNU 

XNU =XNU+BIGB4/ (VE-3L I T )**4+B6*EXAL CEP 

WR ITE(6*100) XNU 


XNL M=TF*X J*XNU**2 
*'R I TE ( 6 • 100 ) XNUM 
C C *P =XNUM/DEN 
WRITE(6fl00) CEP 
1 00 FORMAT ( 3E20. 8 ) 

RETURN 

END 

SUBROUT I NE EQU3 ( TE * VE * A NS3 ) 

ENTHALPY 1 

COMMON A2«A3*A4*A5* B2 * B3 #B4 # B5 * C2 «C3»C4 »C5i XKK «R * TWEF i VREF <dL I T « X J 
COMMON D4 , Z4 ,ESS * TUFF tCLODi GEE • AJA t KX • P2 I < 30 ) « U2 I (30 ) • H2 I ( 3 0 ) 
COMMON PT2 (3 n ) • RHO ( 30 ) * TE2 I (30 ) ♦ VE2 I ( 30 ) ♦ CV2 < 30 ) ♦ CP2 (30 ) * GAMS ( 30 ) 
COMMON ACC2 ( 30 ) *XMACH2 ( 30 ) * ROR AT ( 30 ) • TRAT ( 30 ) * PRAT ( 3 0 > * PTR AT ( 30 ) 
COMMON 5 I GB4 ,B6 ♦ A 6 * AL 

EQUIVALENCE (C3 »G3 ) * ( C5 * X J3 ) ♦ ( C2 « D3 ) 
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APPENDIX C - Continued 


Program B - Continued 


PO (XLiAT ) = 1 ./EXP(XL#AT } 

ANS3=Z4* (TE-TREF ) + B4* ( Tt**2-TREF**2 ) /2 . + C4* ( TE**3-TREF **3 ) /3 . + D4 * 
3 (T£**4-TREF**4 ) /4 .+ (PO ( XKK , TE ) * (XKK*TE+ 1 • ) -PO ( XKK , TREF )* ( XKK #TRE F 
3+1 . ) )*(XJ*D3/(VE-BLIT)+XJ*G3/(2.*<VE-BLIT)**2)+XJ*XJ3/<4.*(VE- 
4BL IT ) **4 ) )+XJ* ( VE*R* ( TE-TREF )/ ( VE-3L I T ) +VE* (B2* (TE-TREF )+C2* (PO 
5 (XKK ,TE ) -PO ( XKK , TREF ) ) ) / ( VE-Bl I T ) **2 > + XJ* ( VE* (B3* (TE-TREF ) + C3* ( PO 
fc (XKK « TE ) -PO ( XKK , TREF ) ) ) / ( VE-3L I T >**3 )+XJ* (VE# (B5* (TE-TREF ) + C5* (PO 
7 (XKK , TE ) -PO (XKK , TREF ) ) ) / ( VE-BL I T )**5 )- <R*TREF* ( ALOG ( ( VE-BL IT )/ ( VR 
OEF-BL IT)) +BL IT/ ( VREF-BL I T ) -5L I T / (VE-BL IT ) ) ) *XJ- (2.# ( A2+B2*TREF+C2 
9*P0 (XKK , TREF ) >* ( 1 ./ ( VREF-3L I T >+BL I T/ ( 2 . * ( VREF-BL I T )**2 >-l . / (VE- 
ABL IT)-BLIT/(2.*( VF-BL I T )#*2 ) ) ) *XJ- ( 3. * ( A3+B3# TREF+C3*PO ( XKK , TREF 
R))*(l./(2.*( VREF-BL IT )**2 J+BLI T/ ( 3 . * ( VREF-BL I T ) **3 )- 1 ./ (2. * < VE- 
C3L I T )##2 )-3L IT/ (3.* ( VE-BL IT )**3 ) ) )*XJ- (4.*A4* ( 1 ./ (3.# ( VREF-dLI T ) 
D**3 )+BL I T/ (4 .* (VREF-BL I T >**4 >-l ./ ( 3.* (VE-BL IT )**3 )-BL IT/ (4.* ( VE- 
FBLIT )**4 ) ) )*XJ 

ANS3=AnS 3- ( ( 5. * ( A5+B5*TREF+C5*P0 ( XKK, TREF )))*(l./(4.« (VREF-dLI T ) 
2# #4 ) +3L I T/ (5.# (VREF-BL IT)**5)-l./(4.* (VE-BL I T)**4)-BLIT/(5.*(VE- 
3BL IT )**5 ) ) )*XJ+TREF*R*ALOG ( (VE-BL I T )/ < VREF-BL 1 T ) >*XJ-TREF* (32-XKK 


4*C2#P0( XKK .TREF ) )* ( 1 ./( VE-BL IT )-l ./ ( VREF-BLIT) )*XJ-TREF* (B3-XKK* 
5C3*P0( XKK, TREF ) )/2.*< 1 ./(VE-BL I T >**2-1 ./( VREF-BL I T > **2 ) *XJ-TREF * 

6 (B5-XKK*C5*P0 ( XKK , TREF ) ) /4 • * (1 ./ (VE-BL I T ) **4-1 ./(VREF-BLIT ) **4 ) 

7*XJ 

FXAL=EXP ( AL*VE ) EQU3 

FXAR=EXP ( AL*VRFF ) EQU3 

TEMP=X J*VREF* ( (BIGB4* (TE-TREF )/ (VREF-BL I T ) **4 ) +B6*EXAR* (TE-TREF ) )E0U3 

DUM=BIGB4*TE/ (VE-BL 1 T >**4 EQU3 

DUM1 =A6*EXAL+B6*TE*EXAL E0U3 

TEMP=TEMP+XJ*VE* ( 0.5*VE* ( -4. *DUM/ < VE-BL I T )+AL*OUMl )+ DUM+DUM1 ) EQU3 

DUM=BI GB4*TE/ ( VREF-BL I T ) **4 EOU3 

DUM1 = A6*EXAR + E36 + Tt*EXAR EQU3 

TEMP=TEMP-X J*VREF* ( 0.5*VREF# ( -4 .*DUM/ (VREF-BL I T >+AL*DUMl ) + EQU3 

1 CUM+DUM 1 ) EQU3 

TEMP = TEMP+ X J+ ( B IGB4+TE* ( 1 ./ ( 3.* ( VREF-BL IT)**3)-l./(3.*< VE-BL IT) EQU3 

1**3 ) ) +B6* TE* (EXAL-EXAR )/AL ) EQU3 

ANS3=AnS3+TEMP EQU3 

RETURN 
END 

FUNCTION PART(V,T) 

PARTIAL P W.R.T. RHO PART 


COMMON A2, A3 ,A4 ,A5, B2 ,63.34 ,B5 » C2 . C3 , C4 ,C5. XKK , R , TREF . VREF ,dL I T , XJ 
COMMON 04, Z4. ESS. TUFF. CLOD. GEE. AJA ,KX,P2 I ( 30 ) . U2 I ( 30 ) , H2 I < 30 ) 
COMMON PT2 (30 > , RHO ( 30 ) . TE2 1(30) , VE2 1(30) . CV2 ( 30 ) , CP2 (30 ) .GAM2( 30 ) 
COMMON ACC2 ( 30 ) , XMACH2 ( 30 ) , RORAT ( 30 ) , TRAT ( 30 ) .PRAT { 30 ) « P TR AT ( 30 ) 
COMMON B1GB4.B6.A6.AL 

EOU I VALENCE (C3.G3) , (C5.XJ3) , (C2.D3) 


APPENDIX C — Continued 
Program B — Continued 


PO (XL . AT ) = 1 ,/EXP ( XL*AT ) 

PART=R*T*V**2/ ( V-3L I T ) **2+4 . * A4*V**2/ ( V-BL I T ) **5+2. *V* *2* ( A2 + B2* 

2T + C2*P0 (XKK iTl )/ (V-BL I T ) **3 + 3 . * V**2* ( A3 + 83*T + C3*P0 ( XKK . T ) )/ ( V-BL 
3 I T >**4+5.*V**2* < A3+B5*T+C5*P0 (XKK ,T ) )/ (V-BL IT )**6 
EXAL=EXP (AL*V> PART 

PART=PART- ( V**2 ) * ( — 4 » *B I GB4* T / ( V-6L I T ) **5+A6*AL*EXAL+B6*T*AL*EXAL JPART 
RET JRN 
END 

FUNCTION FONDX(VE) 

C PARTIAL H W.R.T. T EONDX 

COMMON/PRTS/ PXTV * PVT V 

COMMON A2 . A3 * A4 .A5.S2 .B3.B4 «35.C2 *C3 . C4 « C5 . XKK . R , TREF . VREF .BL I T . XJ 
COMMON D4 , Z4 ,ESS . T'JFF , CLOD . GEE . A JA , KX . P2 I (30 ) *U2 I (30 ) «H2 I ( 30 ) 

COMMON PT2 (30 ) . RHO ( 30 ) . TE2 1(30) .VE2I (30 ) . CV2 ( 30 ) . CP 2 (30 ) .GAM2<30 ) 

COMMON ACC2 ( 30 ) . XMACH2 ( 30 ) * RORA T ( 30 ) . TRAT (30 ) .PRAT (30 ) .PTRAT (30 ) 

COMMON BIG64.B6.A6.Al 
CCMNON/DOPE/A ,B ,C. .PRAY 

EQUIVALENCE (C3.G3> . ( C5 . X J3 ) . ( C2 .03 ) 

PO ( XL. AT > = 1 ./EXP ( XL*AT ) 

TE=PRAy 

FONDX = Z4+B4*TE + C4*TE**2-XJ* ( ( D3*TE* (XKK** 2 ) *PO (XKK ,TE ) )/ ( VE-BL IT)) 


I -XU* ( ( G3*TE* ( XKK * * 2 )*PO ( XKK . TE ) )/ ( 2.0* ( VE-BL IT) **2 > ) 

F — XU* ( (XU3*TE*( XKK* *2 )*PC(XKk,TE) )/ (4 .0* ( VE-BL I T ) **4 ) ) 

3+VE*PXTV+D4*TE**3 

EX AR = E xP ( AL* VR FF ) FONDX 

FXAL =EXP ( AL* VE ) FONDX 

DUM = B I GB4/ ( VE -BL I T ) ** 4 FONDX 

CUM1=R6*FxAl FONDX 

TEMP=XU* VE * (DUM+DUM 1 +VE* ( -4. *DUM/ ( VE-8L I T )+AL*DUMl ) ) FONDX 

TEMP = TEMP-0. 5* XU* VREF* VREF* (-4. *B I GB4/ ( VREF-BLIT ) **5+b6*AL *EXAR ) FONDX 
FONDX = TEMP+FONDX + XU* ( 3 I GB4* ( I ./ (3.* (VREF-BLIT ) **3)-l ./ (3.* (VE- FONDX 
1 BL I T ) **3 ) )+ ( B6/AL )* (FXAL-EXAR ) ) 

RETURN 

END 

FUNCTION FONCX(VF) 

C PARTIAL OF P(V.T) W.R.T. T FONCX 

COMMON/PRTS/ PXTV.PVTV 


COMMON A2 . A3 . A4 . A5. 82 .63.84 ,B5 . C2 . C3 «C4 .C5. XKK ,R , TREF .VREF ,3L I T . XU 
COMMON D4 ,Z4 .ESS . TUFF .CLOU. GEE . AUA ♦ KX . P2 I ( 30 ) • U2 I ( 30 ) » H2 1 ( 30 ) 
COMMON PT2 ( 30 ) . RHO ( 30 ) . TE2 I ( 30 ) . VL2 I ( 30 ) . CV2 ( 30 ) »CP2 ( 30 > ♦ G A M 2 ( 30 ) 
COMMON ACC2 ( 30 ) . XMACH2 ( 30 ) . RORAT ( 30 ) . TRAT ( 30 ) .PRAT (30) . PTRAT ( 30 ) 
COMMON 3IGB4 .B6.A6, AL 

ECU I VALENCE (C3.C3) . (C5.XU3) , (C2.D3) 

COMMON/DOPE/ A , B , C , PRA Y 
PO < XL , AT ) = I ./EXP ( XL*AT ) 
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APPENDIX C - Continued 
Program B — Continued 


EL = 1 44. 

TE = PR A Y 

FONCX=R / <VE-BLIT)+(B2-C2*XKK*P0<XKK,TE) )/ ( VE-BL IT >**2 
1 + ( B3-C3*XKK*PO ( XKK iTE) ) / ( VE-BL I T ) **3 + < B5-C5*XKK#P0 < XKK • TE ) ) 

2/ ( VE-BL 1 T ) **5 

FONCX=FONCX+B I GB4/ < VE-BL I T > **4+B6*EXP ( AL*VE ) 

RETURN 
END 

FUNCTION FONDV(VE) 

PARTIAL H W • R • T • V 
COMMON/PRTS/ PXTV « PVTV _ 

COMMON A2»A3»A4 »A5«B2»B3«B4 iB5' C2 *C3*C4«C5« XKK 1 R 1 TREF i VREF .6LIT.XJ 
COMMON D4 , Z4 . ESS . TUFF , CLOD . GEE , A JA , <X , P2 I < 30 ) *U2 1(30 ) . H2 1 * * 
COMMON PT2 ( 30 ) , RHO ( 30 > . TE2 I < 30 ) . VE2 I < 30 > . CV2 ( 30 ) *CP2 ( JO > * GAM2 ( 30 ) 
COMMON ACC2 ( 30 I . XMACH2 < 30 M RORAT ( 30 > , TRAT ( 30 > < PRAT ( 30 ) .PTRAT (30 ) 
COMMON BI G94 *B6 • A6* AU 
COMMON/ DOPE/ A »B «C » PR A Y 

EQUIVALENCE (C3«G3)« (C5.XJ3) . (C2.D3) 

PO (XL.AT) = 1 ./EXP(XL*AT ) 

TE-PRAY 

fondv=ve*pvtv +te*pxtv +(C**2)*VE/< GEE*AJA) 

EXAL = EXP <AL*VE ) 

DUM=BIGB4*TE/ (VE-BL IT )**4 
DUM] =A6*EXAl 

?EMP=I!oIxJ*VE*(-4.*DUM/(VE-BLIT) + AL*(DUMl + DUM2)> + XJ*(2.0*(DUM + 

1 FONDV=FONDV+TEMP+0.5*XJ*VE*VE*(20.*DUM/(VE-BLIT)**2+(AL#*2>* 

1 (DUM1 +DUM2 ) ) 

RETURN 

END 

FUNCTION FONCV(VE) 

PARTIAL OF P(V»T> W*R*T* V 
COMMON/PRTS/ PXTV i PVTV 


FONCX 


FONDV 


FONDV 

FONDV 

FONDV 

FONDV 

FONDV 

FONDV 

FONDV 

FONDV 


FONCV 


COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
EQUI VAL 


P RTS/ rX I V ♦r V I V - l/Drc n, TT Vj 

A2«A3»A4 »A5« B2 » B3 ♦ B4 ♦ B5 « C2 »C3«C4«C5» XKK «R » TREF « VREF *BLIT«XJ 

D4,Z4.ESS,TUFF,CLOD,GEE,AJA,KX.P2I(30).U2I(30).H2I(30) 

PT2 < 30) ,RH0(30 ) » TE2 1(30) * VE2I ( 30 ) «CV2(30 ) .CP2 (30 ) «GAM2(30 ) 
ACC2 ( 30 ) » XMACH2 ( 30 ) » RORAT (30 ) «TRAT (30) « PR AT (30) «PTRAT ( 30 ) 
B!GB4fE6»A6»AL 

NCE { C3 * G3 ) « ( C5 * X J3 ) * (C2*D3> 


CO MM QN /DOPE /A ♦ B , C*PRAY 
PO ( XL t AT ) =1 */EXP ( XL* AT ) 

EL = 1 */l 44 * 

FONC V=I-R*TE )/ ( VE-BL I T ) **2-2 . * ( A2+B2*TE+C2*PO < XKK . TE > ) / < VE-BL I T ) 
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APPENDIX C - Continued 


Program B - Concluded 


1 **3-3 • 0* ( A3+B3*TE+C3*P0 (XKK t TE ) )/ < VE-BL I T ) **4-4.0*A4/ (VE-BL I T ) **5 
P-5 .0* ( A5+35*TE+C5*P0 ( XKK «TE ) )/ C VE-BL I T ) **6 + < EL/GEE ) *C**2 


POT=PO ( XKK t TE ) 

FONC V= ( -R*TE )/(VE-0LIT)**2-2.O* ( A2+B2*TE+C2 *POT ) / ( VE-BL I T ) **3 


FONCV 
FONC V 


1- 3.0*(A3+B3*TE+C3*POT>/ (VE-BL IT } **4 -4 • 0* < A4 +B I GB4 *TE )/( VE-BL I T ) **5F0NCV 

2- 5.0* ( A5+B5*TE+C5*P0T )/(VE-BLlT)**6+AL* ( A6 + B6*TE )*EXP (AL*VE ) FONCV 

3+ < EL /GEE ) *C * *2 FONCV 

RETURN 

END 


-00 

SD A T A A2 = -2. 162959tA3=4.04057E-3tA4=l .921 072E-4 # A5 = -4 • 48 I 049E-6 t 
XKK=.976E-2tB2=2. 1351 14E-3te3=l . 2828 1 8E-5 , B4 = 3 . 00892783E- 4 t 
05=9.06231 8E-9 * XJ=. 1 8505 3* 

C2=- 1 8.94 1 131 tC3=0.539776«C4=-l *30237441E-7*C5=-4. 836678E-5 t 
R = 0 . J 2 1 9336 .BL I T = C . 1 5E-2 # 04 = 1 .96802 894E- 1 1 ,PEE=0 • 1 026 t 
I END = 5 « 

VI (1 ) =3000 * 3500 t 4000 # 4500 t 5000 • 

TC ( I ) = 1 344 * 1610* 1 91 0 « 2250 *2600* 

VC ( 1 )=l4.o, 13. 2*12*1 .11. 7, M . 3 t 
PEI ( 1 ) =5*3 * 87E- l t 
A 1 = 1 5 « 

TEE = 260 t GEE = 32 #174* A JA = 778.3 t A= . 008 *B = 1 2000 t DELTX = 2 ♦ 

TREF = 82 0 t PREE= 1 t T Z = 4 09 . 5 t Z4 = 1 .904 580S4E-2 t 
TE I =5*532 t El = #00001 tE2=.0000 1tMAXr =30 ♦ NUMB =5 ♦ I NT = 0 t 
A6=5.83882E7 t B I GB4= -3 . 9 1 8263E-7 t AL=-6.61 l99E2t B6=-9 • 26392E4 t 
PEF= .387$ 
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n | ; Joi 7 PROGRAM NC 

LAR 


Program C 

NASA - LANG L EY RESEARCH CENTER 
COMPUTER PROGRAM ABSTRACT 


10/68 


01 JO TITLE OF PROGRAM f&l CHARACTERS MAXIMUM 

Total Conditions Behind Shock 


PARENT PROCRAI . 1 

"oThi 1 07 15 SITE |0? IB PROGRAM HO 

CATEGORY 


0? 76 

CATEGORY 

0? ?7 

language 
NO \ 

0? 3? 

LANGUAGE 
NO 2 

ED 

F«5R6H 

WHO 


K. A. Smith 


"oTmiTt Ob 3 X ORGN CODE P5 39 PROJECT NO 

^ cente 

L.AR 11.170 RGL-206 


OS 46 STATUS 

| \ A UNDER DEVELOPMENT 

1 \ B OPERATIONAL 

m c COMPLETED 


(ft 49 

QA THIS PROGRAM 
IS NOT FOR 
SHARING 


DATES 

05 so initiated Ids 54 completed” 


0k/6B 10/68 


(ft to REVISION CODE 

| 1 A REVISION 

I 1 B CANCELLATION 


[ TIME AND COST TOR DEVELOPMENT 

0S 59 MANMONTHS 

05 64 MACHINE 
HOURS 

1 111 

05 69 COMPUTER TYPE 

6000 

05 TOTAL COST 

* DOLLAR 

1111 J 1— 

59 SO 61 67 63 

64 65 66 67 66 

74 75 76 77 78 79 60 

ELITE MARGIN 1 

PICA MARGIN 


CARD NUMBER 
06 


O ABSTRACT — I— 

yj-i-.h the static conditions behind a given shoc k In CFi^ 
kno>m.^the~FotaT~c onditi ons b ehin d this shock are obtained ty — 

' 'bringing the_floy t o rest , i sent ropi c^lly * __^e governing 
mRnt in this cOE Tpression process is that the 
— • - — - 1 — ^ ins constant. 





APPENDIX C - Continued 
Program C - Continued 


PROGRAM PCH ( I NPuT. OUTPUT .TAPES® INPUT, TAPE6=0UTPUT,TAPE8 ) 

COMMON A2, A3.A4 .A5.B2.B3 »B4 *05 » C2 » C3 «C4 • C5.XKK.R ♦ TREF « VREF «BL IT, XJ 
COMMON 04 , Z4 « ESS * TUFF , CLOD » GEE . A JA ,KX ,P2 I (30 > .U2 t ( 30 ) »H2 1 ( 30 > 
COMMON S2.H0 

COMMON PT2 (30) « RHO ( 30 ) » TE2 1(30) .VE2I (30 ) ,CV2(30> ,CP2(30 ) «GAM2(30 ) 

COMMON ACC2 ( 30 ) » XMACH2 ( 30 ) .RORAT (30 ) «TRAT (30 ) .PRAT (30 ) .PTRAT (30 ) 

COMMON BIGB4.B6.A6.AL 

COMMON/CLUX/TC .VC 

COMMON/DESP/VOODOO . ZULU 

EXTERNAL FOP 

DIMENSION TC ( 30 ) ♦ VC ( 30 ) 

DIMENSION Pi ORAT ( 30 ) .T10RAT(30) .REN(30> 

COMMON /RATS/PEE. TEE. PT20RT (30 ) .P1T2RT(30 ) .T20RT (30 ) 

DIMENSION NO (30) 

EQUIVALENCE (C3.G3). (C5.XJ3), (C2.D3) 

PO (XL « AT ) * 1 ./EXP ( XL* AT ) 

DIMENSION TEI (30 ) . VEI (30 ) .PEI ( 30 ) . VI (30 ) .GAM (30 ) .SPEED (30 ) 

D I MENS I ONXMACH ( 30 ) . CP ( 30 ) . CV ( 30 ) . H ( 30 ) ♦ V I S ( 30 ) 

DIMENSION P2RT(30).P2T(30).P2T1 ( 30 ) « T2TT ( 30 ) . V2TV ( 30 ) .S2(30) 
DIMENSION AMAT (2.2 ) .BMAT (2 ) . IPIV0T(2) , INDEX (2 .2) 

NAMELIST /DATA/A2. A3, A4 .A5. A6.B2.B3.BIGB4.B5.B6.B4.AL. 
1C2.C3.C4,C5,XKK.XJ.R.BL I T, D4, PEE .TEE. GEE. A JA, NUMB. INT.Z4.A ,B, 
2DELTX. IEND.TEI .El .E2.MAXI ♦ TREF « PREF. TZ. TC , VC . THETA , 

3VEE . HO , PRA T ♦ XM ACH , NO , VE I . VREF , S2 .U2 I 
99 READ (5, DATA) 

WRITE (6. DATA ) 

WRITE(6.aB)A2,A3. A4.A5.B2.B3.B4.B5 
88 FORMAT ( 1 0H0A2 5 B2 F//0E13.5) 


.C3.C4.C5, XKK , XJ ,R ,BL I T 

89 FORMAT ( 20H0C2 5 K J R B/8E13.5) 

WR I TE ( 6 . 90 ) D4 , PFE , TEF . GEE . A JA . NUMB , I NT . Z4 

90 FORMAT (23H0D4 PO TO G J N I /5E 1 3.5 , 2 1 1 3 .El 4 • 5 ) 
CALL E0U3 (TFF , VFE , ANS3 ) 

DO 1 I»1,IEND 

S2 ( I )=ENT(VEI ( I ) , TEI ( I ) ) 

PRINT 9.S2 
9 F ORMAT ( 4E 16.8) 

H0 = ANS3+200.+ (U2I ( I )**2 >/ ( 2 . 0»GEE*A JA ) 

PRINT 9. HO 

DO 22 J»t .50 

PRINT 9. VC ( I ) . TC ( I ) 
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Program C - Continued 


C H(T.V) 

HTVaFOND (VC ( 1 ) » TC ( I ) ) 

C S { T ♦ V ) 

STV*ENT(VC U)tTC(I) >-S2 ( I ) 
C HSUBT(T.V) 

HST =FONT (VC ( I ) . TC ( I ) ) 


C 

c 

c 

c 


c 


HSUBV ( T ♦ V ) 

HSV*FONV(VC < I ) .TC ( I ) 1 
S SUBT(T.V) 

SST»ENTT(VC( ! ) .TC( I ) ) 

5 SUBV(T.V) 

SSV*ENTV (VC ( n »TC ( I ) ) 

SET UP MATRICES FOR NEWTON-RAPHSON 

AMAT { 1 . 1 >=HST 

AMAT (2.1 ) *SST 

AMAT ( 1 » 2 ) *HSV 

AMAT (2.2) *SSV 

BMAT ( 1 ) = — HTV 

BMAT ( 2 ) *—STV 

CALL MAT INV ( AMAT. 2 . OMAT . 1 . DETERM , I P I VOT , I NDEX .2.1 SCALE ) 
DELI «ABS (BMAT ( 1 )/TC < I ) > 

DEL2=AB5(BMAT(2)/VC( I I ) 

PRINT 9 . DEL 1 . DEl2 . BMAT 

IF (DELI «LT. El .AND. 0EL2.LT. E2 ) GO TO 33 
TC ( I ) = TC ( I )+BMAT ( 1 > 

VC ( I )=VC ( I ) +BMAT ( 2 ) 

22 CONTINUE 
PRINT 101 

101 FORMAT (24H0MAX. ITERATIONS EXCEEDED) 

33 PRINT 100 .NO ( I ) .PRAT ( I ) *VEI ( I ) *TE1 ( I ) .XMACH ( I ) . 

1 VC ( I ) ♦ TC ( 1 ) 

100 FORMAT (IS. 6E 18.8) 


P(V.T) 

P2T € I > bFONC ( VC ( I) . TC ( I ) ) 
P2RT ( 1 1 *P2T ( I ) /PEE 
P2T1 ( 1 ) *P2RT ( I ) /PRAT ( I ) 
T2TT ( I ) *TC ( I l/TFE 
V2TV ( I )»VC ( I )/VFE 


1 CONTINUE 

PRINT 102. HO. SO 
102 FORMAT (4H1H0=E15 


8.5X.3HS0=E15.8//1H0. 10X.3HT2T. 17X.3HV2T ♦ 17X, 


13HP2T.18X.2HS2) 

PR I NT 1 03 . (TC ( I ) . VC < I ) . P2T ( I ) . S2 (I ) . I = 1 . I END ) 
103 FORMAT (4E20. 8) 

PRINT 104 
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Program C - Continued 


c 


c 


1 04 FORMAT ( 1 HI * 1 0X.6HP2T/P0 . 1 AX « 6HP2T/P 1 • 1 4X. 6HT2T/T0 . t 4X.6HV2T/V0 ) 
PRINT 1 03, (P2RT ( I > »P2T1 ( l ) , T2TT ( I ) «V2TV( I ) . 1=1 • I END ) 

WRITE(8) I END 

WRITE (8) ( XMACH { I ) ,P2RT( I ) .P2TI ( I > , 1 = 1 , I END) 

END FILE 8 
GO TO 99 
END 

FUNCTION ENT ( VEE » TEF ) 

ENTROPY 

COMMON A2 , A3 « A4 , A5, 82 »B3 ,B4 ,B5 . C2 • C3 « C4 » C5 . XKK , R , TREF « VREF .BL I T , XJ 
COMMON D4 .Z4.ESS. TUFF ♦ CLOD « GEE * A JA ,Kx , P2 I (30) ,U2I (30) ,H2I (30) 
COMMON S? « HO 

COMMON PT2 (30) ,PH0(30 ) * TE2 I (30) ,VE2l (30 ) . CV2 ( 30 ) . CP2 ( 30 ) ,GAM2(30 ) 
COMMON ACC2 ( 30 ) ,XMACH2 ( 30 ) .RORAT (30) «TRAT(30 ) .PRAT (30 ) .PTRAT(30 ) 
COMMON BIGB4.B6.A6.AL 
DIMENSION S2 ( 30 ) 

EQUIVALENCE (C3.G3) . (C8.XJ3) , (C2.D3) 

PO ( XL « AT ) = 1 «/EXP ( XL* AT ) 

ESS 1 =Z4*AL0G (TEE/TREF >+B4* ( TEE— TREF ) +C4* ( TEE* ♦ 2- TREF* *2 )/2. 
ESS2=ESS1 +D4* (TEE**3-TREF**3 ) /3 . +X J*03*XKK/ ( VREF-BL I T )* (PO ( XKK, TEE 
2 > -PO ( XKK , TREF ) ) 

ESS3=ESS2+X J#G3*XKK/ ( 2 . * ( VREF-BL I T ) ** 2 ) * ( PO ( XKK , TEE ) -PO ( XKK , TREF ) ) 
ESS4 = ESS3+X J*X J3*XKK/ ( 4 .* ( VREF-BL I T ) **4 > * (PO ( XKK .TEE ) -PO ( XKK « TREF ) 
2 ) 

ESS5=ESS4+ (R*ALOG ( (VEE-BL I T )/( VREF-BL !T ) ) )*XJ 

ESS6=ESSS- (B2-XKK*C2*P0 (XKK, TEE) )*( 1 ./ (VEE-BL I T )-l •/( VREF-BL I T ) ) *X 
2J 

ESS7*ESS6- CB3-XKK*C3*P0 (XKK, TEE ) )/2.» ( 1 ./ (VEE-BL I T ) **2- 1 ./ (VREF-BL 
2 I T ) **2 )*XJ 

ENT =ESS7- ( B5-XKK*CS*P0 ( XKK . TEE ) ) /4 . * ( 1 . / ( VEE-BL I T ) **4- 1 . / ( VREF-BL I 
2T)**4)*XJ 


EXAL*EXP(AL*VEF) 

EXAR=EXP ( AL4VREF ) 

ENT =ENT+XJ* (B I GB4* ( 1 ./ ( 3.* ( VREF-BL I T ) **3 )- 1 •/ ( 3 .* ( VEE-BL I T ) 4*3 > ) 

1 +B6* (EXAL-EXAR ) /AL ) 

ENT =ENT+0 .848 

RETURN 

END 

FUNCTION FONO(VF.TE) 

ENTHALPY FOR NFWTON— RAPHSON 

COMMON A2, A3. A4.A5.B2. B3.B4.B5.C2. C3.C4.C5. XKK, R, TREF. VREF, BLIT, XJ 
COMMON D4.Z4. ESS. TUFF. CLOD, GEE. AJA.KX.P2I (30 ) « U2 I ( 30 ) . H2 1 { 30 ) 
COMMON S2.H0 

COMMON PT2 ( 30 ) .RHO ( 30 ) . TE2 1(30) » VE2 1(30) .CV2 ( 30 ) <CP2 ( 30 ) ,GAM2 (30 ) 
COMMON ACC2 (30 ) .XMACH2 (30 ) .RORAT (30) , TRAT ( 30 ) . PRAT ( 30 ) . PTRAT ( 30 ) 


ENT 

ENT 

ENT 

ENT 


FOND 
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Program C - Continued 


COMMON BIGB4«B6«A6»AL 
DIMENSION S2 00) 

EQUIVALENCE ( C3 «G3 ) » CC5.XJ3) , (C2.D3) 

PO (XL » AT ) c 1 »/FXP C XL* AT ) 

H6 = Z4* (TE-TREF ) + B4* ( TE**2-TREF**2 ) /2 . +C4* ( TE**3-TREF**3 ) /3 . + D4* 
j> ( TE**4-TREF**4 ) /4 • + ( PO ( XKK « TF ) * (XKK*TE+ 1 • )— PO (XKK. TREF )* (XKK* TREF 
3+1 . ))* (XJ*D3/< VREF-BL IT >+XJ*G3/ (2 •* ( VREF-BL I T ) **2 >+X J*X J3/ (4 •+ 

4 (VREF-BLIT )**4 ) )+XJ* (VREF*R* (TE-TREF )/ (VREF-BL I T )+vREF* (B2* (TE- 
5TREF )+C2* ( PO ( XKK « TE ) — PO (XKK. TREF > ) )/( VREF-BL I T ) **2 ) +XJ* < VREF* ( B3* 

6 ( TE-TREF )+C3* ( PO ( XKK . TE ) -PO ( XKK * TREF ) ) ) / ( VREF -BL IT )**3 ) + XJ* (VREF* 

7 ( B5* ( TE-TREF ) + C5* ( PO ( XKK * TE ) -PO (XKK . T REE ) ) ) / ( VREF-BL I T ) **5 )- (R*TE* 
B ( ALOG( ( VE-BL IT ) /< VREF-BL I T ) >+BL I T/ (VREF^BL I T ) -BL I T/ ( VE-BL I T ) ) ) *X J 
H7=H6— ( 2» * ( A2+B2*TE+C2*PO ( XKK * TE ) )*( 1 * / ( VREF-BL I T ) +BL I T/ (2 . * ( VREF- 


2BL IT )**2 >- 1 •/( VE-BL I T )-BLl T / (2»* (VE-BL I T )**2 ) ) >*XJ 
H8*H7- ( 3.* ( A3+B3*TE+C3*P0 (XKK • TE > ) * ( 1 •/ ( 2.* (VREF-BL I T ) **2 > +BL I T/ ( 3 
2»* (VREF-BL IT) **3 ) — 1 •/ (2 •* ( VE-BL I T ) **2 ) -BL IT/ ( 3«* ( VE-BL I T )**3 ) > > *X J 
H9«H8-(4.*A4*( 1 . / ( 3** (VREF-BL IT)** 3 >+BL!T/(4.» (VREF-BL I T >**4 )- 1 ./( 
23 •* (VE-BL I T )**3 )-BL I T/(4.# (VE-BL IT )**4 > > )*XJ 
HI 0-H9- ( (5.* (A5+B5*TE+C5*P0 (XKK.TE >>)■*(! •/ (4 • * (VREF-BL 1 T )**4 )+BL IT 
2/(5**( VREF-BL I T )**5 )-l • / (4 • * ( VE-BL I T )**4 )-BL I T/ (5»* ( VE-BL I T >**5 ) ) ) 
3*XJ 

HI 1 *H1 0+TE*R*ALOG ( (VE-BL I T ) /(VREF-BLIT ) )*XJ 

H12*H1 1-TE* (B2-XKK*C2*P0(XKK,TE> )* ( 1 . / ( VE-BL I T )- I • / ( VREF-BL I T ) )*XJ 
HI 3*H1 2- TE* (B3-XKK*C3*PO (XKK . TE ) )/2.* ( 1 •/ ( VE— BL IT) **2-1 •/ (VREF-BL I 
2T)**2 )*XJ 

HI 4*H1 3-TE* (85-XKK*CS*P0 (XKK.TE ) )/4 .* ( 1 • / (VE-BL I T )**4- 1 ./(VREF-BL I 
2T)**4 )*XJ 
FOND*H 1 4— HO + 200 » 0 
EX AL“EXP ( AL*VE ) 

EXAR«EXP ( AL*VREF ) 

TEMP*X J*VREF* ( ( B I GB4* (TE-TREF ) / (VREF-BL I T )**4 ) 

DUM-B I GB4*TE/ ( VE-BL I T ) **4 
0UM1 »A6*EXAL+B6*TE*EXAL 

TEMP»TEMP+XJ*VE*(0.5*VE*(-4.*DUM/(VE-BLIT)+AL*DUM1 )+ DUM+DUM1 ) 

DUM»B I GB4*TE/ ( VREF-Ol I T ) **4 
DUM1 «A6*EXAR+B6*TE*EXAR 

TEMP*TEMP-X J*VREF* ( 0 . 5* VREF* ( -4 . *DUM/ ( VREF-BL I T ) +AL*DUM 1 ) + 

1 DUM+DUMl ) 

TEMP*TEMP+XJ* ( B I GB4*TE* (I •/ ( 3.* (VREF-BL I T )**3 )-l •/ ( 3.* (VE-BL I T ) 
1**3) )+B6*TE* (ExAL-EXAR ) /AL ) 

FOND=FOND+TEMP 
RETURN 
END 

FUNCTION ENTV ( V * T ) 

PARTIAL S W.R.T* V 


E0U3 

EQU3 

+B6*EXAR* (TE-TREF ) >EQU3 

E0U3 
EQU3 
EQU3 
E0U3 
EQU3 
EQU3 
E0U3 
EQU3 
EQU3 
FOND 
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Program C — Continued 


c 


♦VREF.BLIT.XJ 
, H2 ! ( 30 ) 

30) ,GAM2(30 ) 

) .PTRAT (30 ) 

DIMENSION S2 (30 ) 

PO (XL* AT ) = 1 » 0/FXP ( XL# AT ) 

VM = V-BL I T 

SUM* (R/ ( V— BL IT) )+(B2-XKK*C2*P0(XKK,T) )/(VM**2) 


COMMON A2 » A3.A4 .A5.B2 .B3.B4.B5 .C2.C3 .C4.C5.XKK.R.TREF 

COMMON D4 , Z4 *ESS * TUFF • CLOD * GEE * A JA .KX.P2 I (30 ) *U2 1 ( 30 ) 
COMMON S2.H0 

COMMON PT2 ( 30 ) « RHO ( 30 ) ♦ TE2 I ( 30 ) • VE2 I (30 ) ♦ CV2 ( 30 ) «CP2 ( 
COMMON ACC2 (30 ) .XMACH2 (30 ) .RORAT (30 ) *TRAT (30 ) .PRAT ( 30 
COMMON BIGB4.B6.A6.AL 


SUM*SUM+ (B3-XKK*C3*P0 (XKK. T ) )/ ( VM**3 ) 

ENTV = SUM+ ( B5-XKK*C5*P0 ( XKK « T ) ) / ( VM**5 ) 
ENTV=ENTV+XJ*(BIGB4/(V-BLIT>*#4 +B6*EXP ( AL*V ) ) 


RETURN 

END 

FUNCTION ENTT(V.T) 

COMMON A2.A3.A4.A5. B2 . B3 *B4 . B5 . C2 . C3 . C4 » C5 »XKK .R » TREF ♦ VREF .BL I T » X J 
COMMON D4.Z4. ESS. TUFF, CLOD. GEE. AJA.KX.P2I (30) .U2I (30) .H2I (30) 
COMMON S2.H0 

COMMON PT2(30) .RHO (30) . TE2 I (30) .VE2I ( 30 > . CV2 ( 30 ) . CP2 ( 30 ) .GAM2 ( 30 ) 
COMMON ACC2 ( 30 ) « XMACH2 ( 30 ) . RORAT ( 30 ) , TRAT (30 ) .PRAT ( 30 ) . PTRAT ( 30 ) 
COMMON BIGB4 ,B6, *6. AL 
DIMENSION S2 ( 30 ) 

EQUIVALENCE (C3.G3) , (C5.XJ3) , (C2.D3) 

PO(XL.AT) *1 «0/FXP(XL*AT) 

XK2*XKK*XKK 

VERB=VREF-BL IT 
V I B*V— BL I T 
POT*PO (XKK.T ) 

DUM* (Z4/T )+B4+C4*T+D4* (T**2 )-XJ4XK2*P0T* 

1 (D3/VERB+G3/(2.*(VFRB**2> ) + X J3/ ( 4 . 0* ( VERB**4 ) ) ) 

DUM=DUM-C2*XK2*P0T* ( ( 1 ./VIB)-( 1 ./VERB) ) 

ENTTsDUM-O.B*C3*XK2*POT*( ( 1 ./(V!B*»2 ) )-( 1 «/(VERB**2 ) ) > 

1 -0 .2S*XK2*P0T*CS* ( ( 1 ./ ( VIB**4 ) )- ( 1 ./(VERB*»4 ) ) ) 


RETURN 


END 

FUNCTION FONT ( VE » TE ) 

PARTIAL H W.R.T. T 
COMMON A2, A3.A4 .A5.B2.B3.B4 
COMMON D4 ,Z4 .ESS. TUFF .CLOD, 
COMMON S2.H0 

COMMON PT2 ( 30 ) »PHO ( 30 ) » TE2 I 
COMMON ACC2 ( 30 ) .XMACHP ( 30 ) , 
COMMON BIGB4.B6.A6.AL 
DIMENSION S2 ( 30 ) 


FONOX 

.B 5 .C 2 . C 3 ,C 4 .C 5 .XKK,R,TREF»VREF .BL I T ,XJ 
GEE , A JA .KX.P 2 I ( 30 ) ,U 2 I ( 30 ) ,H 2 I ( 30 ) 

( 30 ) . VE2 I ( 30 ) » CV2 ( 30 ) ,CP2 ( 30 ) .GAM 2 ( 30 ) 

RORAT (30 ) .TRAT (30 ) .PRAT ( 30 ) .PTRAT ( 30 ) 
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Program C - Continued 


EQUIVALENCE (C3.G3 ) ♦ (C5.XJ3 > , (C2«D3 > 

PO CXL.AT)*1 ./EXP(XL*AT ) 

pot«po(xkk.te> 

VER=VREF-BLIT 

VIB«VE-BL!T 

WVB=1 • /V I B 

WVR*1 ./VER 

B2C2«B2-C2*XKK:*POT 

B3C3«B3-C3*XKK*P0T 

B5C5»B5-C5*Xt<K*P0T 

VR2»WVR*WVR 

VR3»VR2*WVR 

VR5«VR2*VR3 

VB2=WVB*WVB 

VB3»WVB*VB2 

VB5-VB2*VB3 

XK2*XKK*XKK 

FONT «Z4+B4*TE+C4*TE**2-XJ*( (D3*TE* (XKK**2 )*P0 (XKK *TE ) >/(VE-BLIT) ) 

1 - X J* ( ( G3*TE* (XKK**2 > *P0 (XKK , TE ) ) / ( 2 . 0* ( VE-BL I T )**2 ) ) 

2- XJ*( < XJ3*TE* < XKK**2 ) *PO < XKK ♦ tE ) > / ( 4 • 0* (VE-BL IT >**4 ) ) 

3+X J *VREF* <R*WVP+B2C2*VR2+B3C3*VR3+B5C!5*VR5 ) -X J * (R* ( ALOG ( V IB*WVR ) 
4+BL I T* ( WVR-WVB ) ) +2 . 0*B2C2* ( WVR+0 • 5*BL IT* VR2-WVB-0 • S*BL 1 T*VB2 ) 

O+3*0*B3C3* (0.5*VR2+BLlT*VR3/3.0-0.5*VB2-BLlT*VB3/3.0) 

E+5*0*B5CS* ( 0 .25*VR2*VR2+BL I T*0 . 2*VRS-0 . 25*VB2*VB2-0. 2*BL I T*VB5 ) ) 
5+X J * ( R*ALOG ( V IB*WVR ) -B2C2* ( WVB-WVR ) -B3C3* ( VB2-VR2 J-B5C5* ( VB2*VB 2 


6-vR2*VR2 ) ) — XJ*TE*XK2*P0T* ( C2* ( WVB-WVR ) +0 . 5*C3* ( VB2-VR2 ) 

7+0»2S*C5* (VB2*VB2-VR2*VR2 ) >+D4* (TE**3 ) 

EXAR»EXP(AL*VREF) FONDX 

EXAL=EXP(AL*VE ) FONDX 

DUM*B!GB4/ (VE -BLIT)**4 FONDX 

DUM1»B6*EXAL FONDX 

TEMP=XJ* VE * (DUM+DUM1 +VE* ( -4»*DUM/ (VE-BL IT )+AL*DUM 1 ) ) FONDX 

TEMP*TEMP-0 • 5#X J* VREF*VREF* ( -4 « *B I GB4/ ( VREF-BL I T ) **5+B6*AL*EXAR ) FONDX 
FONT *TEMP+FONT +XJ* ( B I GB4* ( 1 • / ( 3»* (VREF-BL I T ) **3 )-l • / ( 3»* (VE- FONDX 

1 BL IT) **3 ) )+(B6/AL)*(EXAL-EXAR) ) 


RETURN 

END 

FUNCTION FONV(V.T) 

PARTIAL H W.R.T* V 1 

PO CXL*AT)*1 ./EXP(XL*AT > 

COMMON A2 * A3« A4 *AS«B2«B3*B4 *B5 * C2 • C3 «C4 « C5 «XKK»R «TREF »VREF *BL I T«XJ 
COMMON D4 *Z4 » ESS « TUFF « CLOD • GEE * AJA »KX« P2 I (30 ) »U2 I ( 30 ) «H2 I ( 30 ) 
COMMON S2 • HO 

COMMON PT2 (30 ) #RHO( 30 ) #TE2I ( 30 ) »VE2 1(30) «CV2 ( 30 ) .CP2 (30 ) »GAM2(30 ) 
COMMON ACC2 ( 30 ) «XMACH2 ( 30 ) ♦ RORAT ( 30 ) »TRAT ( 30 ) «PRAT (30 > «PTRAT (30 ) 
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Program C - Continued 


COMMON BIGB4 »B6«A6« AL 
DIMENSION S2 (30 ) 

VIB*! •/IV-BUIT) 

V2*VIB»VIB 

V3*V I B*V2 

V4*V2*V2 

V5»V2*V3 

V6*V3*V3 

POT*PO(XKK.T) 

DUM»R*T* (VIB+BLIT*V2)+2.0* ( A2+B2*T+C2*POT)*(V2+BLIT*V3 ) 
DUM*DUM+3 • 0* (A3+B3*T+C3*POT )* (V3+BLIT*V4 )+4.0*A4*(V4+BL I T*V5) 
DUM*-XJ * (DUM+5.0*(A«!+B5*T+C<3*P0T)* (VB+BLIT#V6 ) ) 

DUM 1 *XJ *T* (R*V IB+ (B2-XKK*C2*POT >*V2+ (B3-XKK*C3*P0T > *V3 
1 + ( B5-XKK*C5*POT ) *V5 ) 

FONV*DUM+OUM1 
EXAL*EXP(AL*V ) 

DUM=B I GB4*T /(V -BLIT)**4 


DUM1*A6*EXAL FONDV 

DUM2*B6*T *EXAL 

TEMP* 2 # 0*X J* V *(-4»*DUM/<V -BL I T ) + AL* (DUM1 +DUM2 ) )+XJ* (2 .0* (DUM+ FONOV 

1 DUM2 ) +DUM 1 ) FONDV 

FONV *FONV +TEMP+0.5*XJ*V *V *<20.*DUM/(V -BL I T ) **2 + ( AL**2 ) * FONDV 

1 (DUM1 +DUM2 ) ) FONDV 

RETURN 
END 

FUNCTION FONC(VF*TE) 

C EQ. OF STATE FOR NEWTON-RAPHSON FONC 


COMMON A2 « A3 « A4 « A5 . B2 • B 3 . B4 . B5 . C2 « C3 . C4 • C5. XKK .R . TREF ♦ VREF «BL I T « X J 
COMMON D4 «Z4«ESS* TUFF • CLOD « GEE * A JA »KX < P2 I (30) .U2I (30) *H2I (30) 

COMMON S2 * HO 

COMMON PT2 (30) .RHOOO ) ♦ TE2 I (30) .VE2I (30 > «CV2( 30) »CP2 (30) .GAM2<30 ) 
COMMON ACC2 (30 ) *XMACH2 (30 ) .RORAT (30 ) « TRAT (30 ) »PRAT ( 30 ) .PTRAT(30 ) 

COMMON BIGB4«B6«A6*AL 
DIMENSION S2 ( 30 ) 

EQUIVALENCE (C3.G3) « (C5.XJ3) , (C2«03) 

PO (XL.AT)*1 ./EXP(XL*AT ) 

EL * 1 • / 1 44 • 

POT *PO ( XKK « TE ) 

FONC*R*TE/ (VE-BL I T ) + ( A2+B2*TE+C2*POT )/( VE-BL I T )**2+ ( A3+B3»TE+ FONC 

1 C3*POT ) / ( VE-BL I T ) **3+ ( A4+B I GB4*TE )/ (VE-BL I T )**4+ ( A5+B5*TE+C5*P0T ) FONC 
2/ (VE-BL I T ) **5+ ( A6+B6*TE ) *EXP ( AL*VE ) 

RETURN 

END 

SUBROUTINE EQU3 ( TE* VE • ANS3 ) 

C ENTHALPY EQU3 
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Program C - Continued 


COMMON 'A2.A3.A4,A5.B2.B3.B4,BS.C2.C3.e4.C5«XKK,R,TREF.VREF.BLIT.XJ 
COMMON D4 «Z4 . ESS • TUFF . CLOD « GEE * A JA . KX • P2 I (30 ) »U2 I (30 ) »H2 I ( 30 ) 
COMMON S2 « HO 

COMMON PT2 (30 ) ,RH0 (30 ) . TE2 1(30) » VE2 1(30) «CV2 ( 30 ) » CP2 ( 30 ) »GAM2( 30 ) 
COMMON ACC2 (30 ) .XMACH2 ( 30 ) .RORAT ( 30 ) . TRAT ( 30 ) « PRAT (30 ) * PTRAT ( 30 ) 
COMMON B I GB4 « B6 « A 6 « At_ 

EQUIVALENCE (C3.G3) • (C5.XJ3) « (C2.D3) 

DIMENSION S2 (30 ) 

PO (XL ♦ AT ) = 1 »/ExP ( XL*AT ) 

ANS3*Z4* (TE-TREF )+B4* ( TE**2-TREF**2 ) /2 . +C4* ( TE**3-TREF**3 ) /3 .+D4* 
2(TE**4- TREF**4 ) /4 .+ (PO (XKK.TE )* (XKK*TE+ ! • ) -PO ( XKK , TREF )* (XKK*TREF 
3+1 • ) )*(XJ*D3/(VE-BLlT)+XJ*G3/(2.*(VE-BLIT)**2)+XJ*XJ3/(4.*(VE- 
4BLIT)**4 ) ) + XJ* (VE*R* (TE-TREF)/ (VE-BL IT ) +VE* ( B2* ( TE-TREF ) +C2* (PO 
5 (XKK « TE ) -PO { XKK , TREF ) ) ) / ( VE-BL I T ) **2 ) +X J* ( VE* ( B3* (TE- TREF ) +C3* ( PO 
6 (XKK , TE ) -PO ( XKK , TREF ) ) ) / (VE-BL I T )**3 )+XJ* ( VE* ( B5* (TE-TREF ) +C5* ( PO 
7 ( XKK , TE ) -PO ( XKK , TREF ) ) ) / ( VE-BL I T ) **5 ) - ( R*TREF* ( ALOG ( ( VE-BL IT )/ (VP 
8EF-BL I T ) )+8L I T/ ( VREF-BL I T ) -BL I T/ ( VE-BL IT))) *X J- (2 .* ( A2+B2*TREF+C2 
9*P0 (XKK « TREF ) ) * ( 1 . / ( VPEF-BL I T ) +BL I T/ ( 2 • * ( VREF-BL I T ) **2 )-l./(VE- 
ABLIT)-BLIT/(2 .*(VE-Bl!T )** 2 > ) ) *X J- ( 3.* ( A3+B3*TREF+C3*PO ( XKK . TREF 
B) )*( 1 ./ (2.* (VREF-BL IT ) **2 >+BLIT/(3.*(VREF-BLIT)**3)-l ./(2.*(VE- 
CBL IT)**2 ) -BLIT/ (3.* (VE-BL IT )**3) ) )*XJ-(4.*A4* ( 1 ./ ( 3. * ( VREF-BL I T ) 
0**3 ) +BL I T/ ( 4 .* ( VREF-BL IT)**4 j-1./(3.*( VE-BL I T ) **3 ) -BL I T/ (4 . * ( VE- 
EBL IT) *#4 ) ) )*XJ 

ANS3=ANS3-( ( 5. * ( A5+BB*TREF+CB*P0 (XKK ♦ TREF ) ) )*( 1 ./ (4 •* ( VREF-BLI T ) 
2**4 ) +BL I T/ ( 5 .* ( VREF-BL I T ) **5 ) - 1 . / ( 4 . * ( VE-BL IT) **4 > -BL I T/ ( *5 . * ( VE- 
3BL IT) **5 ) ) ) *XJ+TREF' t R* ALOG ( (VE-BL I T )/ (VREF-BL IT ) )*XJ-TREF* (B2-XKK 


4*C2*PO (XKK « TREF ) )* ( 1 ./(VE-BL IT )-l ./( VREF-BL I T ) )*XJ-TREF* (B3-XKK* 

5C 3*P0 ( XKK ♦ TREF ) )/2.*( 1 ./( VE-BL I T ) **2- 1 ./( VREF-BL I T ) **2 ) *XJ-TREF* 
6(B5-XKK*C5*PO (XKK.TREF ) >/4.* ( 1 ./ (VE-BL IT ) **4-1 ./(VREF-BL IT )**4 ) 

7*XJ 

EXAL=EXP(AL*VE ) EQU3 

EXAR=EXP(AL*VREF) EQU3 

TEMP=X J*VREF* ( ( B 1 GB4* ( TE-TREF ) / ( VREF-BL I T ) **4 ) +B6*EXAR* (TE-TREF ) ) EQU3 
DUM*B!GB4*TE/{ VE-BL IT )**4 EQU3 

DUM1«A6*EXAL+B6*TE*EXAL EQU3 

TEMP = TEMP+X J*VE* ( 0 » 5*VE* ( -4 » *DUM/ ( VE-BL I T > + AL*DUMl >+ DUM+DUM1 ) EQU3 

DUM-B I GB4*TE/( VREF-BL IT )**4 EQU3 

DUM1 eA6*EXAR+B6*TE*EXAR EQu3 

TEMP*TEMP— XJ*VRFF* ( 0 • 5* VREF* ( -4 » *DUM/ (VREF-BL I T )+AL*DUM 1 >+ EQU3 

1 DUM+DUM 1 ) EQU3 

TEMP=TEMP+XJ* (B I GB4*TE* ( 1 ./ ( 3. * (VREF-BL I T ) **3 ) -1 ./ (3.* (VE-BL IT) EQU3 
1**3) >+B6*TE* (EXAL-EXAR )/AL ) EQU3 

ANS3*ANS3+TEMP EQU3 

RETURN 
END 


APPENDIX C - Continued 


Program C - Continued 


*DATA A2 = -2 • 162959, A3=4 . 04057E-3 . A4= 1 ,921 072E-4 , A5 = -4 . 48 1 0 49E-6 , 

XKtC* .976E— 2 ,82 = 2,1 35 1 1 4E-3 , B3= 1 . 2828 1 8E-5 «B4= 3 • 008927832-4 , 
B5=9.06231BF-9,XJ=. 185053, 

C2*- 18, 94 1131 , C3=0. 539776, C4=-l , 3023744 1 E-7 , C5 = -4 . 836678E-5 , 

R*0» 1 21 9336 , BL I T = 0 • 1 5E-2,D4= 1 .96802894E-1 1 ,PEE=0.1 026, 

TEE=260,GEE=32. 174, A JA =778,3, A= . 0^8 . B= 1 2000, DELTX* 2 ♦ 

TRFF=820,PREF=l , TZ=409 * 5 « Z4= 1 .90458084E-2 , 

TE I =5*532, El = ,00001 ,F2=. 00001 .MAX 1*30, NUMB = 5 , INT=0. 

I END = 8 , 

PEE* l , 

TEE=532 , 

NO(l 1 = 1 ,2, 3, 4, 5, 6, 7, 8, 9, 

TEI (1 1=8*532. VEI <1 1*8*167,6025. PRAT ( 1 ) *8*, 387, 

U2I (1 1=1500,2000,2500.3000,3500.4000.4500,5000. 

VEE* 1 67,6025 . 

IEND=5. 

XMACHt 1 1=5,091 384,5,939948.6,788512.7,637076.8.485639, 

112 1 (1 1=3000. 3500, 4000 .4500. 5000. 

VE 1(1 1=14.914,12,9277. 1 1 .5794,10.6234,9.91826. 

TE 1(1 1=1338.31 .1583.63. 1 856 . 79 . 2 1 59. 55 . 2492 • 32 , 

PRAT ( 1 1=5*1 . 

VEE= 167,6025. 

VREF =99 . 836 , 

VCll 1=14.28,12,45,1 1.19,10.30.9.63, 

TC (1 1=1343.8.1 589 , 1 862 .2,2165.2. 2498 . 3 , 

A6-5.83882E7. BIGB4*-3.918263E-7,. AL=-6.61 199E2. B6=-9.26392E4, 

PEE= 1 500 » TEE= 1 460 « VFE= 1 .26 1 555E— 1 . 

IFND-1 0. 

XMACH ( 1 1*3.089560. 

XMACH (21* 3.479019,3.858785.4.260383.4.694284.5.174292.5.720277, 

6.191075.6.736399.7.156523. 

U2I (1 1=10*0, 

TE 1(1 1=1447.50,1448.53, 1449.12. 1449.86, 1450.16,1450.52,1451 .19,1451.11. 
1451.28,1451.45. 

VEI <1 1 = • 8249 , 1 .5895.3. 1313.6.5190,14.3189,33.2188.81.6414, 166.849. 
354.500.600.403. 

PRAT ( 1 )*1 • 4 34844 E -2 » “.823559E-3 ,2.3861 69F— 3 .9. 339631 E-4 , 3.475275E-4 . 

1 .221556E-4 .4.019784F-5, 1 . 659644E-5 . 6. 5007«6E-6 . 3. 360228E-6 . 

VC ( 1 )=.70,1.4, 2. 9,6. 0.10. 0.28. 5. 77.0. 160.0,340.0.575.0. 

TC(1 1=1457,9*1456. 

S 

tDATA 

PEE* 1500, TEE* 1260, VEE = 1 • 0 82684 E- 1 , 

I END= 10. 
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APPENDIX C - Continued 


Program C - Concluded 


XMACH (1 >*3.029506, 

XMACH ( 2 ) = 3.489537, 4 ,030964. 4.474 187, 4. 795870, 5. 145822.5.532728. 

5.968807,6.4721 10.7.137680. 

U2I ( 1 )=10*0, 

TEI ( 1 1 = 1244.78, 1 245.«=!5, 1 246.21 . 1 246.63, 1 246.90 ,1247.03,1247.25. 1247.4 3 , 
1247.56. 1 247.7J , 

VEI d 1 = .630 1 .1 .3328,3.3944,7.3098,12.5770.22.2018,40.2732,75.2171 . 

145. 1 13.312.293, 

PRAT ( 1 >=1 .672006E-2.5.892466E-3. 1 .7 1 4 1 49E-3 . 6.400267E-4 ,3 . 2 1 4062E-4 . 

1 .567887E— 4,7.400530E-5,3.361859E-5» 1 .458690E-5.5.446754E-6 , 

VC (1)=. 55. 1.20. 3. 0,6. 9, 9. 5. 18. 5. 36. 0.70. 0.1 38. 0.300.0, 

TCCl 1=1254,2*1253.7*1252. 

S 

*DATA 

PEE= 1500.TEE=l060,VEE=8.976947E-2. 

!END= 10. 

XMACHll 1=3.05391 1 , 

XMACH (21= 3.490262,3.796251 .4. 120834.4.471057,4.856265.5.289648. 

5.791031 .6.392466.7.15091B. 

U21 d 1 = 10*0, 

TE1 ( 1 1=1040.84,1 040.68, 1 040. 77 » 1 040 • 86 » 1 04 1 . 02. 1 04 1 .27.1 041.45,1041.56, 
1041 .76.1 041 .86. 

VEI d }=. 5294.1 .0502,1 .7285,2.9351 ,5.1406.9.2934,17.3738.33.7041,68.2717, 
145.883, 

PRAT ( 1 1=1 ,61 8625E-2 ,6.1 73591 E-3 ,3.1 45939E-3. 1 .559022E-3 » 7. 4891 70E-4 . 
3.471 78E-4 , 1 . 543437E-4 ,6.52045 1 E-5 , 2 . 580923E-5 ,9.349361 E-6 , 

VC < 1 )* *89*1 *6»2*7«4i85»9#0« 16«85t30»0 1 64 • 0 * 140*0* 

TCU 1*1049* 1048*2*1 047*6*1 04 6* 

5 

4DATA 

PEE* 1 500 * TEE=B60 • VEE = 7 . 0 1 4 380E-2 * 

I END* 1 1 * 

XMACHd >*3*002326* 

XMACH ( 2 ) * 3*201 444 * 3. 844078 • 4 .248 1 1 2 ♦ 4 • 706308 * 5* 246 1 22 * 5 • 696985 * 

6 * 072574 *6*417008*6*81 3756 *7*281658* 

U2I d )»11*0* 

TE1 ( 1 >*833* 158*832*331 * 830 • 975 * 830. 705 . 830 * 630 * 830 * 71 0 * 830 * 793* 830 . 844 * 
830.890*830.902*830.965* 

VEI ( 1 ) *.3744 • *5006* 1 .3143*2. 3821 *4.5076 * 8.9480 *14. 9496*22.0079*30*4224 * 
42*6999*61 *0630* 

PRAT d )*1 . 857906E— 2 * 1 . 2 1 5 1 65E-2 * 3 • 1 46945E-3 * 1 . 40 1 97E-3* 5. 931 171 E-4 * 

2 • 350020E— 4 * 1 • l68208E-4*6.B57934E-5*4.366940E-5*2*704147E-5* 1 *616425E-5* 
VC ( 1 ) * . 36 « • 49 * 1 .2*2.1 *4.1 *7.8* 12.0* 19.0*27.0*36*0* 57 * 0 * 

TC(1 >*840*838*836*8*835* 

S 
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APPENDIX C - Continued 


Program D 


NASA-LANGLEY RESEARCH CENTER 


ai 4 

LAR 

01 ? PROGRAM NO 

D 3000 

COMPUTER PROGRAM ABSTRACT 

01 U DATE 

10/29/68 

01 20 TITLE OF PROGRAM (61 CHARACTERS MAXIMUM. 

Hypersonic Nozzle Contours for Tetrafluorome thane 

j PARENT PROGRAM | 

02 M 

CATEGORY 

B 1 

02 15 SITE 

LAR 

02 IJ PROGRAM NO 

DQ700 


02 26 

CATEGORY 

02 27 

LANGUAGE 
NO l 

02 12 

LANGUAGE 
NO 2 

02 37 KEY WORDS B MAXIMUM SEPARATEE) BY COAWAS! 

Tetraf luo rome thane , Nozzle Contour, Displacement Thickness 


F«5R6M 



1 R l 





WHO 

TO CONTACT ABOUT THE PROGRAM 

05 4B STATUS 

05 49 


05 11 CONTACT 


L. R. Boney 


LAR 


05 11 ORGN CODE 


ll.l60 


C5 39 PROJECT NO 


RGL-168 


05 15 
NASA 
CENTER 


1 1 A UNDER DEVELOPMENT 

B OPERATIONAL 
C3 c COMPLETED 


IS NOT FOR 
SHARING 


05 58 REVISION CODE 


05 50 INITIATED 


10/66 


05 54 COMPLETFD 


10/68 


CD A REVISION 
dJ B CANCELLATION 


TIME AND COST FOR DEVELOPMENT 


05 5*> MANMONTH5 


05 64 MACHINE 
HOUR5 


5? SC 61 6? 63 64 65 66 67 68 


05 69 COMPUTER TYPE |05 71 TOTAL COST 
OOLLARSi 


6000 


_» 1 LJ l — 1_ 


74 75 76 77 78 79 80 


ELITE MARGIN 


PICA MARGIN 


CARD NUMBER 


ABSTRACT 


This program calculates the_ turbulent boundary-layer 
disp lac e men t thic knes s al ong a tetra fluorome thane invisci d 


nozzle, contour.. The flow properties at the ed ge of the invis- 


cM^egioa_are supplied Jby Parts, I and H of the program. 


listed .la JIASA TN D-1622 whlch__calculatss. . the JLnvlscid contour 
of a hypersonic nozzle fo r a real ga s. The format of this 


program is identi cal with Part III of the parent. program, 
(NASA TN D-1622) which is restricted to a gas composed of 
diatomic or linear pol yatomic molecules whose vibratio nal 


T 


energy mo des disp lay harmonic oscill a tor charact er istics 

havin g only one fund amen tal fre quenc y. The sha pe parameter 7 1 

the skin- fric tion coefficient , and t he mom entu m thickness are l 

calculated at each, point, along the inviscid^ contour. An l._ 

itera ti ve proce dure d etermines the value of the d i splacement I _ 

thickness which is a dded to the inviscid contour to give a I 

better approximation of _the physic al wal l. Finally a second- j 
order interpolation is employed to locate points on the physi- 
cal wall at desired increments . |. 

■ - -- -t 

— — — - T 

- - 1 

] 

-t 

-.1 


4 



APPENDIX C - Continued 
Program D - Continued 


SIBFTC 700] L I FT , REF , DECK , M94 

C P-706.1 0010 
CLRR 0020 
C PART III COMPUTES BOUNDARY LAYER USING 0030 
C WALL POINTS PUNCHED BY PART II 0040 


EXTERNAL FOFX 
D I MENS I ONTHB I(IO).GIUO) 

COMMONTHBI .GI .NOTH 

D I MENS ION FS (2 ) , SUM ( 2 ) , ANS < 2 ) . X (2000 ) • Y (2000 I »W<2000 ) .RRT <2000 ) ♦ 0050 

1 THF (2000 ) ■» T I (2000 ) » SVY ( 2000 ) . SVDELT (2000 ) .DELS (2000 ) « YY (2000 ) 0060 


COMMON 

OPT 

9 T I 

9 PRI 

9 TO 

9 THB 

9 QSUM 

0070 

COMMON 

QW 

• EN 

9 ERR 

9 UI 

9 HW 

9 R 

0080 

COMMON 

HT 

• CONI 

9 C0N2 

9 CON 

9 T 

9 H 

0090 

COMMON 

0 

♦ FR 

9 FPR 

f H3 

9 HWS 

9 K 

0100 

COMMON 

RHOI 

• Q I 

9 X 

9 Y 

9 W 

9 RRT 

01 10 

COMMON 

THF 

• svv 

9 SVDELT 

9 DELS 

9 YY 

9 RHOT 

0120 

COMMON 

VE 

• PI 

9 Cl 

9 D 1 

9 El 

9 THI 

0130 

COMMCN 

SC 

* OMEGA 

9 TW 

9 XM 

9 TLTG 

9 JL I M 

0140 

COMMON 

DEBUG 

9 N 

9 L 

9 ACASE 

9 ALPHA 

, BETA 

0150 

COMMON 

GAM 

9 TAU 

9 PR 

9 TT1G 

9 XI 

9 Y 1 

0160 

COMMON 

EMi 

9 W1 

9 RRT 1 

9 THF 1 

9 Til 

9 TLTI 

0170 

COMMON 

TT1I 

t THETA 

9 DXB 

9 DU 

9 DR 

9 DY 

0180 

COMMON 

PI 

9 EMU 

9 RETH 

9 TEMP 

9 FT 

9 FTPR 

0190 

COMMON 

HI 

9 DSDEL 

9 THDEL 

t TAW 

9 . alp 

9 CW 

0200 

COMMON 

C2 

9 C3 

9 C4 

9 FTT 

♦ FTTPR 

9 HZ 

0210 

COMMON 

CF2 

9 DSTH 

9 DTDX 

9 YPD 

9 DEL 

9 DP 

0220 

COMMON CP 
DIMENSION FS (2 ) 

.BP , A P .HI 

.SUM (2 ) , ANS(2> .X (2000 ) « Y(2000) ,W(2000 ) .RRT <2000 ) , 

0230 

0240 


C 1 THF (2000 ) ,TI (2^00) , SVY (2000 ), SVDELT ( 2000 ) .DELS (2000 I .YY(2000 ) 0250 

C COMMON OPT 0260 

C COMMONTI »RRI .TO.THB.OSUM.QW.EN.ERR.UI < HW . R «HT .CONI ,CON2, 0270 

C 1 CON » T » H , (5 . FR . FPR • H3 , HWS . K .RHO I .01 . X . Y , W . RRT . THF . SVY . SVDELT . DELS . YY 0280 

C 2 .RHOT.VE.Rl ,C1 .D1 .El .THI ,GC. 0290 

C 3 OMEGA, TW« XM , TLTC. JL I M , DEBUG ,N «L * ACASE .ALPHA . 0300 

C a PET A « GAM « TAU .PR .TT1G.X1 «Y1 »EM1 ,W1 , RRT1 , THF1 *T I I ♦ 0310 

C E TLTI ,TT1 1 .THETA, OXB, DU. 0320 

C 6DR.DY, PI, EMU. RE TH» TEMP .FT , FTPR , HI , DSDEL. 0330 

C 7THDEL , T Aw , ALP .CW.C2.C3.C4, FTT . FTTPR « H2 » CF2 » DSTH. 0340 

C 8CTDX.YPD.DEL 0350 

C 9 , DP, CP. BP. AP, HI 0360 

1 WRITE (6,25) 0370 

2E FORMA T ( 1 H 6X5HACASE 1 3X4HRH0T 1 4X2HVE 1 5X2HR 1 1 5X2HC1 1 5X2HD1 0380 

1 15X2HE1 /7X3HTH I 1 5X2HGC I 5X1 HRl 5X3HTHB 0390 

21 3X5H0MEGA1 4X2HTW1 4X3HERR/8X2HXM1 4X4HTLTG 0400 

31 2X5HALPHA 1 3X4HBET A 1 3X3HGAM 1 4X3HT AU1 5X2 HPR/ 04 1 0 

47X4HTT1 Gl 4X2HT0 1 5X2HAP 1 5X2HBP1 5X2HCP 1 5X2HDP 1 5X2HHT ) 0420 

]01 FORMAT ( 1 H 7E17.8) 0430 

100 FORMAT ( 1 H06EI 7 . 8/6E 17.8) 0440 


! I i 
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APPENDIX C - Continued 


Program D - Continued 


C ALLF I NP ( 4 ♦ THB I *G I « NOTH *NOTHM ) 

CALLF I NP < 39 * RHOT * VE *R 1 • Cl • D1 t E 1 • TH 1 * GC « R « THB ♦ OMEGA ♦ 0450 

1 TW * ERR * XM * TLTG * JL I M * DEBUG iN»L» ACASE • ALPHA iBETA, GAM « 0460 

2TAU * PR • TT 1 G * TO • AP * BP * CP • DP * HT * APP * BPP * CPP * YSTAR t 0470 

30PT*DIFN*D?FD ) 04 80 

T HB=THB I (NOTHM ) 

WRITE (6* 1 03 )ACASE*RHOT* VE* Rl , C t • D t * E I * TH I ♦ GC * R * THB * OMEGA ♦ TW * ERR * X 0490 

1 M*TLTG, ALPHA .BETA .GAM*TAU«PR*TT1G * TO * AP « B 0500 

2P*CP*DP*HT 0510 

26 FORMAT ( 1 H 8X 1 HX 1 6X 1 HY 1 4X5HTHET A 1 3X4HCF/2 0520 

1 1 4X1 HN1 4X6HTH/DEL/5X7HY+DELST 1 I X5HDELST t 3X3HDEL 0530 

21 4X4HRETH1 3X4HTLTI 1 0X9HDELST/DEL > 0540 

J=JLIM 0550 

C SUBROUTINE BIRD READS BINARY CARDS CONTAINING Xt Y«M * Wi RHO/RHOT • 0560 

C THETA « T 0570 

C LIMIT OF 2000 CARDS 0580 

C READ IN WALL POINTS FX I T TO THROAT 0590 

C MULTIPLY X AND Y BY XM AND STORE THROAT TOEXIT 0600 

lOOO CALLBIRD* XI *TI I ) 0610 

X(J)=X1*XM 0620 

Y(J)=Y1*XM 0630 

W(J)=W1 0640 

RRT ( J ) -RRT 1 0650 

THF ( J ) =THF 1 0660 

T I ( J ) =T I 1 0670 

SVY(J)sY(J) 0680 

50 SVDELT ( J ) s 0 0690 

J= J- 1 0700 

I F ( J ) 8 1 *81 *1000 0710 

P 1 TLTI=TLTG 0720 

TT 1 1 =TT 1 G 0730 

THETA=TH! 0740 

WRITE (6*26) 0750 

CALL SL ITE ( 1 ) 0760 

2 DOQO J* 1 * JL I M 0770 

K = J 0780 

RHOI sRRT { J ) *RH0T 0790 

U I = W ( J ) *VE 0800 

PI=R1*TI (J)*RH0I 0010 

EMU= ( C 1 * ( T I ( J) **1 *5) )/( ( T I (J) **D1 ) + E 1 ) 0820 

C EQUATION BIO 0830 

IF (PI -20000, 0 ) 203 * 203 * 204 0840 

203 01=1.0 0850 

GOT04 0060 

204 OI=PI*(PI* (PI* ( TAU ) +GAM ) +BET A ) + ALPHA 0870 

4 I F ( TH0/TT 1 1-24.0) 120*120*101 0880 

1 0 1 FT = VE*VF-QI* (8**TT1 1 *R ) 

ETPP = -0 r *8.*R 

GOTO 102 0910 

120 S E M = 0 
REM®0 
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APPENDIX C - Continued 


Program D - Continued 


DO 1 20 1 ! = ! ,W6TW 


tem^expcthbi u >/tti i ) 

SEM*GI <J >*THBM I )/<TFM-I • V+SE*. 

1201 QEM*GI ( I ) * T HB I ( i )*TH0I ( I ) *TEm/ ( TEM- 1 • )**2+REM 

ft*ve*ve-gi* (8 **tti i#r+2.*r*sem ) 

FTPR=-01*(8**R+2.*R4REM/(TT1 1*TTU > ) 

102 H 1 s-FT/FTPR 
TT1 1 =TT 1 1+H1 

T F ( ABS (HI /TT 1 1 >-ERR)5*5*4 
C EQUATION 816 

5 IF (CPP >400 *401 * 400 

401 TW=TT 1 1 
GOTO402 

400 TW= ( CPP— APP* (5VV ( J ) /YSTAR )** 1 • 0 )/ ( 1 *0+BPP* ( SVY ( J l/YSTAR ) ** 1 #8 ) 

C FQUATION 81 1 

402 JF(Pi-20000*0)205, 205*206 

205 0W*1«0 
G0TO207 

206 QW*P I * ( P I * (P I * ( DP ) +CP ) +BP > + AP 

207 IF (THB/TW-24 *0)104,1041105 

105 HWS*4**TW/T0 
GOTO 1 06 

1 04 SEM*0 

DO 1 04 1 I * 1 , NOTH 
TEM=EXP (THB 1(1) /TW ) 

1041 SEM«GI f I )*TH8I C J >/(TEM-l • )+SEM 
HWS*4**TW/T0+SEM/T0 

106 I F (THB/T ! (J >-24.0)107# 107,108 
108 HI=4.*T1 (J)/TO 

G0T01 09 
10T 8EM=0 

0010711*1 , NOTH 
TEM=EXP < THB I (I ) /T I ( J ) ) 

1071 SEM=G1 ( I )*TH01 ( I ) / ( TEM— 1 • J+SEM 
H1«4**T1 ( J )/TO+SEM/TO 
\09 OSUM* (QI-QW)/<HI-HWS) 

I F ( THB/TW-24 * 0 ) 1 10* 1 10* 1 1 1 

1 1 1 HW*HWS*QW 
GOT 0 1 1 2 

I 10 HW*HWS*OW 

112 RET*(RHOI*UI )/ ( EMU* 1 2 • 0 ) 

507 RFTHeRFT*THFTA 

C EQUATION B5 

EN=1 # 77* • 43429448* ALOG (RETH >-.38-200 *0/RETH 
CALL SLITETM *<000FX) 

GO T0C51 1 *502) *K000FX 
502 IF (ABS (EN-ENP)-DIFN >500*500 *501 
511 CALL 5LITE (1) 

501 RR 1**125 

C SUBROUTINE GAUSS INTEGRATION 5 INTERVALS AND 32 POINTS PERINTERVAL 

A«0. 


0960 
0970 
0980 
0990 
1000 
1010 
1020 
1030 
1040 
1050 
1 060 
1070 
1080 
1090 

1110 


1 1 30 
1 1 50 


1 1 70 
1 180 

1200 

1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 
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no on 


APPENDIX C - Continued 
Program D — Continued 


B* 1 . 

N0FNx2 

CALL GAUSS (N<L< A 1 B 1 NOFN * S * FS ( 1 > . SUM ( 1 ) . ANS ( 1 ) .FOFX ) 1320 

C DELTA STAR/DELTA f33 ° 

DSDEL* 1 » O-ANS ( 1 ) 1340 

C THETA/DELTA 1350 

THDEL=1 .0— DSDEL-ANS (2 ) 1360 

TAW= (PR**.33333333 >* <TT 1 1-T I ( J ))+T I ( J ) 1370 

ALP= ( 1 .O+OMEGA >/ (EN+ 1 . 0 > 1380 

CW=TW/T I ( J ) 1 390 

C2= (TW-TAW )/Tt ( J) 1400 

C3 = ( 20 . 0*EN*THDEL/RETH ) ** ( 1 .0/(EN+l . 0> ) 1410 

C4=(TAW-TI (JM/TI (J) 1420 

7 FTT*CW— C2*C3*(TLTI**ALO)-C4*C3*C3*<TLTI**<2.0*ALP) l-TLTI 1430 

FTTPRx- ALP*C2*C3* ( TLT 1 ** < ALP— 1.0) 1-2 »0*ALP*C4*C3*C3* _ 1440 

1 (TLTI**(2.0*ALP-1 .0) ) - 1 . 0 1450 

H2=— FTT/FTTPR 1460 

EQUATION B 1 4 1470 

TL/TI 1480 

TLT I =TLT I+H2 1490 

IF(ABS(H2/TLTI )-ERR)6.6.7 1500 

EQUATION B 1 3 1510 

CF/2 1520 

6 CF2*<1.0/( 20 .0*EN ) 1 *C3*C3* ( ( TLT I >** ( (1.0+2. 0*OMEGA-EN ) 1530 

I / ( FN+ 1 .01)) 1540 

C DELTA STAR/THETA 1550 

DSTH=DSDEl/THOFL 156 P 

500 DELS < J ) =DSTH*THET A 1570 

YPD=SVY ( J l+DELS ( J > 1580 

DEL “DELS ( J ) /DSOEL 1590 

WRITE ( 6 . 1 00 ) X ( J ) . Y ( J ) . THETA , CF2 . EN , THOEL . YPD . DELS ( J ) . DEL . RETH , TLT 1600 

1I.DSDEL 1610 

YY(J)*Y(J) 1620 

Y ( J ) * YPD 1630 

IF (J- 1 ) 504 .505. 504 1640 

505 DELSP=DELS ( J ) 1650 

G0T05Q9 1660 

504 CXB= (X ( J ) -X ( J- 1 ) > /COS ( THF < J- 1 )) 1670 

DU = VE* (W ( J )-W ( J-l ) ) 1680 

DR=RH0T * (RRT ( J ) -RRT ( J— 1 ) ) 1690 

DYxy ( J )-Y ( J- l ) 1700 

C D THETA/D X 1710 

DTDX=CF2P-THETAP*( ( ( 2 • O+DSTHP ) *DU ) / ( U I P*DX8 >+DR/ (RHOIP*DXB ) 1720 

t +DY/(Y (J— 1 ) *DXB > > 1730 

THETA = THET AP+(DTDX>*DXB 1 74 ° 

CALL SLFTETd .KOOOFX) I 750 

GO T0(508. 506 ) .KOOOFX 1760 

506 IF (ABS ( (DELS (J ) -DELSP )/SVY ( J )) -DIFD > 509.509.508 1770 

508 DELSP=DELS ( J ) 1780 

ENP=FN 1790 

G0T0B03 1800 
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APPENDIX C - Continued 


Program D — Continued 


50 <? )=V< J+n+DELS< J> 16i0 

CF2P=CF2 1020 

THETAP=THETA 1830 

OSTHP=DSTH 1840 

UIP=UI 1850 

RHOIPaRHOI 1860 

51 n CALL SLITE (1) 1070 

80 CONTINUE" 1080 

1001 X I NT = 0 1S90 

WRITF (6,514) 1900 

C USE SECOND ORDER INTERPOLATION TO EIND GIVEN POINTS 1910 

C DX* • 1 0 IF Y LESS THAN 1. 1920 

C DX=.25 IF Y GREATER THAN 1. AND LESS THAN5. 1930 

C DX«.50 IF Y GREATER THAN 5. 1940 

1005 CALL FTLUP (X INT, YINT.2 , JLI M,X ( 1 > » YY( 1 ) > 1950 

WRITE ( 6 » 51 3 )X I NT , Y I NT I960 

I F (Y INT— YY ( JL I M ) ) 1006* 1 * 1 

1006 IFCYINT-1 .0 ) 1002, 1009, 1 009 1980 

1002 X INT=X I NT+ • 1 1990 

GO TO 1005 2000 

1009 IF(YINT-5.0) 1003, 1004, 1004 2010 

1003 XI NT*X INT+ • 25 2020 

GO TO 1005 2030 

1004 X I NT *X I NT + • 5 2040 

GO TO 1005 2050 

513 FORMAT < 1 H 2F12.4) 2060 

514 FORMAT ( 1 H 6X 1 HX 1 1 X 1 H Y ) 2070 

END 2080 

$ IBFTC FOFX LI ST, REE, DECK, M94 

SUBROUT I NEFOFX(S,FS> 0010 

D I MENS I ONTHB 1(10), GI (10) 

COMMONTHBI *GI .NOTH 

D I MENS ION ES (2 ) * SUM ( 2 ) , ANS ( 2 ) , X <2000 ) , Y ( 2000 ) , W <2000 ) ,RRT (2000 ) , 0020 


1 THF (2000 ) ,T I (2000 ), SVY ( 2000 ) , SVDELT (2000 ) .DELS (2000 ) . YY (2000 ) 0030 


COMMON 

OPT 

9 T I 

9 PPI 

9 TO 

9 

THB 

9 OSUM 

0040 

COMMON 

GW 

• EN 

9 ERR 

9 UI 

9 

HW 

9 R 

0050 

COMMON 

HT 

t CONI 

9 C0N2 

9 CON 

9 

T 

9 H 

0060 

COMMON 

0 

t FP 

9 FPR 

9 H3 

9 

HWS 

9 K 

0070 

COMMON 

RH01 

• 01 

9 X 

9 Y 

9 

W 

9 RRT 

0080 

COMMON 

THF 

9 SVY 

9 SVDELT 

9 DELS 

9 

YY 

9 RHOT 

0120 

COMMON 

VE 

• R 1 

9 Cl 

9 D 1 

9 

El 

9 thi 

0130 

COMMON 

GC 

9 OMEGA 

9 TW 

9 XM 

9 

TLTG 

9 JL I M 

0140 

COMMON 

DEBUG 

t N 

9 L 

t ACASE 

9 

ALPHA 

9 BETA 

0150 

COMMON 

GAM 

♦ TAU 

, PR 

t TT1G 

9 

XI 

9 Y 1 

01 60 

COMMON 

EMI 

9 W1 

9 RRT1 

9 THF 1 

9 

T I 1 

9 tlti 

0170 

COMMON 

TT1 1 

9 THETA 

9 DXB 

9 DU 

9 

DR 

9 DY 

0180 

COMMON 

PI 

9 EMU 

9 RETH 

9 TEMP 

9 

FT 

9 FTPR 

0190 

COMMON 

HI 

t DSDEL 

9 THDEL 

9 TAW 

9 

ALP 

9 CW 

0200 

COMMON 

C2 

9 C3 

9 C4 

9 FTT 

9 

FTTPR 

9 H2 

0210 

COMMON 

CF2 

9 DSTH 

9 DTDX 

9 YPD 

9 

DEL 

9 DP 

0220 

COMMON 

CP 

t BP 

9 AP 

9 HI 
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APPENDIX C - Continued 
Program D - Continued 


C D I WENS ION FS (2 ) *SUM (2 ) * ANS < 2 ) * X (2000 ) • Y (2000 ) « W (2000 ) «RRT <2000 ) t 

C 1 THF ( 2000 )• T I ( 2000 ) * SV Y ( 2000 ) * SVDELT ( 2000 ) « DELS (2000 ) ♦ YY ( 2000 > 

C COMMON OPT 

C COMMONT I • RR I • TO t THB «QSUM «QW «EN ♦ ERR • U I « Hw • R *HT « CONI « C0N2 * ^ 

C 1 CON tT*HUQ«FR«FPR*H3«HWS«<* RHO I *QI «X.Y«W*RRT • THF « SVY « SVDELT « DELS « YY 

CONI « C <UI*U! )/ (R*T0 ) ) 

C0N2*S** ( 1 • O/EN ) 

CON = C ONI * (HW* ( 1 • 0-CON2 ) /CONI +HT*C0N2/C0N 1 - • 5*C0N2*C0N2 ) 

C ITERATE FOR RHO/RHOI 

11 T*TI(K)/RRI 

C IF EXPONENTIAL GREATER THAN 24* OMIT TERM 

C EQUATION B9 

IF (THB/T-24* 0) 1 1 3 • 1 1 3 • 1 1 4 

1 1 4 H = 4 • *T/TO 
GO TO 1 1*52 

1 1 3 SEM=0 

D01 1 31 I = 1 * NOTH 
TEM*EXP(THBI <1 )/T) 

1131 SEM*GI ( I )*THBI ( I )/(TEM-l • )+SEM 
H*4 *#T/T0 + 5EM/T0 
1 1 52 IF (OPT > 1 1 5* 1 t 50 1 1 15 

1150 0*01 

GO TO 1151 

115 O* (OSUM )* (H-HWS )+QW 

'c IF EXPONENTIAL GREATER THAN 24* OMIT TERM 

1151 IF <THB*RPI /T I (K )-24.0 )116#U6tll7 
117 FR=Q*4**TI (K)/(T0*RRI ) -CON 

FPR * — 4 • *0*T I (<)/<TO*RRI*RRI j 
GO TO lie 

1 1 6 SEM*0 
RFM*0 

HOI 161 1*1 * NOTH 

TFM=EXP ( THB Id) *RR I /T I ( K ) ) 

SEM«G I ( I )*THBI ( I l/CTEM-l • )+SEM 
1161 REM*G I ( I ) *THB I ( I ) *THB I ( I ) *TEM/ 

1 (TEM-1 • )**2+REM 

FR*Q* ( 4 * *T I (K)/(TO*RRI ) +SEM/T0 ) -CON 
FPR=0* ( -4 **T I (K )/ <TO*RR I*RR I ) —REM/ 

1 <TO#TI (K) ) ) 

1 10 H3*-FR/FPR 
C EQUATION B 1 2 

RQIN=RRI+H3 
IF (RPIN ) l 25* 125, 126 

125 RR I =RP I * • 5 
GO TO 1 1 

126 RPI=RRIN 

I F ( ABS ( H3/RR I ) -ERR ) 10 ' 1 0 * 1 1 
10 FS( 1 )«RRI* (S** ( 1 *0/EN) ) 

40 FS<2)*RRI*(S**(2.0/EN) ) 

HOLD* S 
HOLD?»FS < 1 ) 


0100 

0110 

0120 

0130 

0140 

0150 

0160 

0170 

0180 

0190 

0200 

0210 

0220 

0240 


0260 

0270 

0280 

0290 

0300 

0310 


0340 


0390 

0400 

0410 

0420 

0430 

0440 

0450 

0460 

0470 

0400 

0490 

0500 
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APPENDIX C - Concluded 


Program D - Concluded 


u 

HOLD3=F5<2 ) 

0510 

100 

FOPMAT ( 1H 7E17.8 ) 

0520 


PFTUPN 

0530 

SDATA 

FND 

0540 


J1 = D*N1 =2341 • N 1 127.N3322.N1653.F 
J2=0*N1 =1 .N2.N3.N4.F 
J3-4B35 ♦ J4=2B35$ 

N1 = I 2*451 »N2 = 4 1 93.95 tN3= 1 7. 55 * N4 = 7 . 4 866E-7 t N5= 1 .0tN6=251 • 1 7.N7=.0001 t 
N6»l . »N9=565. %NU = *6*N12 = 66 0. t N 1 3= . 00005 * N I 4 = • 205?7 « N 1 5= 1 • t N2 1 * J * *N22=0t 
N23=0*N24*0*N25=.76.N26=1460. «N2 7=455. .N28=l • t N29 =0 * N30 =0 « N3 1 =0 « 

N32 = 34 . 2 ♦ N33«0 . N34 = 0 • N35= 1374#«N36=1 • . N37= 0 « N38= • 0 1 «N39*.001 ♦ 
N10=0.N16=388B35 tNl 7 = 0 »N1 0 = 32B35«N1 9=5B35 ♦ N20 = 5 * S 
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T 0 ,°K 



(a) p Q = 1034 n/ot 2 . 

Figure 11.- Variation of free-stream specific-heat ratio with Mach number. 
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(C) Pg = 1724 N/W. 
Figure 11.- Concluded. 
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Figure 12.- Variation of ratio of free-stream static temperature to stagnation temperature with Mach number. (Independent of stagnation 

pressure from 1034 to 1724 N/cm^J 
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Mach number. (Independent of stagnation pressure from 1034 to 1724 N/cm2.| 
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Figure 17,- Variation of ratio of total pressure behind a normal shock to free-stream total pressure with Mach number. (Independent of 

stagnation pressure from 1034 to 1724 N/cm2 ) 
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Figure 19.- Variation of vJv, with Mach number. (Independent of stagnation pressure from 1034 to 1724 N/cm2.) 
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(b) p Q = 1378 N/cm2. 
Figure 20. ~ Continued. 
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Figure 2i.- Ballistic- range shock density ratio. 
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Langley Unitary Plan wind tunnel- 
2 < M < 5 


Langley 22- inch 
■ helium tunnel 
M = 20 
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Figure 25.- Density-ratio simulation capability of several facilities. 



Figure 26.- Variation of shock density ratio with shock-wave inclination angle. 
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Figure 27.- Schematic diagram of CF 4 tunnel. 


Pilot CFj tunnel 



Figure 28.- Variation of area ratio with Mach number. 




(a) Pilot CF 4 tunnel. Mi = 6.2; P 2 /P 1 ~ 12.1. 



(b) Langley Mach 8 variable-density tunnel. M^ = 8; P^i - 5.56. 


Figure 29.- Shock shape for 120° cone. L-69-1323 



Figure 30.- Drag coefficient of 120° cone. 





Figure 32.- Drag coefficient of tension shell. 




(a) M - 6.14 in air. 



(b) M = 6.61 in CF 4 . 

Figure 33.- Shock shapes for tension shell. L-69-1324 
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Figure 35.- Sketch of CF 4 nozzle. 




